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Summary 


In this article a novel method for the design of a bass-reflex loudspeaker enclosure is 
introduced. It is shown that simplification and transformation of the electrical analogue 
circuit makes the latter susceptible to the discipline of network synthesis. By this method 
the design constants of the system can be found in a straightforward manner. The resulting 


~ data, when grouped into general forms, show the implications for actual construction of the 


system very clearly. 

Important properties like bandwith (lower cut-off frequency), size and total damping 
required are shown graphically. Assets as well as basic limitations, both of the method and 
the realization, are discussed. In the appendices some details of the synthesis procedures and 
of motional feedback (to increase the damping) are outlined. 


Sommaire 


Dans cet article on expose une nouvelle méthode de détermination des données de con- 
struction d’un haut parleur dit «bass-reflex». On montre que la simplification et la trans- 
formation du circuit électrique équivalent lui permettent d’étre traité par les méthodes 
applicables aux réseaux électriques. A l’aide de cette méthode les constantes caractéristi- 
ques du systéme peuvent étre trouvées d’une facon directe. Les résultats obtenus groupés 
dans des formes générales montrent trés clairement les applications réelles 4 la construction 
du systéme. 

D’importantes propriétés, telles que la largeur de la bande (fréquence du cut-off infé- 
rieur) ainsi que l’amortissement total, sont représentées graphiquement. On discute l’effet 
des restrictions de base et de fait imposées par la méthode ou par la réalisation. Dans 
Yappendice on expose en détail les procédés de calcul ainsi que la réaction motionelle (afin 
d’accroitre l’amortissement) . 


Zusammenfassung 


In dieser Arbeit wird eine neue Methode zur Entwicklung eines BaBreflexsystems ein- 
gefiihrt. Es wird gezeigt, da nach einer bestimmten Vereinfachung und einer Transforma- 
tion das elektrische Ersatzschaltbild vierpol-theoretisch behandelt werden kann. Die Groen 
der elektrischen und akustischen Elemente kénnen dann auf direkte Weise bestimmt wer- 
den. Wenn wir diese Daten zu allgemeinen Formen kombinieren, werden Vor- und Nachteile 
des Bafireflexsystems besonders deutlich. 

Wichtige Daten sind die Bandbreite (untere Frequenzgrenze), die GroBe des Luftvolu- 
mens und die gesamte Dampfung; diese werden in Kurven dargestellt. Die Anwendungs- 
grenzen dieser Methode und die physikalischen Beschrankungen des Systems werden kurz 
diskutiert. SchlieBlich werden Einzelheiten der Synthesemethode und der zur Erreichung 
gréBerer Dampfung benutzbaren dynamischen Gegenkopplung in drei Anhangen angegeben. 


1. Introduction ii. Approximation. 


The response function is 
adapted to the response requirements. The coeffi- 
cients are now prescribed to have specific values. 


design of a bass-reflex loudspeaker enclosure should 
theoretically lead to the following three steps: 


i. Analysis. The electrical analogue is analyzed 
as to its response for sinusoidal excitation. The 
response will be the quotient of two polynomials in 
the complex frequency s. The coefficients in the 
response function are combinations of circuit ele- 
ments. 


iii. Synthesis. The design is completed by deter- 
mining the circuit elements in such a way that the 
desired coefficients in the response function are 
obtained. 

In the case of the bass-reflex system this proce- 
dure cannot be usually closely adherent to. This may 
partly be due to an unsuitable starting point being 
chosen. Another reason is that the final step neces- 
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sitates the solution of simultaneous equations of a 
high degree, which is a difficult problem indeed. In 
conjuction with this there exists the problem of 
constraints imposed upon the solution by reasons 
of physical realizability. It is therefore not surpris- 
ing to find that the usual procedure of bass-reflex 
design is based on a rather weak theoretical basis. 
Design nomographs are consulted to make the reso- 
nance frequency of the enclosure equal to that of the 
loudspeaker. The cabinet damping is adjusted to 
equalize the peaks in the voice-coil electrical impe- 
dance. In this method no attempt is made to pre- 
dict the acoustic sound pressure response. At most 
the behaviour of the mechanical impedance pre- 
sented to the diaphragm is discussed. In exceptional 
cases a comprehensive analysis is given [1], [2]. 
The present article deals with the use of network 
synthesis methods in the design of the bass-reflex 
enclosure. With such a method the system can be de- 
signed with relatively little computational effort. 
‘The success of any method in this field depends, of 
course, on the starting point of the analysis. We 
will utilize here the most simplified network that 
still represents bass-reflex action, i.e. efficient sound 
radiation at maximal diaphragm impedance in a 
certain frequency region. This report will mainly be 
concerned with the basic theory so that questions 
relating to physical realization will not be fully 
discussed. 


2. Analysis 


A fairly comprehensive electrical analogue of a 
bass-reflex system is given in Fig. 1. In this circuit 


Diaphragm Port 
Radiation Radiation 
1,=Up Map Map Ih=Up 
Rp —_> 
A\\—-~ 000 


Co Mp 


1® Rem 


Enclosure 


Fig. 1. Electrical analogue of bass-reflex loudspeaker 
system. 


currents represent volume velocities and voltages 
represent sound pressures. Electromagnetic and 
frictional damping are represented by resistors Rey 
and Rp, respectively. The quantities Mp und Cp 
are the mass and the stiffness of the diaphragm. The 
impedance Z,p (represented by a parallel combina- 
tion of a resistance Rap and an inductance Map) 
is the radiation impedance of the diaphragm open- 
ing. The radiation impedance Zap of the port (the 
auxiliary aperture) is represented in a similar way. 
The reactance of the enclosed air volume is given 
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by the capacitance Cy. Finally, the effective mass 
of air contained in the duct, if present, is shown. 
The mutual inductance m represents the near-field 
interaction between port and diaphragm [3]. 

The sound pressure at a large distance from the 
enclosure consists of two parts. One originates from 
the diaphragm, the second contribution is due to 
radiation of the port. If we regard both sources as 
point sources, we may combine them into a single 
source as long asthe wavelength of the sound is 
sufficiently large. The volume velocity of the com- 
bined sources is the sum of the volume velocities 
Up and Up as represented by the currents J, and J, 
in the circuit of Fig. 1. Signs taken into account, 
the effective source strength is then proportional 
to the current /; [3]. This means that the acoustic 
sound pressure response, in which we are ultimately 
interested, can be found from the electrical analogue 
as well as the more directly evident quantities. 
Here is need to study more than the total mechanical 
impedance to the diaphragm as is usually dure. 

For a successful analysis the electric, circuit has 
to be simplified considerably. Essentially it contains 
a series resonant circuit (the loudspeaker unit), 
connected in series with a parallel circuit (the en- 
closure). At lower frequencies the radiation resist- 
ances do not affect the response much. Leaving them 
out altogether, and combining the remaining ele- 
ments, one obtains the circuit of Fig. 2a. The 
mutual inductance m is shown in series with capa- 
citor D, due allowance is made‘in the other induct- 
ances to ensure full equivalence. 

All damping is seen to be lumped into one re- 
sistor, the enclosure proper is no longer damped 
either internally or by acoustic radiation. This 
will have some peculiar consequent. At the fre- 
quency of resonance of the enclosure, the impedance 
presented to the diaphragm will be infinite. Yet 
the current /, will not be zero so that the response 
may well be regular in this frequency region. 

The acoustic response is governed by the transfer 
admittance 


Y,3(s) =I5/E. 


In fact the sound pressure is proportional to this 
quantity and to frequency [3]. The acoustic sound 
pressure response function is easily found to be 


ees (1) 


l+as+bs?+cs?+d5* 


where a, b, c and d are determined by the following 
network constants: 


a=RC, 
b=(M+N—2m)C+ND, 
c=RCDN, (2) 


d=CD(MN—m?). 


nn 2 
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3. Approximation 


The form found for the response function is suf- 
ficiently simple to render the approximation prob- 
lem easy. For simplicity the response will be re- 
quired to be of the Burrerwortu type. The response 
has then to satisfy the relation 

| 2G) P_ 1 
[ok | 1+ (ofer)® 


When the cut-off radial frequency is normalized to 


unity, the function z(s) becomes 


k a: (3) 
1+ 2.613 5+ 3.414 5s? + 2.613 s?+ s4 


g(s) = 


_It is only slightly more difficult to prescribe a re- 


sponse of the Cuesysuev form. The conversion from 
the w?-plane back to the s-plane causes some diffi- 
culties, however [4]. For simplicity the computa- 
tions will only be carried out for the Burrerwortu 
case. Results for Cuesysuey-type responses will be 
given also. 


4. Synthesis 


To bring eq. (1) and eq. (3) into agreement, a 
set of equations like eq. (2) has to be solved for the 
unknowns M, N, m, C, D, and R. The equations 
are of the fourth degree. Since there are four equa- 
tions two parameters can be chosen freely. One of 
these degrees of freedom is irrelevant because it 
corresponds to the over-all impedance level. The 
second degree of freedom gives a means to obtain 
the same response with various sets of parameters. 
This freedom is, of course, very useful. 

The process of solving the set of equations (2) 
with a, b, etc. given as in equation (3), is quite com- 
plicated, to say the least [1]. We will show pres- 
ently how the network synthesis approach elimina- 
tes this problem. In a certain sense the synthesis 
method is the inverse of the solution as described 
above. Instead of adjusting the network constants 
to the specifications, one tries to construct a whole 
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Fig. 2. (a) Simplified version of Fig. 1. 
(b) Circuit transformed according to recipro- 
city theorem. 
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network that has the required response. If the re- 
sulting network happens to be of the correct geo- 
metry, the problem is solved, because network syn- 
thesis not only yields a network configuration but 
also the element values. In the application of net- 
work synthesis one can in general exert a certain 
amount of control over the initial steps in the pro- 
cedure. These steps can be aimed at the desired 
network thus increasing the chance of success. 

The network of Fig. 2 a is susceptible to the syn- 
thesis approach. The quantity in which we are inter- 
ested is the transfer admittance Y,,. It is given by 
the following equation: 

8 


4 ak : 3 
138) =k 3 54 3414 00613 Pat 


(4) 
In fact Y,3 implies that we apply a voltage source at 
the input and measure the output current J, . Follow- 
ing the principle of reciprocity the positions of 
voltage source and current meter can be inter- 
changed; the transfer admittance of the new con- 
figuration is then the same as before. The result of 
this procedure shows Fig. 2b. This network has a 
very special structure. It is a lossless coupling 
network between a voltage source and a load re- 
sistor. Such a configuration is well studied in net- 
work theory. Its synthesis is known as the Daruinc- 
ton problem [5]. 

Appendix J describes the synthesis of a lossless 
coupling network that in conjunction with a load 
resistor has the transfer admittance given by eq. (4). 
The result is the network of Fig. 3. The output re- 
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Fig. 3. (Upper part) Result of primary synthesis proce- 
dure. 


(Lower part) Desired network configuration. 


sistor has been normalized to one ohm, which pro- 
cedure eliminates the irrelevant degree of freedom. 
The resulting network has, however, not yet the 
required geometry. It contains two inductances and 
two capacitances, corresponding to the minimum 
number of reactance elements necessary to produce 
four poles in the admittance functions. The desired 
network of Fig. 2b has three inductances and two 
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Fig. 4. Final result of synthesis. 


capacitances. One element is superfluous, which 
feature corresponds with one free parameter. 

This difficulty can be removed by applying a so- 
called internal impedance transformation [6]. The 
procedure is described in Appendix II, the result 
shows Fig. 4. The parameter A can be chosen freely 
within the range 


0.585 ara sti 


At this point the circuit configuration coincides 
with that of Fig. 2b. Without solving non-linear 
equations the network elements are found directly. 
It remains to convert the absolute impedance levels 
to the correct order of magnitude and to scale the 
frequency of resonance of the loudspeaker to the 
proper value. This is, however, not to be discussed 
here. 


5. Results 


For a general presentation of the results it is con- 
venient to express relevant parameters of the system 
in terms of corresponding ones of the loudspeakers 
unit. This represents the usual problem in which the 
loudspeaker constants are known and the cabinet 
parameters are to be determined. However, the loud- 
speaker constants cannot be used as independent 
variables. Hence we assume the resonance frequency 
«,/2 of the loudspeaker unit to be known: 

o2= a 
MC 


Its damping, as expressed by the generalized 


Q-factor 
1 1/M 
ee alee 


will be determined by the bass-reflex design. 
Instead of A we propose to use as independent 
variable the near-field interaction parameter 


y=m/M. 


The interaction inductance depends on the distance 
between diaphragm and port but is independent of 
the port area [3]. Values of y up to unity have been 
used in the calculations. Since M contains diaphragm 
mass as well as radiation mass, such high values of 
y are physically impossible. A method to increase m 

A most important datum is the cut-off frequency 
@ 9/2 of the enclosure for a given resonance fre- 
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Fig. 5. Cut-off frequency of enclosure versus the inter- 
action parameter y. 


quency @,/2x of the loudspeaker. In Fig.5 the 
quotient « /@, is plotted against y . The continuous 
line relates to the Burrerwortu response, the dashed 
lines pertain to CHEBysHEV responses with various 
ripple factors as indicated. We conclude that allow- 
ing peaks in the response (CuEBysHEv) provide us 
with an advantage in the cut-off frequency as ex- 
pected. 

Fig. 6 shows the value of V/M as a function of y. 
Only one line is drawn which pertains to the Burrer- 


0.2 
001 O02 004006 04 02 
y¥=m/M —r 


0406 1.0 
Fig. 6. Inductance ratio versus y. 


worTH case. The values for the CuEBysHEV cases are 
approximately the same. This quantity indirectly 
determines the relative sizes of port and diaphragm 
areas. The relation is not unique, however. The 
arrow marks the region where both areas are appro- 
ximately equal. Larger values of N are obtained 
either by decreasing the port area or by providing a 
duct connecting to the port. From a theoretical point 
of view both methods are equivalent. 


04 06 10 
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Fig. 7. Capacitance ratio versus y. 
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The enclosure volume is shown as the ratio D/C 
in Fig. 7. In effect the volume is here expressed in 
terms of a unit volume determined by the loud- 
speaker stiffness. Such a unit volume, if employed 
as a closed box, would cause the loudspeaker’s res- 
onance frequency to be increased by a factor of 

V2. From Fig. 5 we infer that large values of y are 
not advantageous, whereas from Fig. 7 just the 
opposite conclusion can be drawn. To solve this 
ambiguity, the trading relationship between fre- 
quency range and size is shown in Fig. 8. 


TT 


Closed box 


2 
O04 0.2 O04 06 08 4 2 4 6 810 
0/C —~ 


Fig. 8. Frequency range (lower cut-off frequency (Wo) 
versus volume of enclosure (capacitance D). 


On the abscissa is plotted the enclosure volume, 
on the ordinate the cut-off frequency of the resulting 
bass-reflex structure. Both quantities are normalized 
as before. The line marked closed box serves as a 
landmark showing the relation that would be obtain- 
ed with a critically damped loudspeaker, in a com- 
pletely closed cabinet. The graph shows very clearly 
that a bass-reflex enclosure has advantages over the 
closed box. The improvement increases for the 
Chebyshev responses. The value of y determines the 
“working point” on the graph. Large values cor- 
respond to the smaller volumes, so that y should be 
made as large as possible. Such large values of the 


‘ Fig. 9. Suggested improvement to increase ¥ artificially. 


interaction inductance m can artificially be obtained 
by an internal duct (see Fig. 9). This expedient pro- 
vides an extra inductance in the enclosure branch, 
to be included in m. The parameters M-m and N-m 
now change their meaning somewhat since the 
terms-m here arise from true acoustic interaction 
only. 

Fig. 10 finally shows the required value of the 
generalized Q of the loudspeaker unit. It is seen to 


‘0.01 002° 004006 Oo” 0.2 O4 06 08 1.0 
y=m/M — 


Fig. 10. Quality factor of the loudspeaker’s resonant 
system required to provide the necessary damp- 
ing. 


be exceptionally low, which means damping should 
be very high. This is reasonable since in the case 
under study all damping must arise from the loud- 
speaker unit (Fig. 2a). One may regard this as the 
“cost” that has to be paid for incurring the advan- 
tage of the bass-reflex system. Such a high damp- 
ing can be obtained in several ways. It is paramount 
to use a loudspeaker with high electromagnetic 
damping (high efficiency), if possible provided with 
extra frictional damping. The remaining amount of 
damping can be obtained by applying motional feed- 
back (see Appendix III). In this case the sensitivity 
of the power amplifier to be used with the system is 
prescribed by the bass-reflex design as well. 


6. Conclusions 


With this method bass-reflex design is found to be 
remarkably simple when linear network synthesis 
is utilized. The values of the circuit elements are 
found without solving simultaneous nonlinear equa- 
tions. The resulting structure is geometrically realiz- 
able for not too large values of y. The procedure 
shows that true bass-reflex action is possible, i.e. 
a smooth response can be obtained via a highly res- 
onant enclosure. At the frequency of resonance of 
the enclosure the diaphragm radiates via the port 
under optimal load conditions, the impedance being 
infinite. As regards non-linear distortion this is the 
most advantageous situation. 

On each side of the cabinet resonance frequency 
there appear regions where the mechanical impe- 


4 
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dance is low. In the undamped case these cause the 
familiar maxima in the voice-coil electric im- 
pedance. They correspond in general to maxima in 
the sound pressure response. If the proper damping 
is applied, the response will be uniform without 
requiring abnormally large values of diaphragm 
velocity. The latter statement follows from a study 
of the input admittance of the network of Fig. 3. 

The large damping required can be obtained by, 
amongst other methods, motional feedback. In that 
case an important design datum of the amplifier, the 


sensitivity, is included in the acoustic design. Loud- 


speaker and amplifier are then to be regarded as 
one unit. 

The results from this analysis are purely theoreti- 
cal data. In this respect they are to be regarded as 
tentative. However, the procedure stands on a firm 
theoretical basis. The data found can be used as the 
starting point in the physical design of a bass-reflex 
enclosure. One interesting feature would be the 
effect of deviations from the correct values on the 
acoustic response. This would provide the designer 
with a measure of the required accuracy in the con- 
struction. The simplest method in this respect is the 
use of the electric analogue, realized as an adjust- 
able network, as an analogue computer. The re- 
sulting acoustic response can be obtained by a 
simple current measurement (or a voltage measure- 
ment in the dual circuit). It is then relatively simple 
to increase the complexity, e.g. to study the effects 
of radiation resistances, etc. Elsewhere it has been 
pointed out that in a more elaborate analysis, start- 
ing from a less simplified electric analogue, the dif- 
ficulties encountered will soon be prohibitive [3]. 
A step by step analysis will be called for [2] which 
tends to obscure the overall uniformity in the 
theory. 


Appendix I. The Daruineton problem 


Let the lossless coupling network of Fig. 2b be de- 
scribed by its admittance matrix (y). The elements y;;, 
Yio etc. are the coefficients in the following set of equa- 
tions 

ty =Y11 C1 + Ys Co, 

ig = Yo €1 + Yoo Co - (5) 
These equations describe the four-terminal network in 
the most general case. The network is loaded by a re- 
sistor which can be normalized to one ohm without loss 
of generality. Substitution of e,= —i, into eq. (5) and 
elimination of i, yields as the value of the transfer 
admittance 

eet Yor x 


aon 1+Ys0_ 


The zero’s of Y;5(s) are the zero’s of ys; (or the poles 
of y9). Hence the numerator of ys, is s®. Putting 


Yi =k 83/q(s) , 
Yoo = p(s)/q(s), 
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we get 

ae 
p(s) +4(s) 
Since p(s) and q(s) pertain to a lossless network, one 
(in this case gq) must be an even polynomial, the other 
an odd polynomial in s. Hence p(s) and q(s) are to be 
identified as the odd and even parts, respectively, of the 
denominator polynomial of Y;3(s). We then find 


Y13(s) = 


33 
143.414 52454’ 
2.613 s+2.613 s° 


Yee 73 Ald 24st | 


Ya =Yr=hk 


If we synthesize an LC-network with the correct driving- 
point admittance yy it is relatively easy to determine a 
second pair of terminals to which the correct shont-cir- 
cuit transfer admittance yj; pertains. This is especially 
true when Yp is realized as a ladder network. 

Any inductive shunt branch or capacitive series 
branch creates a transmission zero at s=0. The para- 
meter 1; is irrelevant in this discussion. However, it 
plays a role in the total input admittance of the loaded 
network. The synthesis part is carried out from the out- 
put side on, by realizing the output admittance yo. 
(with the input terminals of the coupling network short- 
circuited) . 

From Fig. 2b we see that the network has to start 
with a series element. Hence we form 1/yo9(s) and note 
that this for small frequencies behaves as 1/2.613s, 
and for large frequencies as s/2.613. Either of these can 
serve as the starting point of the synthesis. The correct 
choice is the first one. The remaining elements result 
from the continued fraction expansion: 


1 1 1 
ee PPAR ie 2 ae 
0:22 91 oem 
0.635 s 


The steps are illustrated in Fig. 11 a, b, and c. Fig. lid 
shows tthe result when the ladder is developed from the 
alternative starting point. It is easily seen that the net- 
work of Fig. llc provides for the correct number of 
transmission zero’s at s=0. 
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Fig. 11. Steps in the Daruieron synthesis. 


Appendix II. 
Internal impedance transformation 


Comparison of Figs. 2b and 3 shows that in the syn- 


thesis result one inductance is still lacking. We want to 
transform the inverted L-structure of two inductors into 


pee ea, SO ere "hola 
Repeat 
‘pan i 


ACUSTICA 


Frases E. DE BOER: BASS REFLEX LOUDSPEAKER ENCLOSURES 7 


a T-configuration. Let the inverted L-structure be de- 
scribed by its impedance matrix (z). As shown in 
Fig. 12 this four-terminal network is thought to be 
driven by an impedance z,4 and loaded by an impedance 


La 2p 


[2] [75 


Fig. 12. Representation of a circuit to perform an im- 
pedance transformation. 


zp. Then the transfer admittance Y,3(s) can be expres- 
‘sed in terms of z4, zg, and the elements of the matrix 
(z) as follows: 


— 712 ; (6) 


j= 
i (244241) (zB-++ 229) — 219” 


Suppose now that we want to change the impedance 
level of the first mesh with respect to the second mesh 
by a constant factor A?. 
The overall behaviour of the network should be in- 
variant, or at most be changed by a constant. Formula 
(6) shows that this can be obtained in the following 
way: 
i. 2,2 is multiplied by A, and 
ii. z4 and 2, are multiplied by A”, whereas 
iii. Zo9 and zg are not changed. 
By this procedure Y,3(s) is changed by a factor 1/4 
which is not important as far as the ultimate response 
is concerned. However, the transformation relieves the 
constraint Z12=Zo5 from the matrix (z). It then is no 
longer restrained to relate to a two-element network. 
The z-network can now be realized in a general T-con- 
figuration. The result of the transformation gives Fig. 4. 


Appendix III. Motional feedback 


In motional feedback the voltage induced in the moy- 
ing voice-coil winding is intended to be used as feedback 
voltage. It is then necessary to separate this induction 
voltage from the total voice-coil voltage. This is realized 
by a bridge circuit as shown in Fig. 13 a. Suppose the 
bridge is in equilibrium for all frequencies when the 
voice-coil is blocked, i.e. restrained from moving. If the 
voice-coil is allowed to vibrate, the voltage appearing 
across the horizontal diagonal of the bridge will be 
exactly proportional to the induction voltage. The latter 


19-9; 


is equal to the voice-coil current times the “motional 
impedance” zayn, the latter being the electrical image of 
the mechanical impedance of the diaphragm. Motional 
feedback implies the use of the output voltage of the 
bridge as feedback voltage. This tends to make the dia- 
phragm velocity constant, thus preventing large dia- 
phragm excursions such as arise at the resonance fre- 
quency. 

Since the dynamics is largely governed by zayn, one 
may regard this as focus. The whole system may be 
described as an “amplifier” that is loaded only with 
Zdyn- This hypothetical amplifier not only includes the 
original driving amplifier but also elements of the 
bridge and the voice-coil resistance. 

The resulting situation is depicted in Fig. 13 b. Mo- 
tional feedback is now equivalent to conventional volt- 
age feedback applied to this revised amplifying system. 
The voltage feedback will decrease the effective internal 
resistance of the amplifier, thus increasing the damping 
exerted on the loudspeaker. For the sake of clarity we 
may best depart from an amplifying system which a pri- 
ori does not damp the loudspeaker at all. The hypothe- 
tical amplifier acts then as a pure current source J, , 
driven by an input voltage e;. The sensitivity of the am- 
plifier is described by its mutual conductance g’ : 


/ /, 
Iv =e &. 


The effective g’ of this hypothetical amplifier may be 
lower than that of the actual driving amplifier due to 
bridge losses. If the amplifier cannot safely be regarded 
as a current source, the complete representation includ- 
ing an internal resistance R; can be used (shown dotted 
in Fig. 13:b). The addition does not materially affect 
the discussion. 

In the representation as a current-source amplifier, 
the effect of feedback is easily assessed. Application of 
voltage feedback cannot change the short-circuit current 
I). The only change in the equivalent circuit brought 
about by the feedback is the appearance of a shunt 
resistance. In the ideal case of full voltage feedback 
the newly introduced internal resistance is 1/g’ [7]. 
The resulting equivalent circuit is shown in Fig. 13 c. 
If the bridge causes a loss given by the factor k in 
Fig. 13a (k<.1) the resistance is 1/k g’. We see that 
the application of voltage feedback to a current-source 
amplifier causes a finite internal resistance to appear in 
the equivalent circuit. The damping effect of this 
resistance can be found when the electromagnetic con- 
version factor of the loudspeaker is known. For a re- 
quired amount of damping the value of g’ for the hypo- 


©) 


(fo) 1/q' 


Fe 


Zy = Ri 
; Fig. 13. Schematic representation of the 
rath it motional feed-back problem. 
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thetical amplifier is determined. In its turn this corre- 
sponds to a minimum value of the sensitivity of the 


actual driving amplifier. 
(Received 30th May, 1960.) 
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WELLENAUSBREITUNG IN UNBEGRENZTEN SCHEIBEN 
UND IN SCHEIBENSTREIFEN 


von K. Tamm * und O. Wets* 
IIL. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Es wird gezeigt, daB die Wellenausbreitung in unendlich ausgedehnten Scheiben und in 
Scheibenstreifen (ebener Spannungszustand) analog ist zur Wellenausbreitung im unend- 
lichen Raum bzw. in Platten (ebener Verzerrungszustand). An die Stelle der Dilatations- 
welle tritt bei Scheiben die Plattendehnwelle; die Schubwelle wird nicht verandert. 

Die Anwendung dieser Analogie erméglicht zum Beispiel die Beschreibung der Schrig- 


reflexion von Plattendehnwellen und die Berechnung von ,,Rayiricu-Wellen* 


auf freien 


Scheibenrandern aus der Kenntnis der analogen Vorgange an einer Halbraumoberfliche. 
Kine eingehende Untersuchung befat sich mit den Wellenformen in einem Scheibenstreifen 
aus gummielastischem Material beziehungsweise — in Analogie — mit denen in Platten 
der Querkontraktionszahl 1/3. Die Beriicksichtigung von Materialverlusten ergibt neben 
Aussagen tiber den Dampfungsverlauf einzelner Wellenformen auch eine Erklarung, wes- 
halb bei hohen Frequenzen yon der Vielfalt der berechenbaren Wellen neben den Rayiercn- 
Wellen nur die Schubwellen und die Dilatations- beziehungsweise Plattendehnwellen ge- 


messen werden. 


Summary 


It is shown that there is an analogy between wave propagation in infinitely extended disks 
or disk strips (plane state of stress) and wave propagation in infinite space or in plates 
respectively (plane state of distortion). In the case of disks the dilatational wave is replaced 
by an extensional plate wave; the shear wave type remains unchanged. 

The application of this analogy renders possible the description of the oblique reflection 
of extensional plate waves and the calculation of “RayteicH waves” at free boundaries of 
disks from the knowledge of the analogous processes at a semispace surface. A detailed 
investigation deals with the wave types in disk strips consisting of rubber elastic material 
or —by analogy —with those in plates with a Poisson’s constant 1/3. From consideration of 
losses follow the attenuation of various wave types and an explanation of the experimental 
finding that at high frequencies — besides Rayteren waves — only shear, dilatational or exten- 
sional plate waves respectively are found from the multitude of wave types which are theo- 


retically possible. 


Sommaire 


On montre que la propagation d’ondes dans des disques infinis ou dans des bandes (effort 
plan) est analogue 4 la propagation dans l’espace infini ou dans les plaques (allongement 
plan). Au lieu de l’onde de dilatation dans les disques on a l’onde d’allongement des pla- 


ques. L’onde de cisaillement n’est pas changée. 


L’application de cette analogie permet, entre autres exemples, d’expliquer la réflexion 
oblique des ondes Pallonpement des plaques et de calculer les ondes de Raytetex sur les 
bords libres des disques lorsqu’on connait les phénoménes analogues 4 la surface de l’espace 
semi-infini. Une intéressante recherche est celle de la forme des ondes dans une bande de 
disques dont la matiére a l’élasticité du caoutchouc comparée 4 celle des plaques dont le 
coefficient de concentration transversale est de 1/3. La considération des pertes dans les 
substances renseigne sur les pertes par amortissement des ondes particuliéres et explique 
pourquoi en haute fréquence parmi toute la variété des ondes possibles, on ne peut mesurer 
que les ondes de Rayzeten, les ondes de cisaillement et les ondes de dilatation ou d’allonge- 


ment des plaques. 


* Zur Zeit Institut fiir Angewandte Physik der Universitat Heidelberg. 


2h 


ACUSTICA 
Vol. 11 (1961) 


1. Einleitung 


In der Statik der Festkorper gibt es zwei tiber- 
sichtliche und relativ einfach zu berechnende Feld- 
verteilungen: ebene Verzerrungsfelder und ebene 
Spannungsfelder. 

Beim ebenen Verzerrungsfeld treten nur Verzer- 
rungen in einer ,,Ebene“ auf; das heiSt in der Be- 
zeichnungsweise des Bildes 1, da} nur die Verzer- 
rungen €;,, €,, und €,, vorhanden sind. é;, , €yz 


und ¢,, sind Null. Dies hat zur Folge, da auch die 


'Schubspannungen t,, und t,, verschwinden und da 


weiterhin keine Teilchenverschiebung in z-Richtung 
sz auftritt. Eine Normalspannung 6;, ist zugelassen. 


7 oe 


Eyy Ouy af 
Exy Try Try 
Exx Oxy 
= |b, On, cl. 
Ex, 
Eyy Oyy 


Bild 1. Gegeniiberstellung von ebenem teat aes 
(a) und ebenem Verzerrungszustand (b). 


. Wie beim ebenen Verzerrungsfeld sind auch beim 
ebenen Spannungsfeld alle FeldgroBen von der 
z-Koordinate unabhangig. Ebene Spannungsfelder 
besitzen nur Spannungen in einer ,,Ebene“; ist dies 
die xy-Ebene, so ist also 6,,, t;; und t,, durchweg 
Null. Eine Verschiebung s, ist moglich, Ein solcher 
ebener Spannungszustand tritt in den sogenannten 
»ocheiben“ auf; das sind diinne Platten, die in der 
Plattenebene mit tiber den Querschnitt konstanten 
Kraften belastet werden. Der Unterschied zwischen 
beiden Feldern liegt also praktisch darin, daB beim 
ebenen Spannungszustand das beanspruchte Mate- 
rial in Querrichtung (z-Richtung) ungehindert aus- 
weichen kann !, wahrend beim ebenen Verzerrungs- 
zustand diese Querkontraktion verhindert wird. 

Ein solches Bestreben, seitlich nachzugeben, tritt 
in unseren beiden Fallen nur dann auf, wenn die 
angreifenden Krafte eine Kompression und damit 
eine Volumenanderung des Materials bewirken — 
es sei denn, da das Material von sich aus keine 
Querkontraktionszahl besitzt, das heift die Querkon- 
traktion o Null ist. Haben wir fiir jedes Volumen- 
element nur eine Formanderung ohne gleichzeitige 


1 Es sind nur zur Mittelebene symmetrische Verschie- 
bungen zugelassen, da eine Biegung eine antisymme- 
trische und damit nicht konstante Spannungsverteilung 
iiber den Querschnitt bedingen wiirde. 


oa 
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Volumenanderung, das heift haben wir eine reine 
Schubbeanspruchung, so geht die Korpererstreckung 
in z-Richtung gar nicht ein; o;;, 
s, treten nicht auf. Bei reiner Schubbeanspruchung be- 
steht kein Unterschied zwischen ebenen Spannungs- 
und ebenen Verzerrungsfeldern. Das hat zur Folge, 
daf} die Ausbreitung von Schubwellen (mit Verschie- 
bungsvektoren in der xy-Ebene) im unendlich aus- 
gedehnten Medium und in Scheiben gleich ist. 
Neben den Schubwellen (divergenzfrei) existie- 
ren im unbegrenzten Medium noch die Dilatations- 
wellen (rotationsfrei). Jedes dynamische ebene Ver- 
zerrungsfeld entsteht durch Uberlagerung dieser bei- 
den primaren Wellen, deren Verschiebungsvektoren 
in der xy-Ebene liegen mtissen. Treten spannungs- 
freie Oberflachen oder andere Storungen senkrecht 
zur xy-Ebene auf, so wandeln sich Schub- und Dila- 


€z, beziehungsweise 


tationswelle an diesen Stellen zum Teil ineinander 


um: Es tritt eine Kopplung auf, die zu neuen sekun- 
daren Wellen fuhrt. Solche sekundaren Wellen mit 
ebenem Verzerrungsfeld sind z.B. die Raytetcu- 
Wellen an der Oberflache eines Halbraumes und die 
Vielzahl der Wellenformen in einer unendlich aus- 
gedehnten Platte [1]. 

Eine allgemeine Behandlung der Schallausbrei- 
tung in Scheiben (bei ebenem Spannungszustand) 
scheint bislang noch nicht vorzuliegen, obwohl ex- 
perimentelle Untersuchungen an diesen oft viel ein- 
facher auszufuhren sind als bei entsprechenden Pro- 
blemen mit ebenem Verzerrungszustand. Wir wer- 
den im folgenden sehen, dafs bei den Scheiben die 
gleichen Wellenformen auftreten wie bei den oben 
angeftihrten sekundaren Wellen mit ebenem Verzer- 
rungsfeld. Das liegt daran, daf} bei Scheiben an die 
Stelle der in xyz rotationsfreien Dilatationswelle die 
allerdings nur in xy rotationsfreie Plattendehnwelle 
tritt, die formal in Scheiben dieselben Eigenschaf- 
ten besitzt wie die Dilatationswelle im unendlichen 
Medium. Die Plattendehnwelle selbst ist natiirlich 
als sekundare Welle eine Kombination einer Schub- 
und einer Dilatationswelle. 


2. Wellenausbreitung in einer unendlich 
ausgedehnten Scheibe 


Wir gehen von den allgemeinen Grundgleichun- 
gen eines homogenen elastischen Mediums mit klei- 
nen Teilchenverschiebungen aus und spezialisieren 
diese auf Scheiben, indem wir beriicksichtigen, daB 
in Scheiben alle Spannungen tiber die Scheibendicke 
konstant sind, das heift, da bei ihnen keine Ab- 
hangigkeit von der z-Koordinate auftritt. Da die 
Spannungen 6,,, ty, und T;, an der spannungsfreien 
Scheibenoberflache verschwinden miissen, sind sie 
durchweg Null. 


Die allgemeinen Bewegungsgleichungen (ohne 


ne A 
ce 


mat 


=F Oe anon 
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Beriicksichtigung der Volumenkrafte) 
Qozz Oey A OT ze 0 O57 a. 

Or Oy ' oO Oo? : 

OFpy OG yy Ol ye SHS 

an 8 oy a 


Oz . Oz? 
Otn Oye 6x2 o ois 
Ox Oy Oz 02? 
reduzieren sich also bei Scheiben auf 
Se 
Or | Oy Ohi. 

Oty OGy) 
poy ser 


Die Materialgleichungen 


On (la) 


a(n 


+ 4 + 


(1b) 


pose = E &y7 = Oz — 6 (Oyy + Gzz), 


Ox 


Qs, 
ES 


=E Eyy = Syy — 6 (O72 + Oza) 


G (See + Se) rs ae (2a) 


=2 GC Eg2=T zz > 


vereinfachen sich ebenfalls: 


Bo tess Os, 
t= o2 (= ee =") - 


Org 


Dabei bedeuten die 0; , 7, Spannungen, die &jj, éix 
Verzerrungen in der Bezeichnungsweise von Bild 1. 
s; sind die Verschiebungskomponenten, @ ist die 
Dichte, o die Querkontraktionszahl, E der Young- 
sche Modul, G=E/2(1+.0) der Schubmodul und 
P=E/(1—o?) der Plattendehnmodul. 
Die noch nicht beriicksichtigte Beziehung 

ie = Ozz — O(Ogr + Syy) 
von Gl. (2 a) ermoglicht uns noch, die Verschiebung 
s, zu berechnen. Da diese Verschiebung antisym- 
metrisch zur Plattenmitte sein mu, kann man 
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S, = const*z ansetzen und erhalt damit 
Tata 
Sai E (Gzz + Syy) Zz. (3) 


Setzt man Gl. (2b) in Gl. (1b) ein, so ergibt ) 
sich die Wellengleichung 


sed g 
+65, (ae = oe oe 
cee 2 


Diese Wellengleichung kann man einfacher in Vek- 
torform schreiben: 


(4b) 


P grad, div, 8 —G rot, rot, §+0 w? 8=0, 


wobei der Index 2 anzeigen soll, daB die Vektor- 
operationen nur im Zweidimensionalen, das heift nur 
beziiglich x und y vorzunehmen sind. Eine Aufspal- 
tung des Verschiebungsvektors 8 = (s,,s,) in einen 
in xy rotationsfreien Anteil ss und einen in zy 
divergenzfreien Anteil 8p fiihrt auf zwei voneinan- 
der unabhangige Wellengleichungen. Es ergibt sich 
fiir div. $= 0 die Wellengleichung 


ea 
Aga, BOS PG (5a) 


As Tee eae (5b) 
Losungen der Wellengleichung (5a) sind Schub- 
wellen (divergenzfrei) mit einer Ausbreitungs- 
geschwindigkeit cs = VG/o im verlustlosen Medium. 
Lésungen der Wellengleichung (5b) sind Platten- 
dehnwellen (lediglich in xy rotationsfrei) mit der 
Ausbreitungsgeschwindigkeit cp = VP/e. 

Formal unterscheiden sich die Wellengleichungen 
(5) und die Materialgleichungen (2b) der Scheiben 
von den gewohnten des unendlich ausgedehnten 
Mediums in zwei Dingen. Erstens sind Gl. (5) und 
(2b) auf zwei Dimensionen (2, y) beschrinkt, was 
jedoch bei den ebenen Verzerrungsfeldern durch 
die spezielle Art des Wellenfeldes auch zutrifft und 
zweitens — und das ist der entscheidende Unter- 
schied — tritt an die Stelle des Dilatationsmoduls 
D=E(1—o)/[(1+0)(1—20)] der Plattendehn- 
modul P=E/(1—o?). Die Dilatationswelle mit 
der Ausbreitungsgeschwindigkeit cp—VD/o (ver- 
lustfrei) wird also durch die Plattendehnwelle mit 
cp=VP/o (verlustfrei) ersetzt. Diese Ersetzung 
wird beinahe selbstverstandlich, wenn man sich die 
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Verformungen und Beanspruchungen vergegenwar- 
tigt, welche die Moduln D, P und E£ beinhalten und 
mit denjenigen beim ebenen Verzerrungs- und beim 
ebenen Spannungszustand vergleicht: Haben wir an 
einem Volumenelement eine einzige Normalspan- 
nung, zum Beispiel o,, und eine zugehérige Deh- 

“nung €,,, so wird der Zusammenhang zwischen 
Spannung und Dehnung — wie man leicht nachrech- 
nen kann — bei allseitig verhinderter Querkontrak- 
tion 


(unendliches Medium) 


bei einseitig verhinderter Querkontraktion 


Oz2 = D Ezz » 


(unendliche Scheibe) OP Peres 
bei unbehinderter Querkontraktion 
(dunner Stab) One =F Ebner 


Der Ubergang vom ebenen Verzerrungsfeld zum 
ebenen Spannungsfeld mit gleicher Berandung und 
gleichen auferen angreifenden Kraften erfolgt also 
einfach durch eine Modultransformation, indem 
G/D ersetzt wird durch G/P . Damit ist der Verlauf 
der Feldgréfen in Scheiben mit der Querkontrak- 

_tionszahl o analog dem Verlauf im ebenen Verzer- 
rungsfeld mit der Querkontraktionszahl G6, wenn 
das Verhaltnis G/P = (1 — 0) /2 mit dem Verhiltnis 
G/D = (1—26)/[2(1—G)] ibereinstimmt?. Das 
ist der Fall, wenn die Beziehung 


oO 


ee (6) 


C= beziehungsweise G= 


1—G 


besteht. Besitzt also z. B. das Scheibenmedium die 
Querkontraktionszahl o=0,5 (0,25), so liegt die 
gleiche Feldverteilung vor wie beim ebenen Verzer- 
rungszustand in einem Material der Querkontrak- 
tionszahl =1/3 (0,2). Bei verschwindender Quer- 
kontraktionszahl o=G=0 stimmen die Felder, wie 
schon friher erwartet, tiberein. Die Feldverteilun- 
gen der Scheiben haben also ein Analogon bei den- 
jenigen der ebenen Verzerrungsfelder. Umgekehrt 
gilt dies jedoch nicht immer, da der Bereich 
o6=0,5...0 der Scheibe auf den kleineren Bereich 
6=1/3...0 des ebenen Verzerrungsfeldes ab- 
gebildet wird. 


3. Reflexion von Platten-, Dehn- und 
Schubwellen an einem freien Scheibenrand 


AurenBerc [2] hat eingehend die Reflexion von 
Schub- und Dilatationswellen an einer spannungs- 
freien Halbraumoberflache untersucht und Dia- 
gramme in der Art des Bildes 2 und 4a angegeben. 


? Bei Scheiben tritt noch die zusitzliche Verschiebung 
sz (nach Gl. (3)) auf. Sie wird nur durch die Platten- 
dehnwelle hervorgerufen, kann also zur Ermittlung die- 
ses Anteils benutzt werden. 
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Mit Hilfe unserer Analogie nach Gl. (6) konnen 
wir nun leicht auch die entsprechenden Diagramme 
fiir die Reflexion von Schub- und Plattendehnwel- 
len an einem freien Scheibenrand finden (Bild 3 
und 4b). Wir sehen, dafS§ bei der Schragreflexion 
einer Plattendehnwelle im allgemeinen gleichzeitig 
eine Schubwelle mit entsteht und umgekehrt — 


Os 


03 


0,5 
0° (Ope 20 age sO ee On 25am nOOL Ope seO a= 90: 
Sy ee od 


Bild 2. Reflexionsfaktor R (bezogen auf die Verschie- 
bungsamplitude beziehungsweise auf die 
Schnelle) in Abhangigkeit vom Einfallswinkel # 
einer Schubwelle (a) und einer Dilatationswelle 
(b) bei der Reflexion an einer spannungsfreien 
Halbraumoberflache. Parameter ist die Quer- 
kontraktionszahl. 

Ist |R|<1, so geht ein Teil der einfallenden 
Energie in die andere Wellenart iiber. Die 
Schubwelle besitzt einen Grenzwinkel der Total- 
reflexion; bei gréferem Einfallswinkel ist die 
mit angeregte Dilatationswelle quergedémpft 


und wird nicht abgestrahlt (Nahfeld). 
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POOP eiZO- 30% Oa Sn SCOw th 10 mee COs eR Os 


>— 
Bild 3. Reflexionsfaktor R (bezogen auf die Verschie- 
bungsamplitude beziehungsweise auf die 


Schnelle) in Abhangigkeit vom Einfallswinkel 0 
einer Schubwelle (a) und einer Plattendehn- 
welle (b) bei der Reflexion an einem freien 
Scheibenrand. Parameter ist die Querkontrak- 
tionszahl. 

Ist |R|<1, so geht ein Teil der einfallenden 
Energie in die andere Wellenart iiber. Die 
Schubwelle besitzt einen Grenzwinkel der Total- 
reflexion; bei gré8erem Einfallswinkel ist die 
mit angeregte Plattendehnwelle quergedimpft 
und wird nicht abgestrahlt (Nahfeld) . 


ganz entsprechend, wie es von Dilatations- und 
Schubwellen her bekannt ist. 


So, wie es an einer Halbraumoberflache Rayiercu- 
Wellen gibt, die sich aus einer quergedaémpften 
Dilatationswelle und einer quergedémpften Schub- 
welle zusammensetzen, so gibt es auch RayLEIGu- 
Wellen an einem Scheibenrand. Sie bestehen hier 
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Bild 4. Wechselwellenwinkel (= Einfallswinkel, bei 
denen eine vollstandige Umwandlung der einen 
Wellenart in die andere auftritt) bei der Refle- 
xion an einer spannungsfreien Halbraumober- 
fliche (a) und einem spannungsfreien Scheiben- 
rand (b) in Abhadngigkeit von der Querkontrak- 


tionszahl. 
Wechselwellenwinkel fiir einfallende 
Schubwelle, 

— — — Wechselwellenwinkel fiir einfallende 
Dilatations- beziehungsweise Platten- 
dehnwelle. 


AuBerdem sind noch die Grenzwinkel der Total: 
reflexion einer einfallenden Schubwelle 
eingetragen. 


aus einer quergedampften Plattendehnwelle und 
einer quergedampften Schubwelle, die wiederum 
beide zu einer Kinheit verkoppelt sind. Bild 5 zeigt, 
wie gummielastisches Material bei diesen Rayercn- 
Wellen verformt wird. Bild 6 gibt die Abhangigkeit 


et at 'y 
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Bild 5. 


Ubertrieben gezeichnete De- 
formation yon gummielasti- 
schem Material (o = 0,5) beim 


Auftreten yon Ray.eicu-Wel- 


len an einem Scheibenrand. 
Die Verschiebungsamplituden 


senkrecht zur Scheibenebene 
in den Schnitten AA, BB, CC 


Teese Ee 
ogee ee a 
Dassen seeesee=s i 
Bee ee a 
eee a 
PEE EERE EEE ECE = 
PERE rr eer 2 
Hoes conan a 
PP im 
Eee Es 
Pays cay eegueagneiee 
Bi Repel ak ‘a ee 
ches ecanoeee a 


sind gegentiber den anderen 


0,95 
f Halbraumoberflache 
o Scheibenrand 
090 
=) 
0.85 | 


fae) 01 02 03 04 05 
Querkontraktionszahl —> 


Bild 6. Abhingigkeit des Verhaltnisses Rayercu-Wel- 
¢ lengeschwindigkeit zu Schubwellengeschwindig- 
keit yon der Querkontraktionszahl. 


der Phasengeschwindigkeit der Rayiercu-Wellen 
von der Querkontraktionszahl wieder. 


4. Wellenausbreitung in einem Scheibenstreifen 
(Scheibe mit zwei parallelen freien Randern) * 


Das Problem der Wellenausbreitung in einer 
Scheibe, die durch zwei freie, zueinander parallele 
Rander begrenzt ist, (kurz ,,Scheibenstreifen“ ge- 
nannt), ist analog der Wellenausbreitung in einer 
unendlich ausgedehnten Platte. Dabei entspricht die 
Scheifenstreifenbreite der Plattendicke (vgl. Bild 7). 


Bild 7. Zur Bezeichnungsweise bei Scheibenstreifen und 
bei Platten. 


3 Experimente zur Wellenausbreitung in Scheiben- 
streifen wurden bereits in den Jahren 1946/47 von 
W. Kunt und E. Meyer [4] ausgefiihrt. Sie gaben den 
Anla8 zu dieser theoretischen Untersuchung. 


ee 


Verschiebungen 10-mal gro- 
er gezeichnet, um sie tiber- 
haupt sichtbar zu machen. 


Wenn wir die Wellenausbreitung in einem Schei- 
benstreifen untersuchen, so hei®t das mit anderen 
Worten, wir untersuchen die Ausbreitung in einem 
mehr oder weniger breiten (Erstreckung in y-Rich- 
tung in Bild 7), unendlich langen Stab mit recht- 
eckigem Querschnitt in Abhangigkeit von l/A;, das 
hei®Bt vom Verhaltnis Stabbreite zu Schubwellen- 
lange; jedoch soll die Stabdicke (Erstreckung in 
z-Richtung in Bild 7) immer klein gegen die Schub- 
wellenlange bleiben. 

In Analogie zur Platte sind in y symmetrische 
und antisymmetrische Wellenformen zu erwarten, 
die fiir sehr groBe 1/As-Werte, das heiBt sehr breite, 
flache ,,Stabe“, in. die beiden Wellenarten der un- 
endlich ausgedehnten Scheibe tibergehen miissen, 
namlich in die Schub- und Plattendehnwelle. Wir 
werden sehen, da diese Ubergange jedoch erst bei 
Beriicksichtigung der Materialverluste voll sichtbar 
werden. 

In einer fritheren Arbeit [1] wurden alle mog- 
lichen, in einer unendlich ausgedehnten, verlustlosen 
Platte auftretenden Wellenformen ausfihrlich dis- 
kutiert. Spater [3] wurden Materialverluste beriick- 
sichtigt und ihr Einflu8 auf die Plattenmoden im 
Falle eines gummielastischen Materials (0 =0,5) 
untersucht. 

Zur Bestimmung der Wellenausbreitung in einer 
Platte, die in 2-Richtung proportional e~% e'” er- 
folgen soll, liefern die Bedingungen der Spannungs- 
freiheit auf der Plattenoberflache eine Beziehung 
zwischen der Fortpflanzungskonstante ok le oR 
und dem Verhiltnis //As. Dabei spielen die Mate- 
rialkonstanten G/D, der Verlustfaktor fiir Schub 
tan Os und fiir Kompression tand0x die Rolle von 
Parametern. Mit der zweckmaBigen Abkiirzung 


2a —— = (7) 
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lautet die Bestimmungsgleichung * 
tanh V/z?—G/D(1 —itan4/2-ds) x ig) - 
te V2—1 (1—itan4/2-ds) n1/ds| 
(2? — 0,5)? x 
ze Vel V2z2—G/D 
ES | tanh V22—1 (1 —itan4/2-ds) x1/As| 
| tanh Vz? — G/D(1 —itan !/2-ds) x i/asf 


(8) 


Um den Real- und Imaginarteil von a beziehungs- 
weise von z und damit die Dampfung und Phasen- 
geschwindigkeit der einzelnen Wellen zu erhalten, 
mufX{ man Gl. (8) nach z auflosen, das heift in die 
Form 


z=2(I/ds , G/D, tan dg, tan dx) 


bringen. Die Auflosung ist nicht explizit moglich; 
sie kann mit einer elektronischen Rechenmaschine 
 geschehen (vergleiche [3]). 

Fir die Berechnung der Wellenausbreitung in 
einem Scheibenstreifen konnen wir wieder von der 
Beziehung (8) ausgehen; nur ist jetzt der Dilata- 
tionsmodul D=K+4/3-G (K = Kompressions- 
modul) durch den Plattendehnmodul 


P=4G6(3K+G)/(3K+4G) 


zu ersetzen. 

Im Bild 8 ist die Dampfung und Phasengeschwin- 
digkeit aller Wellenformen in Scheibensireifen aus 
verlustlosem, gummielastischem Material (o=0,5) 
in Abhangigkeit von J/As;— das hei8t praktisch von 
der Frequenz — dargestellt. Zugleich ist dieses Bild 
auch giiltig fiir die Wellen in einer verlustlosen 
Platte der Querkontraktionszahl o = 1/3. Dieses auf 
den ersten Blick recht verwirrende Bild der Lésun- 
gen wird etwas tbersichtlicher, wenn man gering- 
fiigige Verluste mit beriicksichtigt. Dann ergibt sich 
eine zwangslaufige Zuordnung der einzelnen Lo- 
sungsaste, wie sie in den Bildern 9 und 10 dar- 
gestellt ist. 

Eine Diskussion der einzelnen Losungen, das 
hei®t von komplexen Nahfeldern, von rein exponen- 
tiellen Nahfeldern und periodischen Wellen, wurde 
bereits frither [1], [3] durchgefiihrt..Die gewahlte 
Bezeichnungsweise kniipft an die Eigenschaften bei 
tiefen Frequenzen (kleinen l/As-Werten) an. Biege- 
welle 4), Dehnwelle Sy und 1. hohere antisymme- 
trische Mode A, stellen tiber den gesamten l/As- 
Bereich jede fir sich eine physikalisch vollstandige 
Loésung dar, die wir mit Mode bezeichnen wollen. 
Demgegeniiber beginnen alle anderen Lésungen bei 
kleinen J/As-Werten als komplexe Nahfelder; setzen 
sich also aus je einer hin- und riicklaufenden Welle 
zusammen, die in ihrer Uberlagerung erst physika- 


4 Die Bestimmungsgleichung (8) wurde in [3] ver- 
sehentlich nicht richtig wiedergegeben. 
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Bild 8. Dampfung (a) und Phasengeschwindigkeit (b) 
der einzelnen Wellenformen in einem Scheiben- 
streifen mit der Querkontraktionszahl o=0,5 
beziehungsweise in einer Platte mit o=1/3 in 
Abhangigkeit von l/As, dem Verhiltnis Strei- 
fenbreite zu Schubwellenlange beziehungsweise 
Plattendicke zu Schubwellenlinge. Keine Mate- 
rialverluste.  ( symmetrische, 
antisymmetrische Lésungen fiir die rein periodi- 
schen Wellen und rein exponentiellen Nahfel- 
der; symmetrische, — — — antisymme- 
trische Lésungen der komplexen Nahfelder.) 


lisch sinnvoll werden [3]. In jenen J/As-Bereichen, 
wo ein komplexes Nahfeld die Losung bildet, muh 
man unbedingt den einen Modenanteil (Bild9) — 
der eine Phasenausbreitung gleichsinnig dem Am- 
plitudenabfall besitzt und mit gleichgerichteten Pfei- 
len bezeichnet wird — mit dem anderen Moden- 
anteil (Bild 10) — der eine Phasenausbreitung ent- 
gegengesetzt zum Amplitudenabfall besitzt und mit 
gegensinnigen Pfeilen bezeichnet wird — mit glei- 
cher Amplitude im verlustfreien Fall superponieren, 
um-so zu einer physikalisch sinnvollen Mode zu 
kommen. Bei Beriicksichtigung groBerer Verluste 
(Bild 11) wird sich das Verhaltnis der Amplituden 
bei der Uberlagerung der Modenanteile in Abhan- 
gigkeit von I/As andern; jedoch wird nach wie vor 
eine Superposition notig sein. 

Jede Mode und jeder Modenanteil besitzt fiir 
jeden //As-Wert eindeutig einen Dampfungswert und 
eine bestimmte Phasengeschwindigkeit. Der Verlauf 
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Bild 9. Dampfung (a) und Phasengeschwindigkeit (b) 
jener Moden und Modenanteile des Bildes 8, die 
sich bei beliebig kleinen — aber endlichen Ver- 
lusten — ergeben und bei kleinen //As-Werten 
eine Phasenausbreitung in Richtung des Ampli- 
tudenabfalles besitzen. 


von Phasengeschwindigkeit und Dampfung in Ab- 
hangigkeit von J/As wird durch die Verluste zum 
Teil stark verandert. Die Diagramme des Bildes 11 
zeigen, da auch bei verlustbehafteten Platten mitt- 
lerer Querkontraktionszahl und bei Scheiben die 
gleichen Erscheinungen auftreten, wie sie friiher [3] 
fur verlustbehaftete Platten mit o=0,5 diskutiert 
wurden. Dazu gehort der gleiche charakteristische 
Dampfungsverlauf einzelner Moden und Moden- 
anteile, der Austausch von Teilen der Loésungs- 
kurven zwischen den einzelnen Moden und Moden- 
anteilen und die Unendlichkeitsstelle der Phasen- 
geschwindigkeit von A,. Daneben behalten nun 
auch noch alle jene Modenanteile ihre Unendlich- 
keitsstelle der Phasengeschwindigkeit, die bei klei- 
nen Frequenzen eine Phasenausbreitung entgegen- 
gesetzt dem Amplitudenabfall zeigen. Diese Un- 
endlichkeitsstellen miissen erhalten bleiben, weil die 
entsprechenden Lésungen bei hohen 1/As-Werten 
eine mit dem Amplitudenabfall gleichsinnige Pha- 
senausbreitung besitzen. Dies sind in Platten gerade 
jene Moden, die, wenn man die Querkontraktions- 
zahl gegen 0,5 (As/Ap gegen Null) laufen 1a8t und 
Ay endlich halt, in die Kompressionswellenmoden 
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l/s —> 
Bild 10. Dampfung (a) und Phasengeschwindigkeit (b) 
und (c) jener Moden und Modenanteile des 
Bildes 8, die sich bei beliebig kleinen — aber 
endlichen Verlusten — ergeben und bei kleinen 
1/As-Werten eine Phasenausbreitung entgegen 
der Richtung des Amplitudenabfalles besitzen. 
In (c) ist Amplitudenabfall und Phasenaus- 
breitung gleichsinnig, in (b) gegensinnig; es 
tritt also ein Wechsel der Phasenausbreitungs- 
richtung — vyerglichen mit der Richtung des 
Amplitudenabfalles — bei groBer werdendem 
l/As-Verhaltnis auf. 
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einer _ schubspannungsfreien, 
schallweich begrenzten Schicht 
iibergehen. 

Am Bild 11 kann man noch 
gut studieren, wie es kommt, 
daB bei hohen Frequenzen 
(groBen J/As-Werten) in der 
Praxis neben der RAyLEIGH- ~ 


Welle nur noch die Schubwel- 


2 len und die Dilatations- bezie- 


hungsweise _Plattendehnwelle 
ubrigbleiben. Das Auftreten 
der RayuetcH-Welle bedarf 
kaum  weiterer Erklirung, 
nachdem die Biegewelle Ay 
und bei nicht allzu hohen Ver- 
lusten die Dehnwelle in sie ein- 
miunden. Eine Erklarung er- 
iibrigt sich auch fiir die Schub- 
welle, in die normalerweise 
alle anderen Moden einmiin- 
den. Anders liegt es bei der 
Dilatationswelle und entspre- 
chend fiir Scheiben bei der 
Plattendehnwelle. Zwar verwei- 
len bei verlustfreier Rechnung 
mit steigendem //As die einzel- 
nen Moden immer langer bei 
der entsprechenden Phasen- 


Bild 11. 

Diagramme zur Wellenausbrei- 
tung unter Beriicksichtigung von 
Materialverlusten in einem Schei- 
benstreifen (o = 0,5; tan dg = 0,2; 
tan Ox beliebig) beziehungsweise 
in einer Platte (o=1/3; tan ds 
=tan 0x =tan dp=0,2). Dia- 
gramm (a) gibt die Dampfung 
aller Wellenformen an, Dia- 
gramm (b) die Phasengeschwin- 
digkeit derjenigen mit Phasen- 
ausbreitung in Richtung des Am- 
plitudenabfalles, Diagramm (c) 
die Phasengeschwindigkeit der- 
jenigen mit Phasenausbreitung 
entgegengesetzt zur Richtung des 
Amplitudenabfalles. ( 
symmetrische, ——— antisym- 
metrische Moden und Moden- 
anteile ohne Verluste; 
symmetrische, — —— antisym- 
metrische Moden und Moden- 
anteile bei Materialverlusten; 
wae a Dampfung einer ebe- 
nen Dilatationswelle beziehungs- 
weise einer ebenen Plattendehn- 
welle; —.-—..— Dampfung einer 
ebenen Schubwelle.) 
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Bild 12. 
Moden und Modenanteile ge- 
ringster Dampfung nach 


3 
Vi 


geschwindigkeit, doch wenden sie sich spater immer 
der Schubwelle zu. Durch Beriicksichtigung der 
Verluste wird das anders. Im Bild 12 sind aus dem 
Bild 11 jene Moden und Modenanteile heraus- 
gezeichnet worden, welche die geringste Dampfung 
aufweisen. Man sieht, da von einem bestimmten 
1/As-Wert ab fiir eine bestimmte Mode kein Abbie- 
gen der Schubwelle mehr auftritt; sie schmiegt sich 
vielmehr asymptotisch von nun an in ihrem Damp- 
fungs- und Phasengeschwindigkeitsverlauf demjeni- 
gen einer Dilatations- bzw. Plattendehnwelle an. 
Zuerst tritt dies fiir die symmetrischen und dann fiir 
die antisymmetrischen Moden auf. Dabei tritt auch 
der bereits erwahnte Austausch von Teilen von 
Losungskurven in Erscheinung. Bei den _ relativ 
hohen Verlustfaktoren, die den Bildern 11 und 12 
zugrunde liegen, findet der Ubergang schon bei 
S, und bei A,= statt. Alle anderen Moden zeigen 
eine groBere Dampfung; sie verschwinden rascher 
als jene asymptotisch in die Dilatations- beziehungs- 
weise Plattendehnwelle tibergehenden Moden. Der 
Dampfungsabstand wird um so grofer, je groBer 


Bild 11. 


I/As wird. Sind die Verlustfaktoren kleiner, so tritt 
der asymptotische Ubergang bei entsprechend hohe- 
ren J/As-Werten auf. 

Herrn Prof. Dr. Dr.-Ing. E. h. E. Meyer danken 
wir ftir die stete Forderung dieser Arbeit. Besonde- 
rer Dank gebiihrt auch der Deutschen Forschungs- 
gemeinschaft fiir die Bereitstellung finanzieller Mit- 
tel und der Universitat Gottingen, die zusammen mit 
der Deutschen Forschungsgemeinschaft die Benut- 
zung einer IBM 650 fir die numerischen Rechnun- 
gen ermoglichte. (Eingegangen am 20. April 1960.) 
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POUR L’ETALONNAGE ABSOLU DE CERTAINS 
MICROPHONES 


par G. B. Mapetia 


Istituto di Fisica Tecnica. Facolta di Ingegneria dell’ Universita di Pisa 


Sommaire 


Les excitateurs 4 condensateur plan offrent une méthode trés simple pour l’étalonnage 
absolu aux basses fréquences des microphones modernes 4 front plan. Les diverses sources 
d’erreurs sont analysées. Une précision de l’ordre de 0,3 dB s’obtient facilement dans 
V’étalonnage de l’ensemble comprenant le microphone, l’amplificateur et l’indicateur, par 
l'emploi de simples appareils de laboratoire. L’erreur est méme plus petite si l’on dispose 
d’un voltmétre haute tension plus précis. 


Zusammenfassung 


Die Anregung mit einer ebenen Hilfselektrode erméglicht ein sehr einfaches Verfahren 
zur absoluten Eichung moderner Flachmikrophone bei tiefen Frequenzen. Die verschiedenen 
Fehlerquellen werden untersucht. Eine MeSgenauigkeit von etwa 0,3 dB laBt sich bei der 
Kichung der gesamten aus Mikrophon, Verstairker und Anzeigeinstrument bestehenden 
Anlage mit einfachen Laborgeraten leicht erzielen. Der Fehler wird sogar noch kleiner, 
wenn ein genaueres Hochspannungsvoltmeter zur Verfiigung steht. 


Summary 


Plane condenser electrostatic actuators offer a very simple means for absolute low fre- 
quency calibration of modern flat-topped microphones. Various error-sources are analyzed. 
In the overall calibration of the microphone, amplifier and indicating instrument, for in- 
stance, the total error is easily reduced to some 0,3 dB, by employing ordinary laboratory 
equipment. By use of more accurate high-voltage meters, a further reduction in error is still 
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ANALYSE DES EXCITATEURS A CONDENSATEUR PLAN 


possible. 


1. Introduction 


Dans une communication précédente [1] on a 
envisagé l’emploi de l’excitateur électrostatique a 
condensateur plan pour l’étalonnage absolu de cer- 
tains types de microphones. En effet, les excitateurs 
électrostatiques usuels, qui offrent peut-étre la 
méthode la plus simple pour l’étalonnage d’un 
microphone, sont employés presque exclusivement 
pour l’étalonnage «relatif» (caractéristique de ré- 
ponse en fonction de la fréquence). Cela s’explique 
du fait que les excitateurs usuels donnent lieu a des 
erreurs méme trés importantes, mais généralement 
indépendantes de la fréquence. 

Il est bien connu que l’ensemble microphone- 
excitateur peut étre étalonné a son tour en tension 
continue, par comparaison avec une pression stati- 
que [2], mais ce procédé — bien qu’il soit peut-étre 
excessif de l’oublier tout a fait — est incontestable- 
ment peu aisé. C’est pourquoi, on mesure générale- 
ment l’efficacité absolue d’un microphone 4a une fré- 
quence donnée par d’autres méthodes (a savoir par 
les méthodes de réciprocité ou par la méthode du 
pistonphone) [3]. Comme, toutefois, ces méthodes 


sont, elles aussi, peu aisées, il s’avéra fort désirable 
daméliorer la précision de l’étalonnage absolu par 
excitateur, tout en gardant, autant que possible, la 
simplicité de cette méthode. 


Etant donné que la source d’erreurs la plus im- 
portante, dans l’emploi des excitateurs, est l’incerti- 
tude dans la mesure de la distance d entre le dia- 
phragme et l’électrode, on choisit avant-tout une 
distance suffisante pour pouvoir la mesurer avec 
une précision relative satisfaisante. I] n’est pas pos- 
sible, toutefois, d’éloigner simplement un excita- 
teur usuel du microphone, car on obtiendrait ainsi 
sur le diaphragme une pression non uniforme. On 
peut bien songer a modifier la forme de l’excitateur, 
pour obtenir une pression uniforme, mais cela 
requiert l’usinage de précision d’une surface compli- 
quée. C’est pourquoi, dans le cas de certains types 
de microphones modernes 4 front plan (voir par 
exemple [4] et [5]), les excitateurs 4 condensateur 
plan offrent la méthode la plus convenable pour 
effectuer l’étalonnage absolu aux basses fréquences, 
avec une précision suffisante pour la plupart des 
applications. 


aan a 
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2. Excitateurs a condensateur plan 


Comme on I’a montré en [1], l’excitateur a con- 
densateur plan se compose d’une électrode plane, dis- 
posée parallélement au diaphragme du microphone, 
a une distance d de quelques millimétres au moins. 


Dans le but d’éviter toute incertitude dans le calcul 


de la distance effective, |’électrode n’a pas de trous 
ou de sillons. Enfin, l’uniformité du champ électro- 
statique sur la surface du diaphragme est assurée 
par un anneau de garde. La Fig. 1 montre schéma- 
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Fig. 1. Disposition schématique d’un excitateur 4 con- 
densateur plan. M microphone, D diaphragme, 
E électrode, A anneau de garde. 


tiquement cette disposition. La valeur de la pression 
instantanée sur le diaphragme, p, en fonction de la 
valeur instantanée de la tension, v, est alors, en 
unités Giorgi: 
ae 

ge a 
é) étant la constante diélectrique de lair. Dans ce 
qui suit, on va analyser les sources d’erreurs les plus 
importantes de cette méthode. 


3. Mesure de la distance d 


La distance entre l’électrode et l’anneau de garde 
peut se mesurer facilement avec toute la précision 
requise. I] est quelque peu plus difficile de verifier 
la complanarité de l’anneau avec le diaphragme sans 
endommager ce dernier. En appuyant soigneuse- 
ment une plaque plane sur le diaphragme, ou bien 
par des méthodes optiques, on reduit toutefois sans 
difficulté erreur a moins de 0,02 mm. Si |’on dis- 
pose d’une sonde électrostatique haute fréquence, 
une précision méme meilleure s’obtient facilement. 


4. Irrégularité du front du microphone 


Quelques microphones modernes ont le front pra- 
tiquement plan et s’adaptent trés bien a l’anneau de 
garde. En d’autres cas, le front n’est pas complete- 
ment plan. On cite par exemple le microphone Bruel 
et Kjer Type 4111, pour lequel, lorsque la surface 
du diaphragme est complanaire avec celle de |’an- 
neau, il reste une rainure profonde quelques dixie- 
mes de millimétres autour du diaphragma. Bien que 
Verreur qui en découle soit trés petite, on peut 
’éviter aisément en posant dans la rainure une ron- 
delle conductrice la remplissant complétement. 
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5. Réaction acoustique 


L’électrode plane n’ayant pas de trous ou de sil- 
lons, la réaction acoustique sur le diaphragme est 
plus importante que dans le cas des excitateurs 
usuels. Bien que le calcul rigoureux ne soit pas 
simple, on peut analyser cet effet par approxima- 
tion, en supposant une propagation par ondes cylin- 
driques entre |’électrode plane et l’anneau de garde, 
et en négligeant l’impédance de radiation a l’aréte 
extérieure (rayon rz). La réactance acoustique de la 
surface cylindrique entourant le diaphragme (rayon 
r,) se calcule alors par la théorie des ondes cylindri- 
ques [6]. La Fig. 2 en donne les valeurs en fonction 
de kr, (k=2af/c; f fréquence; c vitesse de pro- 
pagation) pour deux valeur de r./r, . Par comparai- 
son avec |’impédance acoustique propre du micro- 
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Fig.2. Réactance acoustique normalisée que l’excita- 
teur (rayon rs.) présente au diaphragme (rayon 
m4). k=2au fle; m=2arl dXac/oc ; Xac réac- 
tance acoustique effective; 0 masse spécifique du 
milieu; c vitesse de propagation. 


phone, on peut alors calculer les erreurs qui décou- 
lent de la réaction acoustique [7]... [10]. 


6. Type de l’excitation 


De méme que tout autre excitateur électrostati- 
que, l’excitateur plan peut étre employé soit avec 
polarisation continue, soit sans polarisation. Pour 
des étalonnages de trés grande précision, l’emploi 
avec polarisation donne probablement les meilleurs 
résultats. En effet, pour une distance d donnée, la 
sensibilité du microphone s’obtient, par la premiére 
méthode, en fonction de la tension continue de pola- 
risation et du rapport de deux tensions alternatives 
de méme fréquence, grandeurs que |’on peut mesurer 
avec une trés bonne précision, tandis que, par la 
deuxiéme méthode, on l’obtient en fonction d’une 
tension alternative et du rapport de deux tensions 
alternatives de fréquence différente, grandeurs que 
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l'on mesure avec moins de précision. Une compa- 
raison semblable peut se faire si l’on étalonne en 
pression l’ensemble formé par le microphone, 
Yamplificateur et l’appareil indicateur ou enregi- 
streur. Dans plusieurs cas, toutefois, la deuxieme 
méthode s’avére plus simple, tout en donnant une 
précision suffisante, comme on le verra plus loin. 


7. Distorsion de la tension excitatrice 


Les harmoniques éventuelles de la tension excita- 
trice entrainent une erreur, en conséquence de la 
relation quadratique entre la pression et la tension. 
Cette erreur est toutefois négligeable si l’excitateur 
est employé avec polarisation. En effet, on déduit 
de Eq. (1) que lVerreur relative, E,.), causée par 
exemple par un second harmonique, est, au maxi- 
mum 

Lev 


E,a= =~ => 2 
“$7 (2) 


V, étant la valeur de créte de la tension du second 
harmonique et Vy la tension continue de polarisa- 
tion. L’influence des harmoniques supérieurs est 
méme plus petite. 

Dans l’emploi sans polarisation, ’harmonique le 
plus génant est le troisiéme, qui peut donner, au 
maximum, une erreur relative: 


Laie, (3) 


ou V, et Vs sont les tensions du premier et du troi- 
siéme harmonique, respectivement. On doit donc 
veiller a ce que le taux du troisiéme harmonique ne 
dépasse pas la valeur que l’on caicule a l’aide de 
Eq. (3) en fonction de l’erreur admissible *. 


8. Composante continue de pression 


Tout excitateur électrostatique exerce sur le dia- 
phragme une composante continue négative de pres- 
sion. égale pression alternative, la composante 
continue est minimum si l’excitateur est employé sans 
polarisation (emploi en doublement de fréquence), 
étant égale en ce cas a l’amplitude (valeur de créte) 
de la pression alternative. Les erreurs qui en décou- 
lent dépendent des défauts de linéairité du micro- 
phone a étalonner, et pour un bon microphone sont 
généralement petites. 


1 Les expressions (2) et (3) sont valables pour toute 
autre méthode fondée sur une relation quadratique, a 
savoir pour la méthode du thermophone. 

Les données ci-dessus ne se rapportent évidemment 
pas aux erreurs propres des voltmétres dans la mesure 
de la composante de tension désirée, superposée a 
d’autres composantes. Ces erreurs, différentes pour 
chaque type de voltmétre, sont d’ailleurs communes 4 
toute autre méthode d’étalonnage. 
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9. Raideur négative 


Il est bien connu que, par effet de lexcitateur 
électrostatique, la raideur du diaphragme devient 
plus petite et variable en fonction de la pression 
instantanée. On tire de Eq. (1), en effet, que l’on 
doit escompter une raideur négative, égale au double 
de la force instantanée exercée sur le diaphragme, - 
divisée par la distance d de l’excitateur. Cet effet, 
qui est déja petit pour les excitateurs usuels, est 
négligeable pour les. excitateurs en condensateur 
plan, ot la distance d est plus grande pour les rai- 
sons déja mentionnées. 


10. Essais expérimentaux 


On a employé un excitateur 4 condensateur plan 
pour l’étalonnage des microphones Bruel et Kjer 
Type 4111. Les dimensions de l’excitateur sont mon- 
trées dans la Fig. 3. Par substitution des supports 


| 


N 


SN 


SSTILEL DEL OS i= . 


N 


Fig. 3. Excitateur 4 condensateur plan. B trois bou- 
lons isolés 4 120°, S surface s’appuyant sur 
une saillie réglable du microphone a étalonner, 
R supports isolants interchangeables. 


isolants, on peut choisir diverses distances d, a 
savoir 4,76 mm et 9,52 mm (ces deux valeurs ayant 
été choisies car on disposait de deux calibres de 
précision de 3/16” et de 3/8” respectivement). La 
distance dy de lélectrode au diaphragme étant 
mesurée a 0,03 mm prés, on peut calculer la pres- 
sion électrostatique en fonction de la tension excita- 
trice a 0,11 dB prés pour la distance de 4,76 mm, et 
a moins de 0,06dB prés pour la distance de 
9,52 mm. 

L’effet de la réaction acoustique a été vérifié expé- 
rimentalement pour la distance d de 4,76 mm, par 
la disposition de la Fig. 4, o& le microphone M, 
fonctionne en émetteur acoustique, et le microphone 
My, en récepteur. Par comparaison avec un coupleur 
de 19,5 cm?, on a trouvé les valeurs de la réactance 
acoustique, X,p, représentées sur la Fig. 5 par des 
petits cercles”. La courbe en trait plein représente 
les valeurs calculées a V’aide de la Fig. 2. La droite 


2 La distance entre les microphones étant petite en 
comparaison de la longueur d’onde, la tension aux bor- 
nes de My, pour un courant donné en M,, est pro- 
portionnelle a4 Vimpedance acoustique «vue» par les 
microphones, pourvu que ceux-ci aient une admittance 
acoustique négligeable. 
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Lm, 
pee ler) 
suo Dees pacientes 
——— —— 
Bg B2 
nee 
~ Fig.4. Schéma de la disposition employée pour la 
mesure de la réactance acoustique de l’excita- 
teur. B,, By plates semblables 4 l’anneau de 
garde de l’excitateur. M;,;, Ms, microphones 
fonctionnants, respectivement, en émetteur et 
en récepteur. 


Réactances acoustiques X,(unités Giorgi) ——> 


6 
4 
100 200 400 600 Hz 1000 
Fréquence —> 
Fig.5. Diagramme inférieur: réactance acoustique, 


Xag, de l’excitateur. Courbe en trait plein: 
valeurs calculées a l'aide de la Fig. 2. Petits 
cercles: valeurs mesurées a l’aide de la disposi- 
tion dela Fig.4. 

Diagramme supérieur: réactance acoustique, 
Xam (changée de signe), du. «volume équiva- 
lent» du microphone. 


descendente dans la partie supérieure du diagramme 
représente la réactance acoustique, X,,,, du micro- 
phone (changée de signe). On voit que, dans la 
gamme de fréquence explorée, |’impédance acousti- 
que opposée par l’air est au moins 350 fois plus 
petite que celle du microphone, l’erreur qui en dé- 
coule étant donc négligeable. 

La raideur négative, calculée comme on Il’a déja 
vu, méme pour la distance d la plus petite et pour 
les pressions électrostatiques les plus élevées est, en 


3 La réactance X,y a été calculée d’aprés le volume 
équivalent du microphone. Du fait que la caractéristi- 
que de réponse du microphone ne monte pas au voisi- 
nage de la fréquence de résonance, on déduit que le dia- 
phragme est suffisamment amorti pour que son impé- 
dance acoustique, Z,,,, ne soit jamais plus petite, en 
valeur absolue, que | X44]. 
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valeur absolue, de l’ordre de 10~° fois la raideur 
propre du diaphragme, et on peut donc la négliger. 
Le but étant d’obtenir létalonnage absolu 4 une 
fréquence quelconque (la courbe de réponse étant 
donnée par un excitateur usuel), on a choisi pour la 
tension excitatrice, employée sans polarisation, la 
fréquence du secteur, ce qui correspond, pour la 
pression, a une fréquence d’environ 100 Hz. On a 
ainsi, tout simplement, a l’aide d’un transformateur, 


une tension de l’ordre de 4000 V ou de 8000 V, 


_donnant, pour les deux distances respectivement, 


une pression de l’ordre de 2,2 newton/m?, (env. 
101 dB réf. 2:10~5 newton/m?). 

La tension excitatrice a été mesurée, a l’aide d’un 
diviseur de tension a résistances, 4 1% prés. L’erreur 
maximum qui en découle pour la pression est donc 
de 2%, soit 0,17 dB. La composante a 150 Hz était 
inférieure a 0,1%, l’effet en étant donc inférieur a 
0,02 dB. On a donc, méme en employant la valeur 
la plus petite de la distance d, une erreur totale 
maximum d’environ 0,3 dB. 

On doit remarquer qu’une précision de |’ordre 
indiqué ci-dessus, suffisante en bien des cas, s’ob- 
tient par l’emploi de moyens trés simples. De |’ana- 
lyse précédente on déduit, toutefois, qu’une amélio- 
ration peut s obtenir sans difficulté, notamment par 
une mesure plus précise de la tension excitatrice, a 
Vaide par exemple d’un transformateur étalon et 


d’un voltméetre électrodynamique. 
(Regu le 22 Septembre 1960.) 


Annexe 


i. Galcul de Vimpédance acoustique de l’excita- 
teur. On pose, pour les ondes cylindriques sinusoi- 
dales entre l’électrode et l’anneau de garde 


P=AJ)(kr)+BN (kr), (4) 
ae eee Yd 

U= neers 5 
ee dr (5) 


Ou: P pression; U vitesse; Jy, No fonctions de 
BrsseL d’ordre zero, de premiére et de seconde 
espéce, respectivement; @c résistance acoustique 
intrinséque du milieu. Autres symboles comme dans 
le texte. Ayant supposé nulle l’impédance de radia- 
tion a l’aréte extérieure (r =r,), on tire de Eq. (4) 


Afp No (kre) (6) 


J(k rg) 
et par substitution en Eq. (4) et Eq. (5), on calcule 
P/U en ry. Enfin, en divisant par 2xr,d, on ob- 
tient ’impédance acoustique cherchée. Le calcul de 
la dérivée dP/dr se fait le plus aisément par des 
tables des fonctions de Brssex de premier ordre, 
J, et N, (voir par exemple [11]). 

it. Calcul de V'effet des harmoniques de la tension 
excitatrice. Si la tension excitatrice est 


v=V)+V,sinwt+V,cos(2mt+¢9), (7) 
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on a d’ou 
1 2 
Be lex | 3k 0) IP dake Sole el eRe de nernee as at, GD. eae v kan ane 
? dx 2 (dg +2)* do t+2 dy+ax 
+2V )V,snot—V,V,zsin(wt+@)+ (8) (13) 


d’ot |’erreur maximum Eq. (2) sig=0 o0uqm=n. 
Si la tension excitatrice est 


v=V,snot+V,sin(3 wt+¢@) (9) 


En multipliant par l’aire du diaphragme, on a 
enfin la relation mentionnée au paragraphe 9. 


[1] 
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BETRACHTUNGEN ZUR FRAGE DER OPTIMALFREQUENZ 
BEI WASSERSCHALLORTUNGSVERFAHREN 


von H. TurepE 


Atlas-Werke AG, Bremen 


Zusammenfassung 


An Hand der frequenzabhangigen Parameter bei der Schallausbreitung im Meerwasser 
wird fiir Ortungsverfahren untersucht, welche Optimalfrequenzen sich zur Erzielung einer 
bestimmten Reichweite r bei ungestérter Schallausbreitung fiir konstante Strahlerflache 
ergeben. AufSerdem wird der Zusammenhang zwischen Reichweite und erforderlicher Strahl- 
flache bei Optimalfrequenz errechnet und dargestellt. 


Summary 


Assuming undisturbed sound propagation and constant radiation-emissive surface area 
the optimum frequencies are investigated, for sonar systems, which lead to a given range r, 
taking into account the frequency dependent parameters for sound propagation in sea water. 
Furthermore the relation between range and the necessary minimum radiation-emissive sur- 
face area at optimum frequency is calculated. 


Sommaire 


En uitilisant les paramétres, variables avec la fréquence, qui interviennent dans la pro- 
pagation du son dans l’eau de mer, on a fait des recherches de détection pour déterminer les 
fréquences optimum qui permettent d’obtenir une certaine portée de la propagation sonore 
non déformée quand la surface de l’émetteur est constante. On a ensuite expliqué et cal- 
culé la relation qui existe entre la portée et la surface minimum de |’émetteur pour la fre- 
quence optimum. 
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1. Einleitung 


Wasserschallortungsverfahren arbeiten mit Schall- 
impulsen im Frequenzbereich von einigen kHz bis 
uber 100 kHz. 

Um Echos von Objekten zu erkennen, mu der 
Echopegel E aus dem Storpegel erkennbar sein. 

Im folgenden wird untersucht, welchen Einflu8 
die die Schallausbreitung frequenzabhangig beein- 
flussenden Grofien auf das Zustandekommen des 
Echopegels haben und welche Frequenz zur Erlan- 
gung einer bestimmten Reichweite optimal geeignet 
ist. 


2. Optimalfrequenz bei Reflexion 
an Einzelobjekten 


Wir betrachten den Fall ungestérter Schallaus- 
breitung im freien Wasser ohne Begrenzung durch 
Wande mit Reflexion an Einzelobjekten und an 
spiegelnden Flachen. Der erste Fall ist bei der 
Ortung von Objekten, wie Fischschwarmen, Walen, 
Schiffen und ahnlichen im tiefen Wasser gegeben, 
der zweite Fall bei der Echolotung, also der Ab- 
standsbestimmung des Meeresbodens vom Schiff her. 

Wir definieren folgende Grofen: 


S Sendepegel in dB rel 1 wbar, 

E Empfangspegel in dB rel 1 ubar, 
P,,_akustische Leistung in W, 

D BiindelungsmaB in dB, 

s Strahlflache des Wandlers in m?, 

A4-  Wellenlange im Wasser in m, 

Te Zielma8 in dB, 

f Frequenz, 

Af  Bandbreite des benutzten Empfangssystems, 


H(r) Ausbreitungsfunktion, 
r Entfernung in km. 

Zwischen den Groen bestehen folgende Zusam- 
menhange: 

S=10lg P+ D-+ 70,9 dB rel 1 wbar , 

D=10lg4x +5 dB fiir Strahlflachen s> 72. 

Im Falle der Reflexion an Einzelobjekten ist 
E=S+T—2H(r) (1) 

H(r) =20lgr+ar+60. 

a ist die Absorptionskonstante des Seewassers. Fiir 

eine analytische Betrachtung mu8 a, das durch 

Experimente bei verschiedenen Frequenzen gemes- 

sen wurde, als frequenzabhangige Funktion dar- 

gestellt werden. 

Eine solche, die die Verhaltnisse im Frequenz- 
bereich bis rund 100 kHz in guter Annaherung wie- 
dergibt, wurde von SuerHy und Ha.tey [1] auf- 
gestellt. Sie fanden a=0,036 f* dB/km, wobei f in 
kHz einzusetzen ist. 

An der Grenze der Empfangsempfindlichkeit mu8 
nun der Echopegel E im Stérfeld, das durch einen 


mit 
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mittleren Spektralpegel Ni, gekennzeichnet sei, er- 
kennbar sein. Also mu8 


E=Nj,+10lgAf+M—D 


sein, (2) 


wobei M das Erkennbarkeitsma}} (recognition in- 
dex) in dB ist. Aus den Gleichungen (1) und (2) 
folgt 


2H(r) =S+T+D-10lg Af—-M_N,,= 
=2(20lgr+ar) +120. (3) 


Wir untersuchen den Fall, wie sich bei konstanter 
Strahlerflache, die zum Senden und Empfang dient, 
die einzelnen Groen in Abhangigkeit von der Fre- 
quenz andern. Wir setzen dabei voraus, da der 
Wandler bei allen Frequenzen mit gleicher Intensitat 
sendeseitig erregt werde. In Gl. (3) andern sich 
frequenzabhangig die einzelnen Summanden wie 
folgt: 


§ 
oie 104A, + 70,9 


im mittleren Summanden mit 20 lg f. 

Uber das Zielma8 T lassen sich keine allgemein 
gultigen Aussagen beziiglich seiner Frequenzabhan- 
gigkeit machen. Aus diesem Grunde wird ein ge- 
gebenenfalls vorhandener Frequenzgang eines indi- 
viduellen Zieles auSer Betracht gelassen. 

Der Wert D=10lg4a a andert sich mit 20 lg f. 

Die Bandbreite des Empfangssystems Af wird 
man mit Riicksicht auf gute Storbefreiung so schmal 
wie zulassig wahlen. Infolge der gegenseitigen 
Bewegung von Schiff und Ziel bestimmt in der 
Praxis der Dopplereffekt die Bandbreite. Bei einer 
maximal vorkommenden Relativgeschwindigkeit v, 
ist die Frequenzabweichung infolge des Doppler- 
effektes Af = 20m fo. Somit andert sich der Band- 

c 
breiteneinflu8 mit 


—10lgf. 


Der Summand M, das Erkennbarkeitsma8 (recog- 
nition index), ist frequenzunabhangig. Dagegen 
andert sich der mittlere Spektralpegel des im See- 
raum vorhandenen Storgerausches in Abhangigkeit 
von der Frequenz. Messungen ergeben im allgemei- 
nen Abhangigkeiten mit etwa — 20 lg f [2]. 

Die Frequenzabhangigkeit der Gl. (3) laBt sich 
somit folgendermafen darstellen: 


40 lg r + 0,072 f!*r =50 lg f+C 
oder F(r, f) =40lgr+0,072 f* r—50 lg f=C. 


(3 a) 


Um zu ermitteln, fiir welche Werte f die Werte r ein 
Maximum werden, bilden wir 


dfs. <-OF.. <Af OF dr 
ae i 
df cc} Or df 
und setzen als Maximalbedingung dr/df=0. 
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Da die Funktion F(r, f) eine Konstante ist, mu 
das totale Differential von F verschwinden (dF =0), 
also 


dF OF 50 


fear at 2 0,072 f*r=0 werden. 
df Of 


finlO 2 


Daraus folgt mit In 10 = 2,3 
foot = (200/r)"* (4) 


als optimale Frequenz zur Erzielung der Reich- 
weite r im Falle der Reflexion an Einzelobjekten. 
f ist dabei in kHz, r in km einzusetzen. 


3. Optimalfrequenz bei spiegelnder Reflexion 


Fir den Fall der Reflexion an einer spiegelnden 
Fliche -bekommt die Ausbreitungsfunktion H (r’) 
die Form 


H,(r’) =20lg2r’+2ar’ +60 (5) 


und der Echopegel E wird zu 
E=S + Teer ) 


wobei 7’; das Reflexionsmaf} des Spiegels, also ins- 
besondere des Meeresbodens, ist. Auch in diesem 
Falle mu8 E=N;,+10lgAf+M—D sein, so daB 
an Stelle von Gl. (3) jetzt folgendes tritt: 


60+ 20]lg¢27r°+2ar’ = (6) 
=S$+T,+D—10lgAf—M—N;,. 


Die Frequenzabhangigkeit der einzelnen Summan- 
den in Gl. (6) ist dieselbe wie oben mit Ausnahme 
des Bandbreiteneinflusses beim Echolot. Beim Echo- 
lot, also bei Vertikallotung, tritt der Dopplereffekt 
nicht in Erscheinung; infolgedessen bestimmen 
andere, nicht frequenzabhangige Parameter die 
Wahl der giinstigsten Bandbreite. Im Falle des Echo- 
lotes 1aBt sich die Frequenzabhingigkeit der Aus- 
breitungsfunktion daher ansetzen zu 


20lg27r’'+2ar’'=C+60lgf. 


Nach analoger Rechnung erhalt man in diesem Fall 
als Optimalfrequenz 


fort = (240/r’)** 2 (7) 


wobei r’ die Tiefe in km, f die Frequenz in kHz ist. 


4, Mindest erforderliche Strahlflachen 
bei bestimmten Reichweiten 


Weiter wird untersucht, welche Strahlflache s in 
m? bei der Optimalfrequenz f,p; zur Erzielung einer 
bestimmten Reichweite r in km bei ungestoérter 
Schallausbreitung erforderlich ist. Fiir den Fall der 
Reflexion am Einzelobjekt schreiben wir dafir 
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Gl. (3) in folgender Form: 
2H(r) =S+T+D—101g (Af) 15 cH — M —Nisys ica, 
+10]¢ f/15 (8) 


und beziehen uns auf eine ,,Standardfrequenz* 
15 kHz, auf die Bandbreite des Empfangskanals 
und isotroper St6rpegel umzurechnen waren, und 
beriicksichtigen davon abweichende Frequenzen im © 
letzten Glied. Wir fassen ferner die Glieder, die kei- 
nen Einflu8 auf die Strahlflache und Entfernung 


haben, zusammen und setzen 


By=T — 101g (Af) i510 —M—Nisisigr» (9) 


Damit wird Gl. (8) zu 
2H(r) =101lgP..+D+714+D+10lgf— 
me bi bas ey 3 Pe 
Erregen wir den Wandler sendeseitig mit einer 


Intensitat J W/m? und beriicksichtigen, daf bei einer 
Schallgeschwindigkeit von 1,5 km/s 


2 
D=1lg4a sf = 20 lg f-+101gs+7,5 


ist, dann wird 
2H(r) =10lgs+10lg7+ 
+2(20lg/+10lgs+7,5) +10lgf+ 
+By+59,2= 
=301gs+50lgf+B)+10lg/+74,2. 


(10) 


Die Ausbreitungsfunktion H(r) nimmt, wenn wir r 
in km und f in kHz einsetzen, folgende Form an: 


H(r) =20lgr+60+0,036rf". 


Gl. (10) wird damit zu 
40 lgr +120 + 0,072 r fi? =30 lg s+501g f+ By + 


+10lg7+74,2, (11) 
oder 
301g s =40 lg r+ 45,8 + 0,072 r f' 
— (50lgf+B)+10lg/). (.L2) 


Wir fiihren jetzt f=f.»,=(200/r)** ein und er- 
halten 
30 lgs=40lgr+45,8+ 14,4 
— ($50 lg 200 —2 50 lgr+By)+10lg/) 
= 73,3 lg r—By—10lgJ—16,8, 


oder 


Ig s=2,441gr— (B,/30) —4lgI—0,56. (13) 


Dabei sind s in m?, 7 in km und J in W/m? ein- 
zusetzen. Bild 1 zeigt den Zusammenhang der Gro- 
Ben fopt, r und s bei Schallintensitiéten von 
3000 W/m?, die bei Ortungsanlagen iiblich sind. 
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Bild 1. Zusammenhang zwischen Reichweite r, Optimal- 
frequenz fopt und mindest erforderlicher Strahl- 
fliche fiir verschiedene Werte By bei Optimal- 
frequenz im Falle der Reflexion an Einzel- 
objekten ([=3000 W/m?). 


Fur den Fall der Echolotung nimmt die Ausbrei- 
tungsfunktion (6) die Form an: 


H(r’) =20]g27r’'+2ar’ +60 
=S+T,+D—1lgAf—M—Nisys + 201g 

Hieraus wird, mit By =T,—10lg Af—M — Nig; , 

unter Beriicksichtigung der Werte fiir P,,=SJ und 

D wie oben 

H(r’) =10lgs+10lgJ+ 
+2(20lgf+10lgs+7,5)+20lgf (14) 
+ By +59,2 : 

und mit f = fopp = (240/r’)**. 

H(r’) =301lgs+10lg/+ 40 lg 240 — 40 lg 7’ 
+B, + 74,2. (15) 


Setzen wir fiir die linke Seite 201g 2r’+2ar’ +60 


H. THIEDE: DIE OPTIMALFREQUENZ BEI WASSERSCHALLORTUNGSVERFAHREN 


i) 
S| 
lop) 
Strahlflache — 


Frequent —> 


1 — 0,001 
O1 O02 0406 10 20 406010 20 4Okm 
Tiefe pr’ 


Bild 2. 


Zusammenhang zwischen Tiefe r’, Optimalfre- 
quenz fopt und mindest erforderlicher Strahl- 
flache fiir verschiedene Werte By bei Optimal- 
frequenz beim Echolot (J=3000 W/m?). 


mit a = 0,036 f'? , so wird, wenn wir dort 
f= fopt = (240/r’)* 
setzen, diese zu 
20 lg 2r’ + 0,072-240 + 60 = 20 Igr’ + 83,3. 
Damit wird 
lgs=21g7’ — (By /30) —4 lg J—1,53 
(siehe Bild 2). 


Auch hier ist s in m®, r’ in km und J in W/m? 


einzusetzen. 


(17) 


(Eingegangen am 8. September 1960.) 
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MEASUREMENTS OF ACOUSTIC ATTENUATION 
AND: VELOCITY IN-SAND 


by A. N. Hunter, R. Lecce and E. Marsuxawa 
Physics Dept., University of Leicester, England 


Summary 


Although there is a considerable amount of published work on the.velocity of sound in 
granular materials such as sand, little has been published on the acoustic absorption in 
such media. Experiments show that, whilst little difficulty is experienced in estimating 
velocity using quite small samples, the measurement of absorption is generally difficult 
and is best carried out in bulk media. In the author’s laboratory an attempt has been made 
to measure sound absorption in small samples of granular materials under controlled con- 
ditions. The work has shown that there is a lower limit to the size of the sample for which 
reliable attenuation measurements can be obtained. In the case of sand the upper frequency 
limit of 73 ke/s has been set by attenuation whilst the lower 7 ke/s limit was due to the 
size of the sample. No velocity dispersion or variation of loss factor was found over the 
frequency range investigated. 


Sommaire 


Bien qu’on ait publié un grand nombre de travaux sur la vitesse du son dans les maté- 
riaux se présentant sous forme de grains on a écrit peu de choses sur l’absorption produite 
par de tels milieux. L’expérience montre que, tandis qu’il y a peu de difficultés expérimen- 
tales a déterminer la vitesse de propagation dans le cas d’échantillons de dimensions trés 
petites, la mesure de l’absorption est en général difficile et s’obtient mieux dans le cas de 
milieux formés de masses importantes. On a essayé dans le laboratoire de l’auteur de 
mesurer l’absorption du son sur des échantillons de petites dimensions de matériaux en 
forme de grains et ceci dans des conditions bien déterminées. Ce travail a montré qu’il y 
a une limite inférieure des dimensions de l’échantillon au-dessous de laquelle on ne peut 
plus faire des mesures valables de l’atténuation. Dans le cas de sable la fréquence limite 
supérieure obtenue pour l’atténuation est de 73 kc/s tandis que l’inférieure est de 7 ke/s 
étant donnée la dimension du sable. On n’a pas trouvé de dispersion de la vitesse ni de 
variation du facteur de perte dans la gamme de fréquences examinée. 


Zusammenfassung 


Obgleich eine beachtliche Zahl von Arbeiten iiber die Schallgeschwindigkeit in kornigen 
Materialien, wie z. B. Sand, bekannt geworden sind, ist wenig iiber die Schallabsorption in 
solchen Medien verdffentlicht worden. Wie Experimente zeigen, ist die Bestimmung der 
Schallgeschwindigkeit auch mit sehr kleinen Proben ohne gréSere Schwierigkeiten méglich; 
die Messung der Schallabsorption ist gewohnlich schwieriger und wird am besten im aus- 
gedehnten Medium durchgefiihrt. Im Laboratorium des Verfassers wurde der Versuch ge- 
macht, unter genau bekannten Bedingungen die Schallabsorption in kleinen Proben von 
k6érnigem Material zu messen. Es hat sich dabei herausgestellt, das es eine untere Grenze 
fiir die GroBe der Probe gibt, bei der noch zuverlissige Diampfungsmessungen durchgefiihrt 
werden kénnen. Bei Proben aus Sand wurde die obere Frequenzgrenze von 73 kHz durch 
die Dampfung und die untere Grenze bei 7 kHz durch die Probenabmessungen bestimmt. 
Im untersuchten Frequenzbereich wurde keine Dispersion der Schallgeschwindigkeit und 
keine Anderung des Verlustfaktors gefunden. 


1. Preliminary measurements 


Preliminary measurements were made using ver- 
tical propagation in a box of sand 35 cm square in 
section in order to obtain orders of magnitude for 
the velocity and absorption. The method used was 
a simple one consisting in placing a 20 kc/s magne- 
tostriction transducer, having a surface area of 
about 600 cm”, on a layer of sand in the bottom of 
the box. The active face of the transducer was hori- 
zontal and directed upwards; a layer of sand was 


spread over the transducer and levelled with a 
straight edge. A second similar transducer was 
placed on top of the sand pointing downwards and 
was always adjusted until a maximum signal was 
obtained. Measurements were made by recording 
the setting of a calibrated attenuator in the amplifier 
channel which gave a standard signal output. The 
upper transducer was then removed and the sand 
thickness determined with a wire probe. The experi- 
ment was repeated for different thicknesses of sand 
and a graph was drawn to show the relationship 
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between attenuator setting and thickness of sand 
between the transducers. The slope of this linear 
graph gave the absorption in the sand and velocity 
measurements were made concurrently by pulsing 
the source and measuring the delay time with a cal- 
ibrated time base on a cathode ray oscilloscope. 
Absorption measurements made with both c.w. and 
pulsed carrier gave similar results. Fig. 1 shows 
typical results for the attenuation; comparison with 
the velocity measurements given in Fig. 2 shows the 
greater amount of scatter which occurs with the 


‘absorption measurements. | 
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Fig. 1. Attenuation against sand thickness. 
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Fig. 2. Thickness of sand against pulse transit time. 


Measurements were taken in dry sand and moist 
sand. The velocity in dry sand varied between 275 
and 340 m/s and in moist sand a value of about 
234 m/s was obtained. The values of attenuation 
obtained are summarised in Table I. 


Table I. 
Attenuation | 
Material at 20 ke/s 
dB/em 


Dry sand 1.2 to 2.2 
Moist sand 1.5 to 1.6 
2. Other techniques for measurement 


A number of possible methods were investigated 
in the search for an accurate method suitable for 
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small samples. Unfortunately all had serious draw- 
backs and the best results were eventually obtained 
by a return to a modified version of the experiment . 
outlined above. The difficulties are due mainly to 
two causes: 

i. Inhomogeneities in the sand, arising from very 
small differences of packing, could result in a de- 
structive interference pattern over the surface of the 
transducer. This made the received signal extremely 
sensitive to transducer position; in some cases no 
signal at all could be received over quite small 
paths. 

ii. Variations in packing of the sand against the 
boundary of the container or transducer could re- 
sult in variable mismatch losses which masked the 
true attenuation. 

To minimise these errors it was found that exper- 
iments must be performed with a minimum pressure 
of 4 lb/sq.in on the transducer face (or container 
boundary) and that the transducers should not be 
smaller that about 100 grain diameters across the 
face. In agreement with measurements made else- 
where (WyLLiz, Grecory and Garpner, 1956 [1]) 
it was further found that a constant value for the 
velocity could only be obtained for paths in excess 
of about 100 grain diameters from the transducer 
face. Close to the transducer the measured velocity 
was high and extrapolated to that for bulk silica at 
the face itself. 

The various methods considered are summarised 
below: 


2.1. Reverberation method 


Unfortunately the high attenuation of the order 
of 2 dB/cm taken together with a velocity of about 
250 m/s results in a reverberation time of only 
0.2 milliseconds. At 20kc/s the carrier period is 
0.05 milliseconds so that there are only a few cycles 
in the decay. Furthermore, in order to establish a 
resonant mode, the high attenuation requires a 
sample so small as to be impractical. 


2.2. Waveguide method 


This method seemed promising in theory and 
some time was spent in developing it. The apparatus 
used consisted of a section of electrical conduit 
4 inches square used as a horizontal trough filled 
with water. This trough was lined with sponge rub- 
ber to give a pressure release wall and a second 
lining of polythene sheet to prevent saturation of 
the rubber. The sample was placed in a cell of thin 
sheet metal in the guide and the water level in the 
guide was adjusted to give the required cut-off fre- 
quency for operation in the lowest (1,1) mode. 
Theoretically the attenuation in the sample together 
with the mismatch loss at the walls can be obtained 
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from the standing wave pattern before the sample 
and the signal strength after it. A small barium 
titanate probe was used on a carriage sliding par- 
allel to the axis of the guide and a standing wave 
ratio of about 40 dB could be obtained with an air 
cell as a reflector. 

Although this guide showed promise with solid 
materials and allowed very high reflection coeffi- 
cients to be measured it was abandoned for the 
measurements in sand. This was because of the 
variations of velocity across the cell due to changes 
in packing; tapping the cell could result in a change 
in transmission of 20dB or more. A new type of 
guide was designed using vertical propagation in a 
rectangular tube having thick walls (1.5cm) and 
lined with pressure release material. A thin metal 
diaphragm was used to separate the section con- 
taining sand at the top of the guide from the water 
and the lower part of the guide was completely 
filled with water. Two fixed probes were used and 
the frequency was varied in order to measure the 
standing wave ratio. Fortunately the pattern of 
maxima and minima in a guide is such that, provid- 
ing that the frequency is adjusted to give a mini- 
mum at one probe, the spacing of the probes is not 
critical. The purpose of the vertical arrangement 
was to make possible more homogeneous thrusts by 
weighting a plunger resting on the horizontal layer 
of sand. Also with the guide completely filled with 
water much higher pressures could be used without 
distorting the walls of the container and thus the 
inhomogenuities were much reduced. The guide is 
of particular value for measurement of mismatch 
between water and sand; unfortunately the diffi- 
culties of finding a source transducer of sufficient 
bandwidth and of maintaining the pressure release 
nature of the wall lining under high pressures pre- 
vented further development of the method. 


2.3. Loading on a transducer 


This method is based on the changes in electrical 
impedance to be expected when a transducer is 
placed in contact with different acoustic media. The 
large and variable mismatch between the sample and 
the transducer face, together with the small reactive 
loading of materials such as sand (ratio of attenua- 
tion constant to phase constant about 0.03) ruled 
out this method for absorption methods although it 
may prove useful for mismatch measurements at 
higher pressures. It is to be expected that transducer 
dimensions greater than 100 grain diameters would 
be required. 


2.4, An array of transducers 


In this method four cylindrical transducers of 
Barium Titanate 6 inches long and each about 1 inch 


ACUSTICA 
Vol. 11 (1961) 


in diameter were arranged in a frame so as to form 
the four parallel edges of a rectangular parallel- 
piped. The transducers were radially polarised and 
of 1/8 inch wall thickness; they were buried in the 
sand and measurements of transmission were taken 
between different pairs. From these measurements 
the attenuation in the sand and the mismatch losses 
and insertion losses of the transducers could be ob- 
tained. An experiment in dry sand gave an insertion 
loss of about 35 dB per transducer and an attenua- 
tion coefficient of 6 dB/cm in a cubic metre of sand 
under self weight only. The method was not pursued 
because no means was available of applying suffi- 
ciently large thrusts to such a volume of sand; it 
undoubtedly offers promise for field experiments. 


3. The variation of attenuation with frequency 


In view of the difficulties encountered with the 
methods described in section 2. and considering the 
paucity of information on the variation of attenu- 
ation with frequency, it was decided to make further 
measurements with a modification of the methed 
described in section 1, despite the rather large in- 
accuracies. The following results were obtained by 
containing the sand in a 6 inch diameter can lined 
with 1/8 inch sponge rubber to insulate the sand 
acoustically from the walls. This can was stood on 
a steel disc 3 inches in diameter resting on a ma- 
gnetostriction transducer. The barium titanate de- 
tector was 2 inches in diameter and was mounted on 
a holder fitted with a cathode follower head ampli- 
fier. This detector was loaded with a 2 kp weight to 
maintain uniform contact and its position was 
always adjusted to give maximum signal. The 
amount of sand in the container was varied and the 
signal recorded as a function of thickness of sand 
layer. 

The signal level D, corrected for beam spread, at 
a distance x from the source can be written 


D=Kj dead: 


where K; is a constant for a given frequency and 
amplifier setting. Thus 


K;—D=a2z=z2f (a/f). 


which the values of K; have all been chosen to make 
the lines pass through the origin: it is then found 
that the points fit a single line. It follows that the 
Figure 3. shows plots of (K;—D) against xf in 
slope (a/f) of this line is constant (~10~4dB s/ 
cm) and thus a2 is proportional to the frequency. 
We can express the results in terms of the loss factor 
» (defined as the tangent of the angle of loss) by 
means of the equation 


oo SeGray: dB/cm , 
c 
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Fig. 3. Graph showing variation of (K;—D) with (a f) 


(see text). 


in which 7 is assumed small and c is the velocity of 
propagation. Substituting our value of a/f and the 
convenient value 2.73:10+cms ? for the velocity 
(since it is shown below that there is no significant 
dispersion of the velocity with frequency) we then 
have 

y~O0.1. 


This is a convenient working figure which agrees 
well with the value obtained by Scumipr [2] at 
much lower frequencies. The loss factor in moist 
sand was also found to be approximately 0.1. 


4. The variation of velocity with frequency 


Appreciable dispersion is to be expected when the 
wavelength of the sound approaches the grain dia- 
meter. This occurs at a frequency of about 300 kc/s 
in the case of an average sand having a grain dia- 
meter of about 1 mm; such a high frequency is not 


attainable in practice owing to attenuation. In order 


to test for dispersion at frequencies of practical 
interest a large number of measurements of velocity 
were made in three samples of sand. One sample 
was of coarse natural sand containing granules in 
excess of 1.5 mm in diameter and the other samples 
were obtained by sifting sand to obtain a fairly 
uniform grain size. A constant pressure of 1.25 lb/in? 
was used over the frequency range 7 kc/s to 73 kc/s, 
the lowest frequency being set by standing waves 
due to the low attenuation and limited size of the 
sample. The results are summarised in Table II. 
below and provide convincing evidence for the 


Table II. 
Maximum 
Mean deviation from 
velocity | the mean in any 
single experiment 
m/sec G 


en 


Coarse natural sand 368 4 
Mean grain diameter 

0.86 mm. 344 + 
Mean grain diameter 

0.43 mm. 313 1 
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absence of dispersion in the frequency range in- 
vestigated. 

For these measurements continuous waves were 
used. The phase change in the wave after passage 
through a known thickness of sand was determined 
by squaring and differentiating both the input and 
output signals and measuring the time interval 
between the resulting pair of pulses. 


5. The variation of attenuation, velocity 
and loss factor with pressure 


The apparatus was similar to that described in 
section 3. except that two 3 inch barium titanate 
transducers were used for both source and detector 
in order to provide a more compact arrangement. 
At the lower pressures the top transducer was 
weighted but higher pressures were obtained by 
means of a small hand press with a weighted lever. 
A constant frequency of 25 kc/s was used. 

Two techniques were employed for the velocity 
determinations; one consisted of timing pulses for 
different paths in the sand and the other the tech- 
nique described in section 4. When plotted loga- 
rithmically the results agree within the experimen- 
tal error with the empirical formula 


c=AP*, 


in which P was the pressure and the constant A 
varied between 230 and 350 according to the grain 
size of the sample; k was between 1/5 and 1/6. It 
was found convenient to use “mixed” units and a 
typical set of results is given in Figure 4 a. 

Figure 4b shows the attenuation coefficient a as 
a function of pressure for different grain sizes. 


dB/cm 
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Fig. 4. Variation of velocity and attenuation with pres- 


sure. 
Table III. 

Pressure “ae 10- 
lbs/sq.in. dB/s 
0.284 678.4 
0.568 640.3 
0.852 649.8 
1.136 | 640.7 
1.420 | 695.6 
2.840 655.5 
4,260 641.3 
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Table III. gives values of the product ac calculated 
from the smoothed values taken from Figure 4. It 
will be seen that all the values lie within 6% of the 
mean value of 657-10? db/s which corresponds to 
a loss factor of 0.117. 


6. The effect of grain size 


Four samples of sand of different mean grain 
sizes were prepared by sieving sand samples. The 
mean diameters of the grains were measured by 
counting grains with a travelling microscope and 
are given in Table IV. 


Table IV. 


Mean Grain 


Sample No. Diameter 


mm 


0.86 
| 0.75 
0.43 
0.26 


mw be 


Attenuation 


0 
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08 mm 10 
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Fig. 5. Variation of attenuation with grain size. 
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Fig. 6. Variation of velocity with grain size. 


ACUSTICA 
Vol. 11 (1961) 


The results obtained for attenuation against grain 
diameter at various constant pressures are given in 
Fig. 5 these show that the relation is sensibly linear 
within the range investigated. Fig. 6 shows the 
effect of grain size on the velocity and the result of 
combining the two measurements is to give a varia- 
tion in loss factor between 0.08 and 0.12 as the 
grain size changes from 0.26 mm to 0.86 mm. 


7. The influence of the size of the transducer 
and the length of the acoustic path 


Wve, Grecory and Garpner found that the 
measured velocity of sound in granular materials 
decreases with distance from the source for small 
paths. For infinitesmally small paths the value 
obtained for the velocity corresponded with that for 
the bulk solid, a distance of about 100 grain dia- 
meters being necessary before the measured velocity 
becames constant. This effect is confirmed by our 
measurements and typical results are given in Fig. 7. 
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Fig. 7. Variation of transit time with path for different 
grain sizes. 


It has similarly been found by us that reliable 
attenuation measurements cannot be made with 
transducers smaller than about 100 grain diameters 
in linear dimensions. In the case of smaller trans- 
ducers than this insufficient grains are covered and 
the mismatch loss becomes so critically dependant 
on the details of grain contact with the transducer 
face as to make attenuation measurements valueless. 


8. Summary 


This investigation has produced the interesting 
result that the loss factor in dry sand remains 
sensibly constant and of the order of 0.1 over the 
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pressure range of 0.284 to 4.26 lbs/sq. in and fre- 
quency range 7 kc/s to 73 kc/s. There is some evi- 
dence that the loss factor remains sufficiently con- 
stant for all practical purposes outside these ranges 
as well [2]. No evidence for dispersion of the 
velocity has been found over the frequency range 


~ investigated. It is shown that reliable experimental 


results can only be obtained for paths longer than 
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about 100 grain diameters and with transducers 
which are at least 100 grain diameters in linear 
dimensions. (Received 30th July, 1960.) 
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MESSUNGEN IN ULTRASCHALLFELDERN 
HOMER INTENSITAE IN LUFT 


von P. Bocxser, Braunschweig 


Mitteilung aus der Physikalisch-Technischen Bundesanstalt 


Zusammenfassung 


Mit Hilfe eines Resonanzschwingers wurden in Luft bei 30 kHz sinusformige Ultraschall- 
felder von einer Intensitat bis zu 0,25 W/cm? erzeugt. Bei dieser Intensitat und Frequenz 
bildet sich infolge der Nichtlinearitat der Zustandsgleichung bereits innerhalb eines Lauf- 
weges der Welle von 30cm eine vollstandige Aufsteilung der Wellenfront aus; diese Auf- 
steilung wurde analysiert. In der aufgesteilten Welle wurde der Strahlungsdruck gemessen 
und mit der Energiedichte einer harmonischen Welle verglichen. 


Sommaire 


Au moyen d’un oscillateur.a résonance, on a produit dans l’air, a la fréquence de 30 kHz, 
un champ ultrasonore sinusoidal dont l’intensité peut atteindre 0,25 W/cm?. A cette inten- 
sité et pour cette fréquence, l’équation d’état n’est pas linéaire et il en résulte a l’intérieur 
d’un parcours de l’onde de 30 cm un raidissement complet du front de l’onde; ce raidisse- 
ment a été analysé. On a mesuré la pression de radiation dans cette onde raidie et on l’a 
comparée 4 celle d’une onde harmonique. 


Summary 


Using a resonance oscillator, a sinusoidal ultrasonic field has been produced in air at a 


frequency of 30kc/s and having an intensity up to 0.25 W/cm?. At this intensity and fre- 
quency the equation of state is no longer linear and the wave front within a propagation 
distance of 30cm has become completely steepened. This steepness has been analysed. The 
radiation pressure in this steep-fronted wave has been compared with that of a harmonic 


wave. 


1. Einfithrung 


Ein Schallfeld lat sich durch drei Variable cha- 
rakterisieren: den Schalldruck p, die Dichte @ und 


die Teilchengeschwindigkeit vi, Dementsprechend 
wird es durch drei Gleichungen beschrieben; es sind 
dies 
die Zustandsgleichung des Mediums _ 
p={(9), (1) 


die Bewegungsgleichung 


= dp 
k— grad p=0 — 2 
EE had eg ae P (2) 


1 Es sei vorausgesetzt, daB die Vorginge adiabatisch 
verlaufen und da die innere Reibung des Mediums 
vernachlassigt werden kann. 


und die Kontinuitatsgleichung 


—divov= a (3) 


i sei eine auf die Masseneinheit wirkende aufere 
Kraft. 

In der Akustik sind gewohnlich die Schwingungs- 
amplituden so klein, da} man sich in den nichtlinea- 
ren Gleichungen (1) bis (3) auf Groen erster 
Ordnung beschranken kann, Hierzu vernachlassigt 
man die Nichtlinearitat der Zustandsgleichung (1) 
sowie den konvektiven Anteil der Teilchenbeschleu- 
nigung in Gleichung (2). AuBerdem wird in der 
Bewegungsgleichung (2) und in der Kontinuitats- 
gleichung (3) @ durch Qp ersetzt, soweit @ in Pro- 
dukten mit Groen erster Ordnung auftritt. Bei 
groferen Schalldruckamplituden. werden jedoch 
einige Erscheinungen meBbar, die auf der Nicht- 
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linearitat der exakten Schallfeldgleichungen (1) bis 
(3) beruhen und fiir deren Berechnung deshalb die 
ublichen Vernachlassigungen nicht mehr zulassig 
sind. Dieses sind der durch das Medium hervor- 
gerufene Klirrfaktor des Schalldrucks und der 
Schallstrahlungsdruck. 

In jedem Medium mit nicht streng linearer Zu- 
standsgleichung breitet sich eine Welle endlicher 
Amplitude nicht mit konstanter Phasengeschwindig- 
keit aus; dementsprechend andert sich mit zuneh- 
mendem Laufweg der Welle die Wellenform. Die 
Wellenfront einer urspriinglich sinusformigen Welle 
beispielsweise steilt sich mit wachsendem Abstand 
von der Schallquelle immer starker auf, die Energie 
der Welle wandert immer mehr von der Grundfre- 
quenz in die sich ausbildenden hoheren Harmoni- 
schen [1]. Im dampfungsfreien Fall fihrt die Ver- 
formung schlieBlich dazu, da} die Welle unstetig 


‘wird. Sie besitzt dann einen Klirrfaktor von etwa 


AT%. 


Bei nichtlinearer Zustandsgleichung sind die Volu- 
men- und entsprechenden Druckanderungen bei grofe- 
ren Amplituden einander ohnehin nicht proportional. 
Jeder Wechseldruck, der beispielsweise durch sinus- 
formige Volumendnderung hervorgerufen wird, be- 
sitzt deshalb schon von vornherein einen Klirrfaktor. 
Die Rechnung zeigt [2], dafs dieser Klirrfaktor sogar 
bei einem Schalidruckpegel von 160 dB iiber 2-104 
wbar nicht iiber 2% anwachst, daB er also klein gegen- 
uber dem oben genannten Klirrfaktor infolge der Auf- 
steilung der Wellenfront ist. 


Die Schwierigkeit der Losung der nichtlinearen 
Schallfeldgleichungen (1) bis (3) bei der Berech- 
nung des Schallstrahlungsdruckes wird umgangen, 
indem die Gleichungen von dem raumfesten (,,EuLER- 
schen“) Koordinatensystem in ein solches Koordina- 
tensystem transformiert werden, das auf ein be- 
stimmtes sich bewegendes Teilchen bezogen ist. Fir 
die Bewegungsgleichung in diesen ,,LAGRANGE‘- 
Koordinaten lat sich eine Losung finden, sofern 
man sich bei Medien mit nichtlinearer Zustandsglei- 
chung auf kleine Amplituden beschrankt; bei Medien 
mit konstanter Kompressibilitat ist diese Losung 
streng giiltig, Um den Strahlungsdruck jedoch in 
raumfesten Koordinaten angeben zu konnen, mu 
die Lésung von dem Lacrance- in das Eutrersche Ko- 
ordinatensystem zuriicktransformiert werden. Es 
folgt schlieBlich, daB der Strahlungsdruck, den ein 
Strahl von endlichem Querschnitt in einem damp- 
fungsfreien isotropen Medium als Folge der konvek- 
tiven Beschleunigung auf einen Reflektor ausiibt, 
unabhangig von der Zustandsgleichung dem zeit- 
lichen Mittelwert der Energiedichte im Strahl pro- 
portional ist [3]. 

Es wurden fortschreitende Ultraschallwellen hoher 
Intensitat in Luft erzeugt und in den Schallfeldern 
Analysen der Aufsteilung der Wellenfront und Mes- 


sungen des Schallstrahlungsdruckes durchgefuhrt. 
Hierftir wurde das Anwachsen der nichtlinearen Ver- 
zerrungen des Schalldrucks in Luft bei zunehmen- 
dem Abstand von der sinusformig schwingenden 
Strahlerflache und bei wachsender Schwingungs- 
amplitude untersucht; dann wurde der Schallstrah- 
lungsdruck in den gleichen Schallfeldern bestimmt, 
und es wurden die Ergebnisse der Schalldruck- und ° 
Strahlungsdruckmessungen miteinander verglichen. 


~ 2. Die Eecacaae von Ultraschall 
hoher Intensitat in Luft 


Zur Erzeugung hoher Ultraschallintensitaten fiir 
industrielle Zwecke benutzt man im allgemeinen 
mechanoakustische Wandler wie Sirenen oder Pfei- 
fen. Bei diesen Wandlern befindet sich jedoch das 
Schallfeld vielfach innerhalb des erzeugenden Gas- 
stromes, und die abgestrahlte akustische Energie ist 
meistens auf ein weites. Frequenzspektrum verteilt. 
Deshalb wurde zur Erzeugung von Ultraschall hoher 
Intensitat fiir die vorliegenden Untersuchungen nicht 
ein solcher Wandler verwendet, sondern ein mas- 
siver Duraluminzylinder nach Crarr [4] durch 
elektrodynamische Krafte zu Eigenschwingungen 
angeregt (siehe auch [5]). Dazu taucht ein an den 
Zylinder angedrehter KurzschluBring in das magne- 
tische Feld eines Ringspaltmagneten. Der magne- 
tische Wechselflu8 wird durch eine auf den Kern des 
Ringspaltmagneten fest aufgewickelte Erregerwick- 
lung erzeugt. Die longitudinale Eigenfrequenz des 
Zylinders, die fiir die Frequenz des abgestrahlten 
Schalles maBgebend ist, wird durch die Zylinder- 
lange gegeben; diese ist gleich einer halben Schall- 
wellenlange im Duralumin. 

Besonders wichtig ftir die Erzielung eines még- 
lichst groBen Wirkungsgrades ist eine dampfungs- 
arme Lagerung des Schwingers. Curer [6] hat ge- 
zeigt, daB ein einwandfreies Schwingen des Zylin- 
ders nur dann eintritt, wenn sein Durchmesser d zu 
seiner Lange / im Verhiltnis 


d/l=2 B/x 


steht. Hierin ist 6 eine Nullstelle der 1. Ableitung 
der Besselfunktion 1. Ordnung 


Jy (2) =(0. 


(Infolge eines Druckfehlers in der Arbeit von Cra 
wird dort und in darauf basierenden Arbeiten / als 
Nullstelle der Besselfunktion J, (x) selbst bezeich- 
net.) Danach ergibt sich fiir die erste Lésung 


(d/l) ,=1,178. 
Die Abmessungen des Schwingers, der fiir die 


durchgefiihrten Untersuchungen benutzt wurde, be- 
trugen d=125,8mm und /=106,8mm; er war in 
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seinem Schwerpunkt dampfungsarm aufgehangt. Als 
Ringspaltmagnet diente ein fremderregtes Lautspre- 
chermagnetsystem, mit welchem sich eine Luftspalt- 
induktion von etwas mehr als 20 000 Gauf erreichen 
lief. 

In Bild 1 ist in der Umgebung der Resonanzfre- 
quenz f,., der Ubertragungsfaktor B des Resonanz- 
schwingers, das heift das Verhaltnis von dem in 
einem bestimmten Punkt erzeugten Schalldruck zu 
der Spannung an der Erregerwicklung aufgetragen; 
der Ubertragungsfaktor B ist auf den Wert B,-s bei 
der Resonanzfrequenz bezogen. Infolge des Fehlens 
von Zuleitungen an den schwingenden Zylinder und 
durch die dampfungsarme Lagerung ist die Halb- 
wertsbreite dieser Resonanzkurve sehr gering: Sie 
betragt Af=1,36Hz; das bedeutet bei der Reso- 
nanzfrequenz fres = 30 538 Hz eine Gite 


Q = fres/Af = 2,25- 104. 
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Abweichung von der Resonanzfrequenz —-> 


Bild 1. Relativer Ubertragungsfaktor des Resonanz- 
schwingers. 


Diese scharf ausgepragte Resonanz hat Abstimm- 
schwierigkeiten zwischen dem erregenden elektri- 
schen Generator und dem Schwinger zur Folge, weil 
beide ihre Arbeitsfrequenz im Verlaufe der Mefzeit 
ungleichmaBig andern. Deshalb wurde der Schwin- 
ger in einer Rickkopplungsschaltung durch Selbst- 
erregung betrieben. Bild 2 zeigt den Aufbau, wel- 
cher fiir die Erzeugung von Ultraschall hoher Inten- 


Bild 2. Aufbau zur Erzeugung von Ultraschall hoher 
Intensitat in Luft. 
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sitat benutzt wurde. Im Vordergrund befindet sich 
der Resonanzschwinger, hinter ihm ist das in der 
Selbsterregungsschaltung verwendete Mikrophon zu 
erkennen. Uber einen Vorverstarker (in Bild 2 ganz 
links) und ein Terzsieb (rechts neben dem Vorver- 
starker) ist das Mikrophon an einen Regelverstarker 
(unteres Gerat hinter dem Resonanzschwinger) an- 
geschlossen; bei Schwankungen des Spannungs- 
pegels am Eingang dieses Verstarkers um 24 dB 
andert sich der Spannungspegel am Ausgang um 
nicht mehr als 1 dB. Die geregelte Spannung gelangt 
von dem Regelverstérker an den Eingang eines 
100 W-Leistungsverstarkers (Gerat iiber dem Regel- 
verstarker) und von dessen niederohmigem Aus- 
gang schlieBlich an die Erregerwicklung des Reso- 
nanzschwingers. Mit Hilfe des Rohrenvoltmeters 
(in Bild 2 rechts), das parallel zum Ausgang des 
Leistungsverstarkers geschaltet ist, wird die Span- 
nung an der Erregerwicklung gemessen. 

Das Schallfeld vor der Zylinderstirnfliche des 
Resonanzschwingers ist auf den Bildern 3 und 4 
wiedergegeben. Es wurde mit Hilfe einer von Herrn 
Dipl.-Ing. H.-J. Scuroeper im Laboratorium fiir 
Klanganalyse und Schallaufzeichnung der Physika- 
lisch-Technischen Bundesanstalt entwickelten Appa- 
ratur sichtbar gemacht”. Bild3 zeigt die Kur- 
ven gleicher Schalldruckamplitude. Die Differenz des 


Oe all : 
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Bild 3. Schallfeld vor der Zylinderstirnflache des Reso- 


nanzschwingers. Kurven gleicher Schalldruck- 
amplitude. (Der Abstand zwischen zwei benach- 
barten Kurven entspricht einer Differenz des 
Schalldruckpegels von 5 dB.) 


2 Der Verfasser dankt Herrn Dipl.-Ing. H.-J. Scuror- 
per fiir die Anfertigung dieser Aufnahmen sowie fiir 
die Erlaubnis, die Bilder bereits vor ihrer Verdffent- 
lichung in einer Arbeit tiber die Apparatur zur Sicht- 
barmachung von Schallfeldern publizieren zu diirfen. 
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Schalldruckpegels von Kurve zu Kurve betragt 5 dB. 
In der Nahe der fast nach Art einer kreisformigen 
Kolbenmembran schwingenden Zylinderstirnflache 
sind gut die Nahfeld-Interferenzen und die Neben- 
maxima der Richtcharakteristik zu erkennen (ver- 
gleiche [7]); auBerdem zeigt Bild 3 die scharfe 
Biindelung des Hauptstrahls. Aus den Kurven glei- 
cher Phase auf Bild 4 ist ersichtlich, daB die Wellen- 
fronten in der Nahe der Zylinderstirnflache eben 
verlaufen und da die benachbarten Nebenmaxima 
untereinander jeweils eine Phasendifferenz von 180° 
besitzen. 


Bild 4. Schallfeld vor der Zylinderstirnflache des Reso- 
nanzschwingers. Kurven gleicher Phase. (Der 
Abstand zwischen zwei benachbarten Kurven 
entspricht einer Phasendifferenz von 360°.) 


3. Die Aufsteilung der Wellenfront in Luft 


Nach den ersten allgemeinen Arbeiten von Riz- 
MANN [8] und Hueconror [9] tiber die Theorie der 
Gasschwingungen bei grofen Amplituden wurden 
von Fusini-Gurron [10], [11] und Prrimem [1], 
[12] die Untersuchungen zur Theorie der Wellen 
groBer Schwingungsweite fortgesetzt. Daraus geht 
hervor, daB eine zum Zeitpunkt ¢=0 rein harmoni- 
sche Schallwelle endlicher Amplitude zu allen darauf- 
folgenden Zeitpunkten ¢>0 nicht mehr harmonisch 
ist. Unter Vernachlassigung von kleinen Gliedern 
hoherer Ordnung ergibt sich fiir den Weg x, einer 
urspriinglich harmonischen ebenen dampfungsfreien 
Welle im homogenen Gasraum bis zum Eintritt der 
ersten Unstetigkeit der Wert 

dep ta (4) 
Cy Schallgeschwindigkeit im ungestérten Medium 


(bei kleinen Schalldruckamplituden) , 
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x  Verhaltnis der spezifischen Warmen in dem be- 
treffenden Medium, 

w Kreisfrequenz der harmonischen Schwingung, 

vo Geschwindigkeitsamplitude der zur Zeit t=0 
erzeugten Welle. 

Die Beziehung (4) zeigt, dafi der Weg, innerhalb 
dessen die Wellenform stetig bleibt, mit zunehmen- 
der Frequenz und zunehmender Geschwindigkeits- 
amplitude der urspriinglichen harmonischen Schall- 
welle abnimmt: Im_Ultraschallbereich ist die Auf- 
steilung der-Wellenfront schon in der Nahe eines 
Senders groferer Leistung zu beobachten. 


3.1. Mebanordnung 


Die Messungen wurden in einer reflexionsarmen 
Mefikammer von den inneren Abmessungen 105 cm 
x 78cm x 60cm durchgefihrt, deren Wande in- 
folge einer Auskleidung mit Sillanplatten und einer 
Baumwollwatteschicht im interessierenden Frequenz- 
bereich fast unabhangig vom Einfallswinkel einen 
Amplitudenreflexionsfaktor von etwa 3% besafen. 
Wegen der hohen Ultraschallintensitaten sollte durch 
die Wande auch eine moglichst groBe Schallpegeldif- 
ferenz zwischen dem Kammerinneren und dem 
Auf<enraum bewirkt werden; sie betrug bei Ultra- 
schallfrequenzen etwa 90 dB. 

Zu den Schalldruckmessungen wurde hauptsach- 
lich das piezoelektrische Mikrophon M-213 der 
Firma Massa Laboratories Inc., Hingham, Mass. 
verwendet. Es enthielt einen ADP-Kristall, dessen 
Resonanzfrequenz bei 135 kHz lag; infolge seiner 
sehr hohen akustischen Eingangsimpedanz arbeitete 
das Mikrophon noch bis zu Driicken von der Gré- 
enordnung einer Atmosphare linear. Seine Kapazi- 
tat betrug 13 pF, sein Durchmesser 6,0 mm. 

In Bild 5 ist das Leerlauf-Ubertragungsma8 des 
Ultraschallmikrophons fir Beténung im freien 
Schallfeld senkrecht zur Membran wiedergegeben. 
Fiir das Gebiet des Hérschalls wurde das Ubertra- 


6 810 2 4 


6 8100 kHz 


Freifeld -Leerlauflibertragungsma! 


Bild 5. Freifeld-Leerlaufiibertragungsma des  Ultra- 


schallmikrophons M-213. 


gungsmaS durch Vergleichsmessungen mit einem 
kalibrierten Standard-Kondensatormikrophon We- 
stern Electric 640-AA erhalten; der Frequenzverlauf 
ist in Bild 5 als strichpunktierte Linie eingezeichnet. 
In diesem Frequenzgebiet ist das Freifeld-Ubertra- 
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gungsmaB des Ultraschallmikrophons um _ rund 
59 dB kleiner als das von gewéhnlichen Konden- 
satormikrophonen und verlauft im _ wesentlichen 
frequenzunabhangig. 

Durch die ausgezogene Kurve ist das Ergebnis 
der im Ultraschallbereich durchgefihrten Reziprozi- 
tatskalibrierung dargestellt. Das Ubertragungsmah 
verlauft bei hoheren Frequenzen unterhalb der Reso- 
nanzfrequenz des Mikrophons infolge der Schallbeu- 
gung an der zylinderfoérmigen, senkrecht zur Mem- 
branflache betonten Mikrophonkapsel nicht mehr 
vollig frequenzunabhangig. Die gestrichelte Linie in 
Bild 3 zeigt, wie sich nach Subtraktion der nach 
[13] dem Beugungseffekt entsprechenden Werte 
auch im Ultraschallbereich unterhalb der Resonanz- 
frequenz nahezu ein frequenzunabhangiges Verhal- 
ten des Mikrophons ergibt; diese gestrichelten Kur- 
venziige entsprechen etwa den Ubertragungsmafen 
fir tangentialen Schalleinfall, bei dem der Beu- 
gungseffekt verschwindet. 

Bei der Ausmessung eines bestimmten Schallfel- 
des mit Hilfe von drei unabhangig voneinander 
kalibrierten Mikrophonen ergab sich aus einem Ver- 
gleich der Mefresultate, dafi die Ergebnisse der 
Reziprozitatskalibrierung im Frequenzbereich unter 
100 kHz eine mittlere Genauigkeit von + 0,5 dB be- 
sitzen. 


3.2. Die zunehmende Aufsteilung der Wellenfront 
in Luft mit wachsendem Abstand von der Schall- 
quelle und mit wachsender Intensitiat 


Die durchgefiihrten Analysen der Aufsteilung der 
Wellenfront gliedern sich in zwei Gruppen: die 
Messung der Aufsteilung mit wachsendem Ab- 
stand von der Schallquelle bei unverandert bleiben- 
der Intensitat und die Messung der Aufsteilung mit 
steigender Schallintensitaét bei unverandert bleiben- 
dem Abstand von der Schallquelle. Sie wurden im 
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Bild 6. Oberwellenanteile einer urspriinglich sinusfér- 
migen Schallwelle fiir drei verschiedene Aus- 
gangsschalldruckpegel bei wachsendem Abstand 
vom Resonanzschwinger. 
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Schallfeld des Resonanzschwingers mit Hilfe des 
Ultraschallmikrophons M-213 durchgefihrt, wobei 
der Oberwellengehalt der aufgesteilten Welle mit 
Hilfe eines kontinuierlich abstimmbaren selektiven 
Rohrenvoltmeters bestimmt wurde. 

In Bild 6 sind fir drei Ausgangsschalldruckpegel 
die Ergebnisse der Schalldruckanalysen in Abhan- 
gigkeit vom Abstand vom Resonanzschwinger wie- 
dergegeben. Die Mefabstande betrugen 20 bis 
100 cm; sie reichten vom Ubergangsgebiet Nahfeld- 
Fernfeld bis ins Fernfeld des Schwingers. Bild 6 
zeigt zweierlei: Bei groBeren Ausgangsschalldruck- 
pegeln (rechte Kurvenschar) ist der Oberwellen- 
gehalt groBer als bei kleineren Ausgangsschalldruck- 
pegeln, die Kurven liegen dichter beieinander; 
auBerdem ist ersichtlich, daB bei konstant gehalte- 
nem Ausgangsschalldruckpegel der Oberwellen- 
gehalt der Schallwelle mit wachsendem Abstand x 
vom Resonanzschwinger zunachst zunimmt und dann 
konstant bleibt. Da der Schalldruck im Fernfeld 
einer Schallquelle mit 1/z abfallt und der Schall- 
druck, bei dem die erste Unstetigkeit der Schall- 
quelle eintritt, ebenfalls 1/2 proportional ist (Gl. 4), 
nimmt der Oberwellengehalt einer Schallwelle im 
Fernfeld einer Schallquelle nicht mehr zu. Auch bei 
starkem Oberwellenanteil ist der Effektivwert der 
ersten Oberwelle etwa gleich dem der Gesamtwelle; 
infolge der quadratischen Addition der einzelnen 
Effektivwertanteile einer zusammengesetzten Schall- 
welle beeinflussen die hoheren Harmonischen wegen 
ihrer kleineren Amplituden den Effektivwert der 
Gesamtwelle nur geringfiigig. 

Bild 7 zeigt das Anwachsen des Oberwellengehal- 
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Bild 7. Oberwellenanteile einer urspriinglich sinusfér- 
migen Schallwelle in 40 cm Abstand vom Reso- 
nanzschwinger bei wachsendem  Schalldruck- 
pegel. 


tes einer Schallwelle mit wachsender Intensitat in 
festem Abstand von der Schallquelle. Bei geringeren 
Intensitaten steigt die Amplitude der n. Oberwelle 
mit der n. Potenz des Schalldrucks; schlieBlich miin- 
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den die Oberwellenamplituden in Grenzwerte. Eine 
Behandlung dieser Erscheinung in der amerikani- 
schen Literatur [14] bestatigt dieses experimentell 
erhaltene Ergebnis. AuSerdem ist in Bild 7 zu er- 
kennen, daB im Abstand 40cm vom Resonanz- 
schwinger bei dem Schalldruckpegel 142 dB wber 
2°10~4 wbar die Amplituden der Oberwellen p, bis 
p, zueinander nahezu im Verhaltnis der Fourier- 
koeffizienten einer Sagezahnschwingung stehen: 


Pi: Pe: Ps? Pa=l: te: %/3: 4/4. 


Der Grenzwert der Aufsteilung der Wellenfront ist 
hier also fast erreicht. 
Zur Veranschaulichung dieser Erscheinungen zeigt 


Bild 8 die Schalldruckverlaufe in verschiedenen Sta- 


Bild 8. Schalldruckverlaufe einer urspriinglich sinus- 
formigen Schallwelle hoher Intensitét in ver- 
schiedenen Stadien der Aufsteilung. 


dien der Aufsteilung. Als Empfanger diente bei die- 
sen Aufnahmen ein auf einen Bleiklotz geklebter 
Bariumtitanat-Dickenschwinger, der auf 2 MHz ab- 
gestimmt war; Unebenheiten des Verlaufes seines 
Empfangstibertragungsmafes infolge Schallbeugung 
blieben hierbei unberiicksichtigt. 


4. Der Schallstrahlungsdruck in Luft 


Der Strahlungsdruck, der auf feste Korper in 
einem Schallfeld einwirkt, ist theoretisch schon aus- 
giebig untersucht worden. Nach den ersten Arbeiten 
von Rayieicn [15], [16] haben insbesondere die 
Untersuchungen von Britiours [17], [18] und 
Lancevin [19] wesentlich zur Klarung dieser Er- 
scheinung beigetragen. Von Herrz und Mennve [20], 
[21] wird darauf hingewiesen, da zwei Definitio- 
nen des Schallstrahlungsdruckes zu unterscheiden 
sind: Wird eine Schallwelle von unendlichem Quer- 
schnitt angenommen, so bezeichnen Hertz und 
MenpeE die Erhdhung des Druckes auf einen festen 
Korper gegeniiber dem am Korper bei gleicher mitt- 
lerer Dichte des Mediums im Ruhezustand herrschen- 
den Druck als den ,,Raytercuschen Strahlungs- 
druck“. Er ist von der Zustandsgleichung des Me- 
diums abhaéngig und verschwindet im Giltigkeits- 
bereich des Hooxeschen Gesetzes. Hat der Schall- 
strahl einen endlichen Querschnitt und grenzt er an 
in Ruhe befindliches Medium, so wird nach Hertz 
und Menpe an einem festen absorbierenden Korper 
im Schallfeld der ,,Lancevinsche Strahlungsdruck* 
beobachtet, der unabhangig von der betreffenden 
Zustandsgleichung stets gleich der mittleren Energie- 
dichte ist. Borents [3] macht darauf aufmerksam, 
dai der Strahlungsdruck einer Schallwelle auf einen 
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vollig absorbierenden Korper allgemein durch die 
doppelte mittlere kinetische Energiedichte gegeben 
wird, die nur fiir kleine Amplituden gleich der mitt- 
leren totalen Energiedichte ist. 

Im Gegensatz zu der Vielzahl von Arbeiten iiber 
die Theorie liegen nur wenige experimentelle Unter- 
suchungen des Strahlungsdruckes in Luft und seines 


Zusammenhanges mit der Energiedichte vor: Nach ~ 


ersten Messungen von Zernoy [22] aus dem Jahre 
1906 bestimmten Harrmann und Mortensen [23] 
im Jahre 1947 den Strahlungsdruck eines Schall- 
strahls in Luft mit Hilfe eines Radiometers und die 
zugehorige Energiedichte im Umweg iiber die Schall- 
schnelle mit einer Raytetcuschen Scheibe. Die Mes- 
sungen ergaben, daB im Experiment mit einem 
Schallstrahl von endlichem Querschnitt im ungestor- 
ten Medium der Lancevinsche Strahlungsdruck ge- 
messen wird. Gavreau [24] bestatigte dieses Resul- 
tat; er bestimmte die Energiedichte seines Schall- 
feldes durch Messung der Schwingungsamplitude. 


In jiingster Zeit untersuchten Bupat, Héy und Otson 
im Frequenzbereich 800 bis 8000 Hz den Strahlungs- 
druck auf Kugeln und Kreisscheiben, deren Abmes- 
sungen in der Gré8enordnung der Schallwellenlange 
liegen [25]. Die Resultate dieser experimentellen 
Untersuchungen stimmen in befriedigender Weise mit 
den Ergebnissen der Theorie des Strahlungsdruckes 
auf solche zur Wellenlange kleine Koérper iiberein. 


In der vorliegenden Arbeit wurde die mittlere 
Energiedichte von Ultraschallfeldern hoher Intensi- 
tat in Luft mit Hilfe einer Strahlungsdruckwaage 
und mit einem nach der Reziprozitatsmethode kali- 
brierten Ultraschallmikrophon untersucht und beide 
Ergebnisse miteinander verglichen. 


4. 1. MeBanordnung 


Bei den Strahlungsdruckmessungen wurden als 
Reflektoren Kreisscheiben und Hohlkegel verschiede- 
ner Abmessungen verwendet. (Hohlkegel dienen 
auch in einer Mefeinrichtung fiir den Strahlungs- 
druck in Fliissigkeiten [26] als Reflektoren.) Auf 
Grund der an den MefSstellen vorliegenden Schall- 
felder und nach Vorversuchen mit Reflektoren von 
2cm bis 12cm Durchmesser betrugen die Durch- 
messer der Duraluminreflektoren 3cm, 4cm und 
Sem. Nach Berechnungen von verschiedenen Auto- 
ren (siehe zum Beispiel [25], Fig. 2) weicht selbst 
fiir den kleinsten Reflektor die vom Strahlungsdruck 
darauf ausgeiibte Kraft nur um 7% von dem Wert 
ab, den sie im Idealfall bei einem im Vergleich zur 
Schallwellenlange groBen Reflektor annahme. 

Bei Schalldruckpegeln von 130 bis 150 dB itiber 
2:10‘ ubar sind Krafte auf die Reflektoren von 
etwa 5 mp bis 400 mp zu erwarten. Sie wurden mit 
Hilfe einer WrestpHatschen Waage, wie sie zur Be- 
stimmung der spezifischen Gewichte von Fliissigkei- 
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_ Bild 9. Waage zur Messung des Schallstrahlungsdruckes. 


ten dient, gemessen (Bild 9). An dem Ende des 
auferen Waagebalkens waren an Stelle des 
Schwimmkorpers die Reflektoren mit diinnen Drah- 
ten befestigt; die Waage wurde zunachst mit Hilfe 
von Gegengewichten ins Gleichgewicht gebracht. Da 
der Resonanzschwinger senkrecht nach oben strahlte, 
konnte der Reflektor unmittelbar ins Schallfeld ge- 
bracht und die durch den Strahlungsdruck ausgeiibte 
Kraft durch Aufsetzen von Reitergewichten kompen- 
_siert werden. Die Beobachtungen geschahen stets bei 
schwingendem Waagebalken, wie es bei feinen 
Wagungen tblich ist. 

Eine von der Schallquelle ausgehende Gleichstr6- 
mung (,,Quarzwind“), die in viskosen Medien auf- 
tritt [27], war_zu klein, um das MeBergebnis mefi- 
bar zu beeinflussen: Die Absorption des durchstrahl- 
ten Mediums, die nach Fox und Herzretp [28] die 
Ausbildung des Quarzwindes verursacht, ist bei der 
verhaltnismaBig niedrigen Frequenz von 30 kHz und 
den kurzen durchstrahlten Strecken zu gering. 


Die Bestimmung der Energiedichte im Schallstrahl 
begann jeweils mit der Ermittlung des Schalldrucks mit 
Hilfe des kalibrierten Ultraschallmikrophons M-213. 
AnschlieSend wurde nach Entfernen des Mikrophons 
und Einbringen der Waage in den Schallstrahl der 
Strahlungsdruck bestimmt und danach die Schalldruck- 
messung wiederholt; wahrend der Messung wurde die 
Spannung an der Erregerwicklung des _ Resonanz- 
schwingers konstant gehalten. Die kleinen Anderungen 
des yom Schwinger erzeugten Schalldrucks wihrend der 
Me8Bzeit konnten auf diese Weise bestimmt und im 
Ergebnis beriicksichtigt werden. 


4.2. Ergebnisse 


Die Messungen wurden an verschiedenen Orten 
im Nah- und Fernfeld des Resonanzschwingers 
durchgefiihrt; als Beispiele werden die Ergebnisse 
fiir die MeBabstande 9,5cm und 38cm mitgeteilt. 
An der MeBstelle in 9,5 cm Abstand vom Schwinger, 
auf einer Flanke der Schalldruckverteilung lings 
der Achse des Resonanzschwingers, blieb auch bei 
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den hochsten erreichten Intensitaten der Schalldruck- 
verlauf sinusformig. Bei dem MeR<abstand 38 cm 
hat der Oberwellengehalt der Schallwelle seinen 
Hochstwert gerade erreicht, und der Effektivwert des 
Schalldrucks ist noch nicht stark abgesunken (siehe 
Bild 6); wie Bild 7 zeigt, hat in 40cm Abstand die 
groBte erzeugte Intensitat fast zu dem Grenzwert der 
Aufsteilung der Wellenfront gefiihrt. Die Schall- 
druckverteilung auf die Reflektorflachen war an den 
MeBstellen innerhalb + 0,5 dB gleichmaBig. 

Die Ergebnisse der Messungen sind in den Bil- 
dern 10 und 11 wiedergegeben; hierin ist der Ab- 
solutwert der mittleren Energiedichte des Schallfel- 
des in der Abszisse als Resultat der Schalldruckmes- 
sungen, in der Ordinate als Ergebnis der Strahlungs- 
druckmessungen dargestellt; der Abszissenwert 
wurde mit L,,, der Ordinatenwert mit E, bezeichnet. 
Es zeigt sich, daf} sich mit den Reflektoren unabhan- 
gig von ihrer ebenen oder Hohlkegel-Ausfiithrung 
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Bild 10. Mittlere Energiedichte E des Schallfeldes in 
9,5cm Abstand vom Resonanzschwinger als 
Ergebnis von Schalldruck- (Index p) und 


Strahlungsdruckmessungen (Index r). 
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Bild 11. Wie Bild 10. 38cm Abstand vom Resonanz- 
schwinger. 
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jeweils nahezu die gleichen Werte ergeben. Die Re- 
sultate der Messungen mit den Reflektoren vom 
Durchmesser 3cm weichen infolge des erwahnten 
Beugungseffektes von den tbrigen etwas ab. Bei 
kleinen Intensitaéten nimmt der Streubereich der 
Mefiwerte zu, da sich die zugehorigen geringen 
Krafte mit der Waage nur mit verminderter Genau- 
igkeit messen lassen. Aus der Streuung der Mittel- 
werte, die den Bildern 10 und 11 zugrunde liegen, 
folgt fiir die Ergebnisse eine Unsicherheit von 
+1 dB. 

Unabhangig davon, da im Abstand 9,5 cm 
(Bild 10) der Schalldruckverlauf stets sinusformig 
blieb, wahrend die Aufsteilung der Wellenfront im 
Abstand 38cm (Bild 11) bei den hochsten Inten- 
sitaten ihren Grenzwert fast erreicht hatte, fiihren 
Schalldruck- und Strahlungsdruckmessungen inner- 
halb der MefSgenauigkeit jeweils zu dem gleichen 
. Ergebnis. Daraus folgt, dafS§ auch bei starkem Ober- 
wellengehalt einer Schallwelle infolge der Aufstei- 
lung der Wellenfront die GroBe des Strahlungsdruk- 
kes nicht meBbar von dem Wert abweicht, der der 
mittleren Energiedichte einer reinen Sinuswelle glei- 
chen Effektivwertes entspricht. 

Der Verfasser dankt Herrn Prof. Dr. M. Griz- 
MACHER fiir die standige Forderung dieser Unter- 
suchungen; auSerdem bedankt er sich bei Herrn Dr. 
H. Mrass fiir die Anregung dieser Arbeit und fir 
viele Diskussionen. Ebenso verdankt er Herrn Dr. 


H. G. Dissrex zahlreiche wertvolle Hinweise. 
(Eingegangen am 9. November 1960.) 
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staEpT, W.Rercuarpt, W.A. Rosensuitu, J. F. 
Scuouren, M. Scuiitzensercer, K. Srermpusca und 
N. Wiener. Springer-Verlag Berlin-Gottingen-Hei- 
delberg. Erscheint seit Januar 1961. 

Der von N. Wiener in die Wissenschaft eingefihrte 
Begriff der Kybernetik hat in wenigen Jahren zu einer 
stiirmischen Entwicklung gefiihrt, deren Ergebnisse 
sowohl fiir die Nachrichteniibertragung und Informations- 
verarbeitung als auch fiir die Steuerung und Regelung 
insbesondere bei der Automation von Bedeutung sind. 
Aber auch die entsprechenden Vorgange im lebenden 
Organismus gehéren zum Aufgabengebiet der Kyber- 
netik und gerade der Vergleich zwischen den im Orga- 
nismus yerwirklichten Prinzipien der Steuerung und 


Informationsverarbeitung mit den in der Technik ent- 
wickelten macht das Gebiet so tiberaus reizvoll und an- 
regend. Es ist daher zu begriiBen, da eine neue Zeit- 
schrift erscheint, die sich die Pflege und den Erfahrungs- 
austausch der einzelnen Teilgebiete angelegen sein laBt. 
Die groBe Zahl] hervorragender Wissenschaftler des In- 
und Auslandes, die bei der Herausgabe mitwirken, laBt 
ein hohes Niveau der Beitrage erwarten. Diese Ver- 
mutung wird durch das erschienene erste Heft voll be- 
statigt, das unter anderem Arbeiten von K. Kiprmituer, 
W. Reicuarpr, G. Vosstus, K. Sreinsucu, D. H. Carson 
und I. Konter enthilt. Die Zeitschrift, deren Beitrige in 
deutscher oder englischer Sprache erscheinen, schlieBbt 
eine fihlbare Liicke im Schrifttum und diirfte auch 
dem Akustiker, vor allem dem physiologisch interessier- 
ten, wertvolle Erkenntnisse und Anregungen yermitteln. 

W. Ka.vensacu. 


Pater ATION OF THR AUTOCORRELATION METHOD 
Pei seh eC URAL ANALYSIS OF SOUND RECORDS 
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King’s College, Newcastle upon Tyne 


Summary 


The application of the autocorrelation method of acoustical spectral analysis is discussed 
and examples are given of typical stationary functions suitable for being analysed by this 
method. 

The stationary function is represented in the form of the time samples and its autocorre- 
lation coefficients for different delay times are computed. 

The power density spectrum is next computed as a Fourier transform of the autocorre- 
lation function represented as a set of separate autocorrelation coefficients. An electronic 
computer is used in both cases. 

The results are plotted in graphs and compared with the original time functions. 


Sommaire 


On expose Ja question d’appliquer la méthode d’autocorrélation 4 l’analyse des spectres 
acoustiques et on donne des exemples de fonctions du type stationnaire susceptibles d’étre 
analysées par cette méthode. La fonction stationnaire est représentée sous la forme d’une 
suite d’instants et ses coefficients d’autocorrélation ont été calculés pour divers temps de 
retard. : 

La densité d’énergie du spectre est alors calculée par transformation de Fourier de la 
fonction d’autocorrélation constitutée par l’ensemble des différents coefficients d’autocorré- 
lation. Un calculateur électronique sert 4 faire ces deux calculs. 

Les résultats sont portés sur un graphique et comparés avec la fonction du temps 
initialement proposée. 


Zusammenfassung 


Die Anwendung der Autokorrelationsmethode auf die akustische Spektralanalyse wird 
diskutiert, und es werden Beispiele gegeben von stationaren Funktionen, die fiir eine solche 
Analyse geeignet sind. 

Die stationdren Funktionen werden dargestellt in Form von Zeitfunktionen und ihre , 
Autokorrelationskoeffizienten werden fiir verschiedene Verzégerungszeiten berechnet. 

Das Energiedichtespektrum wird dann als Fouriertransformierte der Autokorrelations- 
funktion berechnet, die dargestellt wird als Satz von einzelnen Autokorrelationskoeffizienten. 
Ein elektronischer Analogrechner wurde in beiden Fallen benutzt. 

Die Ergebnisse werden graphisch dangestellt und mit den urspriinglichen Zeitfunktionen 


verglichen. 
1. Introduction 


The application of the autocorrelation method 
in physical analysis has become more popular in 
recent years, due to the general development of the 
theory of stochastic processes and growing interest 
of physicists in interpreting random phenomena, 
such as noise. The autocorrelation method of analy- 
sis may be applied to the stationary random pro- 
cess, whose distribution function is invariant under 
translation of the t-axis. In physical terms, con- 
sidering for instance noise, this means that the fac- 
tors which affect the time average, such as the dis- 
tance from the source of noise, the damping con- 
ditions ete., should remain constant during the time 
of analysis. 

The application of the autocorrelation function 


* Now at Academy of Music, Warsaw. 


in acoustical analysis of speech was first widely dis- 
cussed at the Speech Communication Conference at 
Massachusetts Institute of Technology in May — 
June 1950, and in some papers puplished in the 
Journal of the Acoustical Society of America [1], 
[2], [3], [4], [5]. 

In later years advanced apparatus has been built 
using the correlation technique for analysing acous- 
tical phenomena [6], [7]. The basic mathematical 
theory of autocorrelation analysis has been mainly 
developped by the American writers Wiener [8], 
[9], Rice [10], [11], Lawson and Untensecx [12]. 


2. Autocorrelation function and power spectrum 


There are two basic ways of approach to the 
problem of acoustical spectral analysis. The first, an 
essentially experimental one, concerns the use of 
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electronic — sometimes mechanical — filters tuned to 
different frequencies, or one filter scanning the 
whole frequency region on the basis of the hetero- 
dyne. No doubt, the sound spectrograph is one of 
the most powerful tools built on this principle. 


The second way, perhaps more traditional, con- 
cerns analysing physical processes as FourtEr trans- 
forms of a time function. Such an analysis is of 
course very simple when dealing with a periodic 
time function and its result is a Fourier series. It 
is also quite easily carried out for many nonperio- 
dic functions with the aid of Fourier integrals. 
There is however a much broader group of statio- 
nary functions for which it is not possible to find 
the direct Fourier transform, this group contains 
random noise and speech. The reason is that Fov- 
RiER’s theory, as it stands, is not applicable to 
functions which are specified in terms of statistics 
and probability and are not representable by spe- 
cific analytic expressions giving their precise values 
for all values of the independent variable. Spectra 
of such functions, is found experimentally, cannot 
be transformed back into the original time func- 
tions. It was proved by Wiener [8] in 1930, that 
the result of transforming such a spectrum is the 
autocorrelation function of the original random 
function. This very important theorem gives the 
extension of the Fourter methods to the harmonic 
analysis of random processes. 


The autocorrelation function of a time function 
may be defined as 


Pi 
eal 
D(z) aor [1 fae) de eM) 


It is an even function of t and represents the aver- 
age value of the product of /(¢) and f(t+ 7). 

It may be useful to investigate the form of auto- 
correlation function for some idealised cases of 
typical random processes. Lez and Wiesner [13] 
give an interesting example of the flat— topped 
wave of alternating positive and negative pulses of 
varying duration (see Fig. 1). It is assumed that 
the distribution of the zero crossings follows the 
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Fig. 1. Function having random zero crossings. 


Poisson law, which means that if the average num- 
ber of zero crossings per second is k, the pro- 
bability P(n,t) that there are n zero crossings in 
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any duration of t of the wave is given by the for- 
mula 


(kr)” 


n! 


P(n,t) = ou (2) 
and the autocorrelation function from eq. 1 may 
be proved to have the form 


P(r) = Bee kel , (3) © 


This autocorrelation function is shown on Fig. 2, 
and its shape is-typical for the autocorrelation of 
all random processes. It is damped and non- 
periodic. 


G(T) 


i Se 


Fig. 2. Autocorrelation function of a random process 
of Fig. 1. 


The definition in eq. 1 may be applied not only 
to stationary random functions, but also to perio- 
dic ones. In the case of a periodic function we need 
only consider one complete cycle, and the auto- 
correlation function is 


T, 
B= 2 [hOhernd. — @) 
1 
0 


Let us put 
f1(t) =A, cos(@, t+), (5) 
and we find from the eq. 4 


@i(x) 3a 
(a= ss COS M,T. (6) 


This relation is shown in Fig. 3, we notice that 
the autocorrelation function of a simple periodic 
wave is an undamped cosine-wave with phase angle 
equal to zero. Similarily assuming f, as an arbi- 
trary periodic function in a general form and 
applying eq. 4, we learn that its autocorrelation 
function is a cosine series preserving all frequen- 
cies and dropping all phase angles of the harmonics 
of the original function. 


A particularly important problem arises when 
we have to deal with a random process including 


t 
f,(t) (Tt) 
t— tT. 


Fig. 3. Sine wave and its autocorrelation graph. 
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hidden periodic components. It may be proved that 
the autocorrelation function of the sum of the noise 
and periodic function is the sum of their individual 
autocorrelation functions. Let us consider the be- 
haviour of such an autocorrelation at compara- 
tively large values of the variable t. The noise auto- 


correlation curve, similarly to this shown on Fig. 2 


approaches theoretically a constant value as t be- 
comes larger. The periodic wave, on the contrary, 
gives the autocorrelation curve in the form of an un- 
damped cosine series, which can easily be detected 
after the larger damped part of noise autocorrela- 
tion has died out (see Fig. 4). 


Te 


Fig. 4. Autocorrelation function of random noise con- 
taining a periodic wave. 


If the periodic part of the autocorrelation func- 
tion is also more or less damped with increasing t , 
this indicates that instead of single periodic com- 
ponents some more or less broad bands of the con- 
tinuous spectrum are carrying a considerable part 
of the energy. We can easily notice in these ex- 
amples the immediate and close relation between the 
autocorrelation function and the spectrum. 

In the case when the time function is given by a 
set of discrete data, samples of a continuous func- 
tion f(t), the autocorrelation function from eq. 1 
may be represented in another form: 

1 ¥ 
DG) oe >» Am Om » (7) 
m=1 
where a», and 6, are amplitudes of f(t) separated 
by the time t. 

Suppose that the time function is represented by a 
set of discrete data 2;, :=1—N. Then the time 
difference t between multiplied values of the func- 
tion is represented as an integer k, and the auto- 
correlation function takes the form 


N=k 
1 
ras 3 ViVi shee 
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@(k) = 
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(8) 

Eq. 8 is true only if 2; and x;,, are samples taken 
on the assumption that their mean value m=O; in 
a general case they should be replaced by 2x;— 2; 
and 2;,%—%j.,. Usually it is more convenient to 
introduce the normalised autocorrelation function, 
dividing ®(k) by the variance o? = E[ (x —m)?]: 
P(k) 
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The normalised autocorrelation function so defined 
has |e,|<1land@=1. 

As stated above, the most important property of 
the autocorrelation function ®(t) is its reciprocal 
relation to the power spectrum G(w), known as 
Wiener’s Theorem. According to this, the auto- 
correlation function and the power spectrum form 
a pair of Fourier transforms: 


@(t) = = [ e() Toneda ce: (10) 


co 


G(@) = fr cosmtdt. 
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(11) 


Eq. 11 indicates that knowing the autocorrela- 
tion function of a time function is in every sense 
equivalent to knowing the power density spectrum. 
We may compute the spectrum from the autocorre- 
lation function, or derive from @(t) some direct 
information which is equivalent to that contained in 
the spectrum. Both ways have been used in the 
investigations described in this paper. 

The concept of a power spectral density was de- 
fined by Rice [10], [11] in electrical terms as the 
power within a frequency band dj that would be 
dissipated in a 1-Ohm load. In the case of mechani- 
cal vibrations the precise interpretation is often 
much more difficult, however sometimes it is enough 
to know the relative spectral distribution of the 
energy, without estimating any absolute values. 
This was the case in the present analysis. For this 
reason, too, the logarithmic scale of frequency for 
plotting the power density might have been used. 


3. The subject of analysis 


The purpose of the investigation was to analyse 
film records representing the amplitude modulation 
of an ultrasonic signal by a turbulent flow. The 
analysed wave had the character of a random func- 
tion, in some cases indicating a certain amount 
of periodicity. The film records were obtained by 
Ricnarpson [14] during his investigations of tur- 
bulence in the wake of a cylinder. The recording 
was carried out in the following way. Two ultra- 
sonic 2Mc/s-quartz transducers were placed on 
opposite walls of a 76cm wide wind tunnel. The 
signal picked up by a receiver was fed to a cathode 
ray oscillograph and photographed, The conditions 
of the flow were next changed by introducing into 
the tunnel a 5.1 cm-diameter wooden cylinder at 
different distances up stream from the ultrasonic 
beam. The row of vortices generated in the stream 
of wind passing an obstacle caused the fluctuations 
of the received signal to be in some degree regular. 
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Records were made with different speeds of wind 
and at different distances from the cylinder. 

Another group of records was obtained by plac- 
ing the two transducers in the boundary layer at 
three different levels not extending more than 5 cm 
above the floor of the tunnel [15]. 


4. Description of the method 


The film records were projected on a graduated 
screen and samples x; were taken every 2 mm along 
the t-axis (see Fig.5). According to the film re- 
cording speed 123 mm/s and approx. enlargement 
the time 


6.4, this corresponds to difference 


Nie 25-108? s 


Fig. 5. Taking the samples of the time function. 


In several cases the enlargement was altered to 
meet the obvious changes in the frequency region of 
modulation. The new value indicated the rate of 
change and an appropriate coefficient had to be 
introduced in computing the standard deviation of 
the data, to reduce the values of the computed spec- 
tral density to the same scale. 

The number of observations NV depended on the 
length of the time recorded in different examples 
and varied from NV = 300 to N = 885. 

The computations of the autocorrelation function 
and of spectral density were carried out with the 
aid of the electronic computer “Pegasus” in the 
Durham University Computing Laboratory, using 
the programmes developed in the Royal Aircraft 
Establishment, Farnborough, Hants. 

The set of data punched on the tape were first 
changed from decimal into binary. Then values of 
the normalised autocorrelation function for diffe- 
rent values of delay time k At were computed using 
the relation 


N-k 


(N—k) > UL +k — > 
T T 
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The computed values of the autocorrelation func- 
tion @;, which are usually called “autocorrelation 
coefficients”, were then plotted on graphs. Later 
they were used as data for computing power spectra. 


The standard deviations o of the sets of data were 
also computed, but in the cases when different en- | 
largements had been used, o had to be multiplied by — 
a coefficient eens in order to make all the - 

Zhe 

results comparable. 
Calculations of the spectral density were carried _ 
out using the relations 


n At o? 
Cis Gy 5 
aU 
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k=1 


r=2nAtf, (15) 


which may be derived from eq. 11 when the auto- 
correlation function is given as a set of separate 
autocorrelation coefficients 0,, k= 1—n. a, repre- 
sents the normalised spectral density which was cal- 
culated by the computer given a set of autocorrela- 
tion coefficients. The real spectral density G, at the 


frequency {= then obtained from 


eq. 13. 


A, is a smoothing factor introduced to compen- 
sate the effect of finite bounderies of integration. 
This problem ist widely discussed by Barrierr 
[16], and J;, in the form 


was 


T 
2nAt 


k 


A,p=1— — (16) 
n 


introduced in our case, is called Banrriert’s 
smoothed periodogram. 

The plots of spectral density are all comparable, 
being drawn in the same scale. Zero decibels corres- 
ponds to the largest maximum density detected in 
all examined examples. 


N-k N-k 
vj By Vi+k 
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which may be derived from eq. 9. 

This computation was carried out in most cases 
for k=1— 40, which usually corresponded to the 
maximum delay time 40 At=40°2.5°3=0.1s, in 
some examples however 0, was computed for 


k=1— 80. 


1 


5. Discussion of the results 


The total number of records analysed with the 
method described was 60, and the discussed ex- 
amples are quoted on the Table I. The most inter- 
esting results are perhaps those obtained from the 
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Table I. List of examples. 
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amber Position 
Speed of Wind of the ultrasonic N At n 
of example b 
eam | 

1 low in the wake, 482 -2.5 ms 40 
near the cylinder 80 

2 medium in the wake, 770 2.5 ms 40 
near the cylinder 

3 | high in the centre 500 2.5 ms 40 
of the tunnel, | 1.3 ms 
without the cylinder | 

4 medium in the boundary 300 2.5 ms 40 

| layer 

5 medium | in the boundary | 300 2.5 ms 40 
layer, nearer the | 
floor of the tunnel 


wake records, as the low fre- 
quency of vortices may be often 
observed here even on the original 
time functions. The ultrasonic fre- 
quency in all the quoted examples 
was 1 Mc/s. 

Fig. 6 illustrates the effect of 
changing the speed of the wind. 
The examples 1 and 2 were re- 
corded in the small distance be- 
hind the cylinder at low and me- 
dium. speeds of the wind. The 
autocorrelation function, falling 
down rapidly, after the value of 
t =0 and with complex periodical 
shape indicates the large band 
noise spectrum with complex peri- 
odicity. This periodicity comes 
out as a set of distinct maxima 
in the power spectrum which 
shows something similar to “over- 
tone structure” with the “fun- 
damental frequency” being that 
of a row of vortices in the wake. 

In the case of the record 1, 
the autocorrelation coefficient 0; 
have been computed for k = 180 
and the spectral density was cal- 
culated in two ways, first for the 
number of autocorrelation coeffi- 


20 50 100 
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cients n=40 (the dotted line), 
and then for n=80 (the plain 
line). In the latter case the long 
time of delay enables us to detect 
the maximum density at the fre- 
quency f=12c/s which roughly 


corresponds to the vortices frequency at this speed. 
It should be emphasised however, that the short 


duration of the example 


(number of samples 
N=482) does not authorise us to derive any firm 


k ——e 
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Fig.6. The time functions (A), autocorrelation graphs (B), and power 


density spectra (C) of the records showing the effect of changing 
the speed of wind in the wind tunnel. 1 and 2: turbulences in the 
wake; 3 turbulences in the tunnel without the cylinder. 


conclusions from the calculations involving such a 
long time of delay. 

In the record 2 the speed of wind was doubled. 
The maxima of the spectral density are transformed 
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floor of the tunnel. They too re- 
present the functions of a random 
character with some amount of 
periodicity. 

The analysing properties of the 
correlation functions were used 
in a number of electronic devices 
built during the last fifteen years 
mainly in the United States of 


America. They are rather compli- 
cated and expensive. The method 
of analysis described above, 


though rather slow and involving 
a number of tedious calculations, 
may be useful where even a small 


Power density —> 


x 


electronic computer is available 
and where the recorded data are 


in a form not suitable for analysis 
by the use of the electric filters, 
e. g. in short series on the film. 


Sue 10, = °20 


Frequency — 


Fig.7. The time function (A), autocorrelation graphs (B), and power 
density spectra (C) of two boundary layer records. 


towards higher frequencies and the total value of 
power density is affected. There is a very strong 
peak of the spectral density at the frequency reg- 
ion 40—50 c/s. Which corresponds to a “second 
harmonic” of the vortex frequency, but may be 
mainly due to some local mechanical resonance. 
This frequency may be also clearly observed both on 
the record itself and on its autocorrelation func- 
tion. 

The record 3 was made without a presence of the 
cylinder in the wind tunnel at high speed of the 
wind, It represents a typical random noise func- 
tion with some high frequency periodicities. 

The random character of record 3 is illustrated 
by the rapid falling of its autocorrelation function 
with t (no correlation between the instanteneous 
amplitudes at points that are spaced more than 
1.5 msec), its later fluctuations however show some 
evident periodicities at high frequencies. As the 
strongest periodicity revealed in the autocorrelation 
function 3 corresponded to the wavelength of the 
doubled distance between two consequent readings 
(2 At), it might have been subject to a compara- 
tively large sampling error. In order to check this 
another set of samples was taken from the record 
with different enlargement (At=1.3 ms) -@®(r) 
obtained in this way shows a great similarity to the 
previous one, and the computed density spectrum 
is shown (the dotted line). 

The records 4 and 5 (see Fig. 7) were made in 
the boundary layer at different distances from the 


The accuracy of the results 
grows with increasing number of 
samples, but it sometimes depends 
on the length of series available 
for study. 


50 100 200c/s 
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IMPROVED DIRECTIVITY USING SYNCHRONOUS 


DEMODULATION 
by D. G. Tucker 


Electrical Engineering Department, University of Birmingham 


Summary 
It is shown that it is possible to eliminate completely all secondary responses in the 


- directional pattern of a receiving line-array by the use of synchronous demodulation in 


suitable circumstances, notably those in which access to both ends of a pulse-transmission 
link is available and interference is due to spurious transmission paths. This is achieved 
by first making all secondary responses have a phase (or polarity) reversed relative to that 
of the main beam; the simplest and most efficient way to do this is to omit a certain propor- 
tion of the array from its centre, but this gives only an approximation to the required con- 
dition. It does not seem possible to obtain all secondary responses in reversed phase without 
the use of amplifiers or gross inefficiency. 


Sommaire 


On montre qu’il est possible d’éliminer complétement les réponses de second ordre dans 
les diagrammes de directivité de groupes de lignes réceptrices en utilisant une démodula- 
tion synchrone dans des conditions convenables; ceci est le cas en particulier quand les 
deux extrémités d’une chaine de transmission d’impulsions sont d’accés facile et quand des 
interférences sont dues 4 de mauvais trajets de transmission. On arrive au résultat en fai- 
sant en sorte que les réponses secondaires aient une phase (ou polarité) inversée par rap- 
port a celle du faisceau moyen; la maniére la plus simple et la plus efficace pour y arriver 
est d’éliminer une certaine partie du groupe autour de son centre mais cela donne seule- 
ment une solution approchée de la condition cherchée. I] ne semble pas possible de réaliser 
Vinversion de phase de toutes les réponses secondaires 4 moins d’employer des amplifica- 
teurs dont l’efficacité soit détestable. 


Zusammenfassung 


Es wird gezeigt, daf’ es unter gewissen Umstanden modglich ist, alle Nebenmaxima im 
Richtdiagramm eines Antennensystems durch phasenempfindliche Gleichrichtung zu unter- 
driicken. Dazu wird zunachst die Phase (oder Polaritét) aller Nebenmaxima relativ zum 
Hauptmaximum umgekehrt; dies kann zum Beispiel einfach durch Fortlassen einiger Ele- 
mente aus der Mitte der Richtantenne geschehen, allerdings werden die geforderten Bedin- 
gungen dadurch nur angenahert. Es scheint nicht méglich zu sein, alle Nebenmaxima mit 
um 180° gegen das Hauptmaximum gedrehter Phase zu erhalten, ohne dazu Verstarker zu 
yerwenden oder einen schlechten Wirkungsgrad in Kauf zu nehmen. 
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[16] Barrterr, M.S., An introduction to stochastic 


1. Introduction 


Much has been written lately on the topic of 
directional patterns of arrays and receiving systems, 
and clearly directional response and directional dis- 
crimination are, in general, separate matters with a 
relationship dependent on the nature of the receiving 
system as a whole. Only when we restrict considera- 
tion to the unprocessed outputs of an array of ele- 
ments added together linearly (i. e. not even rectified) 
does the directional response truly represent directio- 
nal discrimination also. (By directional response we 
mean the variation of output voltage produced by a 


single small signal source as the array is rotated; by 
directional discrimination we mean the ability of the 
array to discriminate against unwanted sources — or, 
more usually, merely against isotropic noise —in 
favour of a source in a particular direction.) The 
“processing”’. of the output signal to enable it to be 
used is inevitably non-linear and leads to the sepa- 
ration of response and discrimination to a greater or 
lesser extent. The matter was discussed in some 
detail by the author and a colleague in a paper on 
multiplicative receiving arrays [1]; these are arrays 
divided into two parts, the output signals from which 
are multiplied together, instead of being added. 
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Recent publications [2], [3] show many attempts 
to obtain © directional responses with particular 
merits (e. g. narrow main lobe, and/or reduced am- 
plitude of secondary lobes) by suitable processing, 
in spite of poorer directional discrimination, and 
one method discussed in the paper on multiplicative 
arrays is particularly interesting: it is to use a rec- 
tifier to exploit the fact that in the multiplicative 
system the wanted output signal is a positive d. c., 
and signals received from directions within the first 
(and largest) secondary lobe give a negative pola- 
rity of d.c. If the rectifier is connected suitably, 
therefore, the positive d.c. from the main lobe is 
passed, but the largest secondary response is elimi- 
nated. As the main lobe is of only half the width of 
that obtained from an ordinary linear (additive) 
array of the same total length, and the residual 
positive secondary responses are very small, a very 
_ good directional response is obtained, — and with 
not too serious a deterioration of directional dis- 
crimination. 

The present paper describes what is really an 
attempt to obtain a similar kind of effect from an 
ordinary additive array. Since the rectifier principle 
can be applied to reduce secondary responses only 
when the wanted output signal is d.c. of polarity 
varying from lobe to lobe, it is clear that synchro- 
nous demodulation (as opposed to ordinary detec- 
tion) of the signal is essential. It will become evident 
that while the system described is of considerable 
academic interest, it is of much less practical impor- 
tance than the multiplicative system owing to its 
more restricted field of application. It will, however, 
be most conveniently described in terms of a prac- 
tical system. 


2. Obtaining improved directivity. The method 


There are practical situations where synchronous 
demodulation [4] can be used because the carrier 
of the transmitted signal is available for use at the 
receiver, and where directional discrimination is re- 
quired at the receiving array to reduce or eliminate 
the effect of a spurious transmission path in the 


medium. An example is shown in Fig. 1, which illu- ~ 


strates an experiment on transmission phenomena 
in an underwater acoustic path. If the transmitted 
carrier can be passed separately to the receiver over 
a cable, then it can be used as the local reference 
tone for synchronous (or coherent) demodulation, 
as shown. If only the direct transmission path exists 
in the water, then the phase adjuster is set to give 
maximum d. c. output from the low-pass filter. Fluc- 
tuations in this output represent fluctuations in 
either amplitude or phase of the received signal; but 
as the amplitude is proportional to cos 0, where 0 
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Fig..1. Schematic diagram of synchronous-demodulation 
system used in conjunction with underwater 
acoustic transmission tests. 


is the phase angle between the received signal and 
the local tone, the output is relatively insensitive to 
phase fluctuations provided they do not exceed say 
+ 20°. Phase and amplitude fluctuations can, how- 
ever, be separated if desired by the use of another 
demodulation path with the local tone in quadrature 
with that in the first path; the d.c. output is then 
proportional to sin O, and using both outputs, am- 
plitude and phase can be calculated separately. 

If, with the single-path transmission, the receiving 
transducer is rotated in the plane of its length, a 
directional response is obtained. If this is measured 
before the demodulator, the result shown in Fig. 2 a 
will be obtained [5], assuming a linear strip array 
of length J, and an abscissa scale of x = (1 2/A) sin 9, 
where @ is the angle of the array from its normal 
position facing the transmitter. Alternate lobes of 
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Fig. 2. Directivity patterns of linear array. 
a Signal envelope as array is rotated, 
b after synchronous demodulator, 
c after synchronous demodulator and rectifier. 
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the response are of reversed carrier phase. After the 
synchronous demodulator, assuming the local tone 
remains unchanged, the output will be as in Fig. 2b, 
where negative d. c. is obtained from the lobes with 
reversed polarity. If the negative outputs are re- 
moved by suitably connecting a rectifier following 
the low-pass filter, then the response of Fig. 2c is 
obtained; this is obviously a much-improved re- 
sponse. 

In these circumstances the question naturally 
arises as to whether it is possible to cause all the 
secondary responses to be of reversed phase, so that 
in the synchronous demodulation system all secon- 
dary responses may be removed by the use of a rec- 
tifier. One simple means by which this may be done 
is shown in Fig. 4; a small (i. e. omnidirectional) 
transducer is fitted at the centre of the array, and 
its output, after amplification, subtracted from that 
of the main array. If the latter (P in the diagram) 
is taken as unity for x=O, then the former (Q in 
the diagram) requires to be of amplitude 0.128 to 
convert all secondary responses to reversed polarity, 
as shown in Fig. 3a. After synchronous demodula- 
tion, the d.c. response of Fig. 3b is obtained, and 


after a rectifier that of Fig. 3c, which has no secon-. 


dary responses at all. 
It must be realized that this improved response 
applies as directional discrimination only when the 
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Fig. 3. Directivity patterns of linear array with additio- 
nal point transducer as in Fig. 4. 
a Signal envelope as array is rotated, 
b after synchronous demodulator, 
c after synchronous demodulator and rectifier. 
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signal to be discriminated against is approximately 
cophasal with the main signal; in fact a tolerance 
of +90° exists since it is only necessary to obtain 
reversed-phase lobes. The new directivity can be 
exploited, therefore, only when the interfering signal 
comes from the same source as (or a source locked 
to) the main signal. In the example of Fig. 1, then, 
where it may be supposed that the indirect path is to 
be eliminated by bringing it outside the main lobe 
of the directional response of the receiver, it is 
necessary for the indirect path to maintain its phase 
constant to +90°; it can be brought initially into 
the correct phase relationship by a small adjust- 
ment of frequency. Obviously if the water surface is 
not steady, difficulties will arise, and the success of 
the system will depend on a large number of fac- 
tors which cannot be discussed here. One thing in 
particular must be pointed out; the system cannot 
discriminate in favour of the main path and against 
the indirect path with signals simultaneously pre- 
sent on both. This is because the rectifier removes 
the negative signals only when they are not can- 
celled out by positive signals. So it is essential for 
the signal to be pulsed, and for the difference in 
transmission time between any indirect path and the 
direct path to be at least one pulse length. The same 
result of eliminating unwanted paths can, of course, 
be achieved in the usual way by gating-out the want- 
ed pulse, but this means adjusting the gating circuit 
to operate at the correct instant and requires a gap 
between wanted and unwanted pulses of consider- 
able duration; the new process is automatic, and 
requires no gap between wanted and unwanted 
pulses. 

An alternative way of achieving the result de- 
scribed above would be to choose a frequency at 
which the indirect signals are approximately in anti- 
phase with the main signal, and to use a type of 
array which has all its responses in the same phase 
as the main beam; an example is the linear-taper 
array [5]. But although this solution may have some 
attractions, yet, for a given length of array, it gives 
a much wider main beam. 


3. Can a single array have all secondary 
responses of reversed polarity? 


The process used above to get all secondary re- 
sponses reversed in polarity relative to the main 
lobe is inelegant and gives a poor signal/noise per- 
formance, since the noise output of the omnidirec- 
tional transducer is amplified before combining with 
the main array. The question naturally arises as to 
whether it is possible to get the same effect by modi- 
fying the distribution of sensitivity with position 
along the length in the main array. As far as the 
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author can determine, this is not possible with a 
reasonably simple and efficient arrangement; it 
appears to him that no single array can have all its 
secondary responses (i. e. from |x|=z to |x|= ~) 
of reversed polarity unless the extreme action is 
taken of simulating the arrangement of Fig. 4 by 
reversing the phase of a very small central portion 
of the array and reducing the sensitivity of the re- 
mainder. This would be grossly inefficient through 
wasting available sensitivity. But it is possible to 
approach the required condition with reasonable 
efficiency, merely by omitting a central portion of 
the array, of length al, where the length of the 
whole array is 1. The choice of a is discussed below. 

The directivity function of the array with the 
central part omitted is that of the full array minus 
that of the central part alone, i. e. 


Dee sing sina _ 


are plotted, with a=0,15 by way of example, in 
Fig. 5. The zeros of A are marked with a dot, and 
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and D=A B=directional response of array with 
central part omitted (a=0.15). 


those of B with a cross. It is clear that if all secon- 
dary responses in the directivity function D (also 
plotted) are to be reversed in phase, then the zeros 
of A and B in each pair of zeros must coincide. The 
condition for this is evidently that 


DRS Sas ppeeee ye Delia 
hv ee Ey: l+a’ 
6x igs 
——— = ——_; etc., 
l-a l+a 


and it is clearly an impossible condition. But it is 

possible to eliminate unreversed polarity in any one 

of the regions concerned. For instance, by making 
Pach: 
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there will be completely reversed secondary respon- 
ses from |x| just below x to just above 4a, i.e. 
over the range included in Fig. 2 and 3. If this 
range includes all the range of real angles, then 
the desired result has been achieved. 

It would seem surprising without recourse to the 
method of Fig. 5 ,that the omission of exactly one- 
fifth of the array is needed to reverse the phase of 
the secondary lobes in the range of | 2| from 2x to 
3x. Fig. 6 curve 1, shows the height of the residual 
unreversed peak in this range if some value of a 
other than 0.2 is taken. Curve 2 shows the residual 
unreversed peak in the range of || from 4x to 5a. 
It is clear that if the minimum unreversed response 
is required up to |z|=6 2x, then a=0.15 is the best 
choice, giving no more than 1% unwanted response 
in the whole of this range; it can easily be shown 
that the maximum unwanted response anywhere is 
then about 3% in the region of |~|=10z. 

A formal approach to the synthesis of a directio- 
nal pattern with reversed-phase secondaries can be 
made by the method of superposing deflected 
(sin x) /x-patterns as described in [5], and it will 
be found that a reasonable approximation to the re- 
quired condition is given by a directional pattern 


sin © sin (a= 2.x); 


— 0.338 
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which is obtained by a sensitivity distribution along 
the length of the array given by 


1—0.776 cos2ra, 


where r is the proportion of the distance from the 
centre to the end of the array, i.e. r=0 at tha 
centre and r=] at the ends. This arrangement gives 
no positive response exceeding 1.25% anywhere 
from |x|=z to | |= co , but is, of course, much less 
convenient than the simple omission of a central 
section of the array. 


4. Conclusions 


It has been shown that in a limited type of appli- 
cation, greatly improved directional responses can 
be obtained by the use of synchronous demodulation 
in conjunction with an array whose secondary re- 
sponses have been caused to be all in reversed phase 
relative to the main beam. The general question of 
whether it is possible to obtain an array with such 
reversed-phase secondary responses over all possible 
angular directions has been discussed, and it ap- 
pears that it is only possible by the use of amplifiers 
or with great inefficiency; fortunately a close ap- 
proach to the condition can be obtained merely by 
omitting a central portion of the array. 

If the results obtained with this system are com- 
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: pared with those of the multiplicative system de- 
scribed in [1], it will be found that they give a 
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Fig. 6. Residual unreversed peak in relation to propor- 
tion (a) of array omitted. 


D. G. TUCKER: IMPROVED DIRECTIVITY 49 


rather better suppression of secondary responses 
but do not give so great a narrowing of the main 
lobe; moreover, the limitations on the field of use 


are greater. 
(Received 5th September, 1960.) 
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Buchbesprechung 


Leo L. Brranex, Noise Reduction, McGraw-Hill 
Book Company, New York 1960, 752 Seiten, tiber 
400 Bilder und zahlreiche Tabellen und Tafeln, 
16 cm x 24cm, Ganzleinen $14.50. 


Dieses neueste Buch des bekannten Autors faft ver- 
schiedene Vortragszyklen zusammen, die itiber das 
Thema ,,Larmbekampfung in der Industrie“, zwischen 
1953 und 1960 am Massachusetts Institute of Tech- 
nology gehalten wurden; daran waren verschiedene 
Dozenten von Hochschule und aus der Industrie betei- 
ligt, deren Beitrage nun von L. L. Beranex, der selber 
einen wesentlichen Anteil daran hatte, herausgegeben 
worden sind. Das Buch wendet sich ganz allgemein an 
Ingenieure und Physiker, setzt also keine besonderen 
akustischen Vorkenntnisse voraus; es enthalt daher als 
Einfiihrung einen ersten Teil ,,Schallwellen und ihre 
Messung™. 

Der mehr als 300 Seiten umfassende zweite Teil 
»Grundlagen der Larmbekémpfung“ ist das eigentliche 
Kernstiick des Buches; er ist in die folgenden Kapitel 
gegliedert: Schallausbreitung im Freien — Schall in 
kleinen Raumen — Schall in groSen Raumen — akusti- 
sche Eigenschaften pordser Materialien — Ubertragung 
und Abstrahlung von Schallwellen in undurchlassigen 
Bauelementen — Schalliibertragung in pordsen Stof- 
fen — Schallschluckstoffe im Bauwesen — _ Schall- 
dimpfer — Erschiitterungsisolierungen. Dieser Ab- 
schnitt enthaélt also neben der Schallausbreitung auch 
sehr vollstandige Grundlagen der Raumakustik, sowie 
die ganze Technik der Schall- und Erschiitterungsisolie- 
rung im weitesten Sinne. Es finden sich dabei nicht nur 


tibersichtliche und klare Darstellungen der theoreti- 
schen Beziehungen und Zusammenhange, sondern auch 
ein reiches Material von Tabellen und graphischen 
Darstellungen, die praktische Berechnungen und An- 
wendungen erleichtern oder sogar erst erméglichen. 


Den dritten Teil des Buches bildet eine Zusammen- 
stellung der besonders in den USA anerkannten und 
bentitzten Grenzwerte der zuldssigen Schallpegel, mit 
besonderer Beriicksichtigung der Verhiitung von Hor- 
schaden sowie fir Gebaude und Fahrzeuge verschiede- 
ner Arten. Der letzte, vierte Teil enthalt sodann aus- 
gewahlte Anwendungen, wobei besonders ausfiihrlich 
und vollstindig die Ventilationsanlagen zur Darstel- 
lung kommen, ferner auch typische Beispiele aus der 
Industrie sowie fiir Birohauser und Wohnungen. Es 
folgt schlieBlich noch ein wohldokumentierter und niitz- 
licher Abschnitt iiber das heute besonders aktuelle Pro- 
blem des Larmes von Diisentriebwerken und endlich 
uber die Isolierung von Eisenbahnwagen und Flugzeug- 
kabinen. 


Fur den wissenschaftlich gebildeten Fachmann ist das 
Buch ohne Zweifel eine Fundgrube zahlloser wertvoller 
Informationen, nicht zuletzt auch dank der ausfihr- 
lichen und durchwegs sorgfaltig redigierten Literatur- 
verzeichnisse. Fiir den kontinentalen Leser mag zwar 
das Nebeneinander, ja sogar Durcheinander englischer 
und metrischer Einheiten oft etwas befremdend und 
storend wirken; vielerorts vermit man auch eine star- 
kere Beriicksichtigung der europaischen Literatur. Der 
hohe Wert des Buches wird aber dadurch in keiner 
Weise beeintrachtigt. W. Furrer 


EINFLUSS VON STREUFLACHE 


UND HALLRAUMDIMENSIONEN AUF DEN GEMESSENEN 


SCHALLABSORPTIONSGRAD 
von U. Karn und W. Kun 


Institut fiir Rundfunktechnik, Hamburg 


Zusammenfassung 


Es werden die bei der Messung des Schallabsorptionsgrades im Hallraum mit und ohne 
Streukorper auftretenden zufalligen und systematischen Fehler untersucht. Dabei wird der 
Einflu8 der Form und Flache der Streukérper sowie der Leerraumabsorption auf diese 
Fehler bestimmt. Ferner werden Korrekturwerte fiir den gemessenen Schallabsorptionsgrad 
angegeben, die folgendes beriicksichtigen: 

1. die Schalldruckerhohung an Kanten und Ecken des Raumes, 

2. Die Abweichung der mittleren freien Weglinge von dem Wert 4V/S,, 

3. das unterschiedliche Absorptionsvermogen bei der Anbringung des absorbierenden 

Materials auf verschiedenen Flachen eines quaderformigen Raumes. 

Bei vollstandiger Belegung einer Raumbegrenzungsflache mit Absorptionsmaterial stimmt 
der im Hallraum gemessene und korrigierte Absorptionsgrad mit dem aus der Impedanz 
berechneten Absorptionsgrad fiir statistischen Schalleinfall innerhalb der Me8genauigkeit 
von wenigen Prozent tiberein. 


Summary 


Incidental and systematic errors occuring in measurements of the sound absorption co- 
efficient in a reverberation room with and without diffusors are investigated. The influence 
of shape and size of these diffusors and of the absorption area of the empty room on these 
errors is determined. Correction terms for the measured sound absorption coefficient are 
given, taking into account: 

1. the increase of sound pressure at edges and in corners of the reverberation room, 

2. the deviation of the mean free path from the value 4V/S, , 

3. the different effective absorptive capacity when the absorbing material is mounted on 

different boundaries of a rectangular room. 

Agreement is found within the accuracy of measurement of a few per cent between the 
measured and corrected sound absorption coefficients and the absorption calculated for 
random incidence, if one of the boundaries is completely covered with absorbing material. 


Sommaire 


On étudie les erreurs accidentelles et systématiques qui se produisent dans les mesures 
du coefficient d’absorption sonore dans une chambre réverbérante quand elle contient des 
corps diffusants ou quand elle n’en contient pas. Pour cela, on détermine l’influence sur ces 
erreurs de la forme et de la surface des corps diffusants, et celle de l’absorption de la 
chambre vide. On donne de plus les corrections 4 apporter aux mesures du coefficient 
d’absorption quand on tient compte des conditions suivantes : 

1. L’augmentation de pression sonore aux arétes et aux coins de la chambre. 

2. L’écart du libre parcours moyen par rapport a la valeur 4 V/S, . 

3. Le pouvoir d’absorption du matérial absorbant variant suivant sa disposition sur les 

différentes parois d’une chambre en forme de cube. 

En revétant complétement les parois d’une chambre réverbérante avec un matériau 
absorbant on trouve que le coefficient d’absorption mesuré puis corrigé coincide a quelques 
pour cents pres avec celui qu’on déduit de l’impédance dans le cas de l’incidence d’une onde 
sonore a distribution statistique. 


1. Einleitung 


Bei der Messung der Schallabsorption im Hall- 
raum treten bekanntlich bei unzureichender Diffusi- 
tat Schwierigkeiten auf. Die zufalligen Meffehler 
sind grob, und der gemessene Absorptionsgrad ist 
bei mittleren und hohen Frequenzen systematisch 
kleiner als der ,,richtige“, d. h. der bei ausreichender 


Diffusitat gemessene Wert. Meyer [1] hat bereits 
1937 durch Messungen in einem Modellraum ge- 
zeigt, daB eine groBe Verbesserung durch Streuk6r- 
per erzielt werden kann. Die von ihm verwendeten 
halbzylindrischen Streukérper waren auf drei zu- 
einander senkrechten Flachen in unregelmabiger 
Anordnung so angebracht, daB sie eine Schallstreu- 
ung in allen Richtungen bewirkten. Bei Untersu- 
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chungen von Sapine [2] war dagegen die Streuung 
nur in einer Ebene besonders wirksam. Er benutzte 
den in amerikanischen Laboratorien tiblichen, recht 
wirksamen rotierenden Doppelfliigel, dessen Streu- 
flache S; (Vorderseite und Riickseite zusammen) 
gleich 23% der Hallraumbegrenzungsflachen Sy war. 


-Zusatzlich stellte er rings um das am Boden an- 


gebrachte Prifmaterial eine Reihe von senkrechten 
Zylindern mit einem Durchmesser von 30cm. auf. 
Obschon ihre Oberflache 22% von S, betrug, fand er 
eine nur geringe Senkung der zufalligen Fehler und 
praktisch keine Anderung des in Anwesenheit nur 
des Fliigels gemessenen Absorptionsgrades. Wenn 
die Streuung nicht nur in der Ebene des Priifmate- 
rials erfolgt ware und dazu noch in der gleichen 


' Ebene wie die durch die Fliigel bewirkte Streuung, 


sondern auch in den beiden dazu senkrechten Ebe- 
nen, ware die Verbesserung sicher groBer gewesen. 


Meyer und Kurrrurr [3], [4] (dort Bild 3 und 


5) und Venzxe [5] stellten fest, da8 eine Schall- 


t 
\ 


streuung in zwei zueinander senkrechten Ebenen 
auch dann nicht ausreicht, wenn die gesamten Streu- 
flache viel groBer als die nach [1] und nach den 
Ergebnissen der vorliegenden Arbeit ausreichende 
Streuflache ist. Nach Meyer und Kurrrurr kann man 
bei einer Streuung in allen Richtungen und bei aus- 
reichender Streuflache, sofern der Absorptionsgrad 
hoch ist (2 >0,9), nur oberhalb einer Grenzfrequenz 
fer richtige Absorptionswerte messen. Die Grenzfre- 
quenz hangt vom Volumen des Hallraumes ab und 
ergibt sich aus der Naherungsformel 


VV fer>1mkHz. 


Die Autoren geben an, dai die Verhaltnisse bei 
geringer Absorption giinstiger seien. Sie diskutieren 
aber nicht die Tatsache, das man nach einigen ihrer 
Untersuchungsergebnisse bei sehr tiefen Frequenzen 
unabhaéngig vom Absorptionsgrad auch ohne Streu- 
korper keine systematisch zu kleinen Absorptions- 
gradwerte erhalt. Das ist aber qualitativ auch von 
anderen Untersuchungen her bekannt [6], [7]. Wir 
sind bei unseren Untersuchungen in einem Hall- 
raum normaler Grofe naher auf diesen Punkt ein- 
gegangen. 

Wenn man die Dicke von 2,5 beziehungsweise 
5,0 cm des im Modellraum [3], [4] verwendeten 
Absorptionsmaterials mafstabgerecht in einen Hall- 
raum tibertragen wurde, kame man in unserem Fall 
(V=125m?) zu Dicken von 16 beziehungsweise 
32 cm, bei dem von der ISO empfohlenen Volumen 
von 180 m* auf 19 beziehungsweise 38 cm. Solche 
Dicken kommen in der Praxis fast nie vor. Wir 
haben unsere Untersuchungen daher mit Mineral- 
faserplatten von 5 cm Dicke ausgefiihrt. 

Im einzelnen haben wir versucht, folgende Fragen 
zu klaren: 
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i. Wie gro sind die zufalligen und systemati- 
schen Fehler bei Nachhallmessungen und bei der 
Messung des Absorptionsgrades bei verschiedenen 
Streukorpern in Abhangigkeit von der prozentualen, 
das heift auf Sy bezogenen Streuflache? 

ii. Sind durch die Kombination verschiedenarti- 
ger Streukorper Verbesserungen zu erzielen? 

iii. Gibt es auch bei Absorptionsmaterialien mit 
einer fiir die Praxis interessierenden Dicke eine 
untere Grenzfrequenz fur die richtige Messung? 

iv. Welche Ursache haben systematische Abwei- 
chungen des richtigen MeBwertes von dem zu er- 
wartenden Absorptionsgrad, wie sie sich zuerst bei 
unseren Messungen klar zeigten, aber auch bei ande- 
ren neueren Untersuchungen zu finden sind [7], 


[8], [9]? 


2. Ubersicht der wichtigsten verwendeten 
Bezeichnungen 


A ____ Absorptionsflache der Probe, 


Ay _ Absorptionsflache des leeren Hallraums, 
A, — Absorptionsflache des Hallraums mit Probe, 
Ay, Absorptionsflache der Lautsprecher, 


Gy | Absorptionsgrad der Hallraumoberflache nach 
Berticksichtigung der Lautsprecher-Absorption, 
%,  Absorptionsgrad nach Eyrine bei ,,ausreichen- 
der“ Diffusitat (Mittelwert) , 
ay »richtiger“ Absorptionsgrad, das heibt ap 
nach Beriicksichtigung aller Korrekturen, 
a;  Absorptionsgrad nach Sasine, 
statistischer Absorptionsgrad, aus der Impe- 
danz berechnet, 
prozentuale Standardabweichung, 
Hallraumabmessung senkrecht zu S, , 
mittlere freie Weglange, 
mittlerer Fehler des Mittelwertes, 
prozentualer, mittlerer Fehler des Mittelwertes, 
berechnete Anzahl der Schallreflexionen, 
erwartete Anzahl der Schallreflexionen, 
Probenflache, 
zutalliger Fehler der Einzelmessung, 
- Hallraumoberflache, 
Raumbegrenzungsflache, auf der die Probe 
angebracht ist, 
Streuflache des Fliigels, 
Hallraumoberflache + Streuflache, 
3 Streuflache der Platten, 
S,  Gesamtstreuflache, 
T) Nachhallzeit im leeren Hallraum, 
T,  Nachhallzeit im Hallraum mit Probe, 
V Hallraumvolumen. 


eo 


nh 


nnn NN ns = 


3. MeBraum 
Der Grundrif8 des IRT-Hallraumes und seine 


Abmessungen sind aus Bild 1 zu ersehen. Bei einer 
Hohe von 2,9 m betragen das Volumen V = 125 m® 
und die Summe aller Raumbegrenzungsflachen 
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Bild 1. Grundri& des IRT-Hallraumes. 


Sy = 164 m?. Alle Flachen sind lackiert, Boden und 
Wande sind eben. Die Decke besteht in der Mitte 
aus sechs konkaven, flachen Zylinderabschnitten * 
und an den Stirnseiten des Raumes aus je einer 
0,85 m breiten, ebenen Flache. Als streuende Kan- 
ten konnen noch Mauervorspriinge auf beiden Langs- 
wanden gewertet werden, die 0,8 m vor der mit dem 
Priifmaterial bedeckten, kleinsten Flache des Rau- 
mes die Querabmessung des Hallraumes um 0,4.m 
verkleinern. Die Tir aus doppeltem Stahlblech be- 
findet sich gegeniiber der belegten Flache in einer 
0,2 m tiefen Nische. Durch den Raum fihren zwei 
isolierte Heizungsrohre von 10cm Durchmesser. In 
einer Ecke ist ein Lautsprecherkasten aufgestellt, des- 
sen Gehduse ein Volumen von 0,4 m® hat. Die Mem- 
bran hat eine betrachtliche Schallabsorption bei tie- 
fen Frequenzen !*. AuBerdem sind an der Decke zwei 
Absorber angebracht, die auf Frequenzen unter 
100 Hz abgestimmt sind. Durch diese drei Schall- 
absorber wurde die Streuung der Nachhallzeiten ver- 
ringert, die ohne sie ziemlich gro ware. Die Verklei- 
nerung des Absolutwertes der Nachhallzeiten, die 
dabei eintritt, kann man ohne Nachteil in Kauf neh- 
men. 

Der mittlere Absorptionsgrad a» der Begren- 
zungsflachen des Hallraums ist in Bild 2 tber der 
Frequenz aufgetragen; dabei wurde die Absorp- 
tionsflache der Lautsprecher A; 1* (vergleiche Ta- 
belle 1) bei der Berechnung von a,» beriicksichtigt. 


O10 


a 8 10? 2 4 


Bild 2. Absorptionsgrad der Hallraumbegrenzungsfla- 
chen in Abhangigkeit von der Frequenz. 


1 Sehnenlange der Bogen 0,8m. Stichmaf 6 cm. 
Achsenrichtung senkrecht zur Langsausdehnung des 
Raumes. Zwei Zylinderflachen werden jeweils durch 
einen 12 cm breiten, ebenen Streifen getrennt. 
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Tabelle I. 


Absorptionsflache A;, der Lautsprecher. 


f Ai ae f Ay 
Hz m? / kHz m? 
50 0,1 | 0,16 0,3 
64 1,0 | 0,2 0,2 
80 1,47 0,25—1,26 0,5 
100 0,77 16 BG 0,4 
126 0,45 6,4 —10,0 0,3 


Das scheinbare Anwachsen von dq bei Frequenzen 
uber 1 kHz ist durch die Schallabsorption der Luft 
verursacht. Bei der Berechnung des Absorptions- 
grads des Priifmaterials wurde der Absorptionsgrad 
der bedeckten Wand beriicksichtigt, und zwar bei 
hohen Frequenzen der gestrichelte extrapolierte 


Wert. 


4. Streukorper 


Die Streukorper, deren Einflu8 auf die Mebf- 
genauigkeit und auf den richtigen MeSwert wir 
untersuchten, waren einerseits gebogene Platten, 
andererseits ein groer Doppelfligel. 


Die gebogenen Platten aus Eisenblech von 1 mm 
Starke hatten die GroBe 1,0m>x 1,5 m. Sie wurden 
auf einer Seite mit einem dampfenden Belag ver- 
sehen, danach lackiert und um ihre kurze Achse ver- 
schieden stark gebogen. Bei den Messungen A, bis 
Ajo (siehe Tabelle II) hangten wir drei, neun, be- 
ziehungsweise zwolf dieser ,,Diffusoren* regellos 
im Hallraum auf. Bei den tibrigen Messungen kamen 
sie stets wieder an ihren urspriinglichen Platz. 


Bei den Me8reihen B und C erhéhten wir die 
Streuflache S, durch einen zusatzlich angebrachten, 
aus zwei ebenen Platten bestehenden Fliigel. Die 
beiden Fliigelflachen (je 2,0m x 2,0m, Streuflache 
S:=16m?) waren um 90° gegeneinander und je- 
weils 60° gegen die Drehachse geneigt. 


Es wurde sowohl mit ruhendem Fliigel (Mef- 
reihe B) als auch mit rotierendem Fliigel (fiinf 
Umdrehungen je Minute) gemessen (MeBreihe C). 


1a Das untere Maximum der Lautsprecherabsorp- 
tionsflache bei 80 Hz ist durch den Tieftonlautsprecher 
in der Ecke verursacht. Es wurde gemessen als die 
Schwingspule, wie immer wahrend des Nachhallvorgan- 
ges, offen war. Wenn die Fremderregung des Lautspre- 
chers abgeschaltet wird, verschwindet die Schallabsorp- 
tion durch diesen Lautsprecher fast vollstaéndig. Sie 
kann nur durch die Annahme erklart werden, da8 die 
Schwingspule eine KurzschluBwindung hat. Wenn die 
Schwingspule auBerdem auS8en durch den Ausgangs- 
transformator des Verstirkers geschlossen wird, erhoht 
sich die Absorptionsflache im Maximum auf 1,8 m?. 
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5. MeBreihen 


Unter einer Messung wollen wir alle Mefwerte 
fiir eine Anordnung der Streuflachen in dem Fre- 
quenzbereich von 80Hz bis 10kHz (22 Terz- 
bereiche) verstehen”. Fiinf derartige Messungen 


_bilden eine Mef{reihe*. In Tabelle II sind die an- 


gebrachten Streukérper, ihre Streuflache (Vorder- 
und Rtickseite zusammen) und die Bezeichnungen 
der entsprechenden Messungen zusammengestellt. 


Tabelle II. 


Bezeichnung der Messungen und Grof8e der Streuflache 


Se . Se = Sp + St . 
Gebogene Platten Fligel, Sy = 16 m? 

TO. Sp ohne ruhend _rotierend 
0 0 m? Ay By Cy 
3 9 m? A, B, C; 
6 18 m2 Ag B, Cs 
9 27 m* Ay By O5 

12 36 m? Ass By Cy. 


6. MeBanordnung 


Das Priifmaterial bestand aus Mineralfaserplatten 
Sillan 1004 (100 kg/m?) von 5 cm Dicke, mit denen 
die kleinste Hallraumwand mit einer Flache von 
16,1 m?, das hei®t 10% der Hallraumoberflache, 
vollstandig verkleidet war. Die absorbierenden 
Platten wurden nicht angeklebt, sondern durch eine 


_ Verspannung zwischen Decke und Boden des Hall- 


raumes fest an die Wand gedriickt. 


Als Priifschall diente Rauschen in Oktavband- 
breite, hinter dem Mikrophon wurden 1/3-Oktav- 
(Terz-) Siebe verwendet. Fiir Frequenzen f < 500 Hz 
wurde ein Lautsprecher mit nur einem System in 
einer Ecke benutzt, von 500Hz bis 10kHz ein 
Kugellautsprecher mit zwolf Systemen, etwa 1,2 m 
entfernt von einer anderen Ecke. Drei dynamische 
Mikrophone waren bei den Messungen bis 500 Hz 
in Raumecken angeordnet. Bei den hoheren Fre- 
quenzen hingen die Mikrophone in verschiedener 
Hohe und unterschiedlichem Wandabstand im Raum. 
Mit jedem Mikrophon nahmen wir drei (f < 100 Hz) 
beziehungsweise zwei (f 2 100 Hz) Nachhallkurven 
auf. Die Raumtemperatur (19,0 bis 22,0 °C) tund 
die relative Luftfeuchtigkeit (61 bis 66%) wurden 


2 Die Ergebnisse der bei den Mittelfrequenzen 50 Hz 
und 64 Hz durchgefiihrten Messungen streuten so stark, 
da8 wir sie nicht mitteilen. 

3 Die MeBreihe C mit rotierendem Fligel reicht nur 
yon 80 Hz bis 500 Hz. 

4 Der Firma Griinzweig & Hartmann AG., Ludwigs- 
hafen, die uns freundlicherweise das Priifmaterial fiir 
diese Messungen kostenlos zur Verfiigung stellte, mdch- 
ten wir bestens danken. 
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bei den Messungen bei f = 1kHz méglichst kon- 
stant gehalten. Zwischen der Leerraummessung und 
der entsprechenden Messung mit Probe anderte sich 
die Temperatur um weniger als 2 °C, die relative 
Luftfeuchtigkeit um weniger als 2%. 


7. Zur Frage der Absorptionsgradbestimmung 


Der Absorptionsgrad wurde zunachst nach der 
Eyrineschen Nachhallformel berechnet. Es erhebt 
sich die Frage, welcher der bei verschiedenen Mes- 
sungen erhaltenen Werte der richtige ist, auf den 
man die anderen Mefiwerte bei der Feststellung der 
systematischen Abweichungen beziehen kann. Mit 
wachsender Streuflache und damit zunehmender 
Diffusitat des Schallfeldes streben die Mefwerte 
einem oberen Grenzwert zu [3], [4]. 


Eine Ausnahme scheinen nur die Werte fir tiefe 
Frequenzen zu machen. Auf diese Erscheinung wer- 
den wir unten (9.1) noch naher eingehen. 


Es zeigt sich, dai die zufalligen Fehler bei den 
Messungen Ay, Ayo, By, Byy, Cy und Cys am klein- 
sten sind (vergleiche 8.1 und 8.2) und daf die Er- 
gebnisse dieser Messungen nicht systematisch von- 
einander abweichen. . 


Der Mittelwert der jeweils sechs MeBwerte fiir den 
Absorptionsgrad wird mit ap bezeichnet und einst- 
weilen als Bezugswert benutzt. a, weicht, wie wir 
spater zeigen werden, systematisch von dem rich- 
tigen, nach Anbringung mehrerer Korrekturen er- 
haltenen Absorptionsgrad a, ab, der mit guter 
Naherung gleich dem theoretischen Absorptions- 
erad fur statistischen Schalleinfall a;4 ist. Letzteren 
kann man aus der in einem Rohr gemessenen Im- 
pedanz mittels des Diagrammes von Morse und 
Botr [10] berechnen. Die Differenz von ap und a, 
beziehungsweise dsta, hangt, wie wir spater noch 
zeigen werden, von der Anordnung des Materials 
und von den Dimensionen des Hallraumes ab. Bei 
unseren Messungen (Priifmaterial auf der kleinsten 
Wandflache) war ay fiir alle Frequenzen systema- 
tisch groBer als ata, obschon Beugungserscheinun- 
gen durch die vollstandige Bedeckung einer Wand 
ausgeschlossen waren. Wir haben daher die ein- 
zelnen Mefiwerte auf den fiir die von uns gewahlte 
Anordnung geltenden Wert a, bezogen und nicht 
wie Meyer und Kurrrurr [3], [4] auf ota. 


In Bild 3 sind folgende Grofen tiber der Frequenz 
aufgetragen: d,,, wurde aus der Impedanz berech- 
net (Kurvea). Der ohne Streuk6rper gemessene 
Absorptionsgrad (Kurve b) weicht nicht so stark 
von ap ab, wie es Meyer und Kurrrurr [3], [4] in 
einem exakten Quaderraum feststellten, und zwar 
infolge der in unserem Fall auch ohne zusatzliche 
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ie | Die Bilder 4 und 5 zeigen das Absinken des pro- 
10 zentualen Fehlers der Nachhallzeiten dp mit der Fre- 
quenz, und zwar bei der Messung in leeren (Bild 4) 
oe beziehungsweise mit Priifmaterial (Bild 5) belegten 
| 06 Hallraumen ohne Streukorper und mit Streukorpern 
a in verschiedenen Kombinationen. Wir verzichten 
ui auf die Wiedergabe der Fehlerwerte der einzelnen 
02 an 
VAC) ra oe 
0 2 3 ae aA Yo 
6 810 2 4 pak 2 4 6 Hz10 100 
Bild 3. Absorptionsgrad von Sillan 100 (Dicke 5 cm) in 80 
Abhangigkeit von der Frequenz. t 
(a) dstat, aus der Impedanz berechnet, 5 60 
(b), (cy), (cg) Priifmaterial auf der kleinsten T ah 
Raumbegrenzungsflache. f 
(d) Prifmaterial auf dem Boden, 20 
(b) Messung ohne Streukorper, : 
fey), (eo), se mit Streukérpern aus- O. rs : ; oy ; mac 


(c,) Absorptionsgrad nach Sasine, 
(cy), (d) Absorptionsgrad nach Eyrine. 


Streukorper bereits vorhandenen Schallstreuung. 
Der Absorptionsgrad ware bei der Auswertung nach 
der Sasrne-Nachhallformel bis zu 8%  gro8er 
(Kurve c,) als bei der von uns zur Auswertung be- 
nutzten Formel von Eyrine (Kurve cy). Auf die 
Kurve d, die den gemessenen Absorptionsgrad des 
gleichen Materials bei der Anbringung auf dem 
Hallraumboden wiedergibt, gehen wir unten naher 
ein (9.4 FuBnote 11). 


8. Zufallige Fehler der Mefergebnisse 


8.1. Zufdllige Fehler der Nachhallzeiten 


Da die Angabe der zufalligen Fehler in der Lite- 
ratur nicht einheitlich ist, werden hier zwei mittlere 
Fehler unterschieden: 

i. die Standardabweichung [11] beziehungsweise 
mittlere quadratische Abweichung, also der zufallige 
Fehler der Einzelmessung 


Soy 
sat \ Bue und (1) 
(Minh 
ii. der mittlere Fehler des Mittelwertes 
a ae s 
M=+ 1/- ee Ns ee 2 
\ N(N—1) VN (2) 


Dabei bedeuten V Anzahl der Einzelmessungen, x 
Mefwert, Z arithmetischer Mittelwert, m= (%—72) 


Abweichung des Einzelwertes vom Mittelwert, 


6=* 100 prozentuale Standardabweichung, 
z 


ji M109 prozentualer, mittlerer Fehler des Mit- 
E 


telwertes. 
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Bild 4. Prozentuale Standardabweichung der Nachhall- 
zeit in leeren Hallrdumen in Abhangigkeit von 
der Frequenz. : 

(a) bis (d) IRT-Hallraum, 

(a) nur Platten (Mittel von Ag und Aj»), 

(b) Platten und ruhender Fliigel (Mittel von 
By und By»), 

(c) Platten und rotierender Fliigel (Mittel von 
Cy und Cys), 

(d) ohne Streukérper (Ay), 

(e) Hallraum der PTB ohne zusiatzliche Streu- 
korper nach [5], 

(f) Mittel der Hallraume 2 und 3 nach [5], 

(g) Hallraum der PTB mit zusitzlichen Streu- 
korpern. 


0 2 
a1—O2tSt”té«S 


6 810° 2 4 
f —— 


Bild 5. Prozentuale Standardabweichung der Nachhall- 
zeit in Hallréumen mit Priifmaterial (Sillan 100) 
in Abhingigkeit von der Frequenz. (a) bis (g) 
siehe Bild 4. 


Messungen und beschranken uns auf die Darstellung 
von drei arithmetischen Mittelwerten der Fehler. 
Diese Mittelwerte wurden jeweils aus den beiden 
Messungen mit maximaler Streuflaiche einer jeden 
MeBreihe, das heiBt aus den Messungen Ag, Ajo, 
By, By. und Cy, Cy. gewonnen, wobei alle auf- 
genommenen Nachhallkurven bei der Berechnung 
der Fehler beriicksichtigt wurden. Um eine Ver- 


Hie 
| 


j 
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gleichsmoglichkeit mit anderen Hallraumen*®® zu 
bekommen, wurden die prozentualen Fehler dy 
nochmals innerhalb einer Oktave arithmetisch ge- 
mittelt. 

Es ist zu erkennen, daB im Leerraum (Bild 4) bei 


hohen Frequenzen (f = 500Hz) durch die Streu- 


_korper nur eine geringe Verkleinerung der Fehler 


zu erreichen ist, da die Nachhallzeiten in diesem 
Frequenzbereich im IRT-Hallraum schon ohne Streu- 
korper eine geringe Streuung aufweisen. Fiir tiefe 
Frequenzen (f <500Hz) ist die Verbesserung der 


MeBgenauigkeit gegenitiber unserem Raum ohne 


Streukérper und den anderen Hallraéumen eindeutig. 
Dabei sind die gebogenen Platten allein giinstiger 
als zusammen mit dem ruhenden Fligel. Der rotie- 


 rende Fliigel verstarkt in dem Frequenzbereich von 


250 bis 500 Hz die giinstige Wirkung der Platten, 
weil deren Streuflache etwas zu klein ist (vergleiche 
Bild 6). 

In Bild 5 ist der giinstige Einflu8 der Streukor- 
per auf die prozentuale Standardabweichung der 
Nachhallzeiten im ganzen Frequenzbereich deutli- 
cher, da hier durch das einseitige Belegen des Hall- 
raumes mit Prifmaterial die Diffusitat des Schall- 
feldes stark gestért ist. Fir hohe Frequenzen zeigt 
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_ Bild 6. Mittelwert der prozentualen Standardabweichun- 


gen der Nachhallzeiten bei tiefen Frequenzen 
(f=100 bis 500 Hz) in Abhdngigkeit von der 
Streuflache bzw. von der auf die Hallraumober- 
flache bezogenen prozentualen Streufliche. 

(a,) bis (cy) leerer Hallraum, 

(ay) bis (cp) Hallraum mit Priifmaterial, 

(a) nur Platten, 

(b) Platten und ruhender Fligel, 

(c) Platten und rotierender Fliigel. 


® In dem von Venzxe [5] untersuchten Hallraum der 
(V =258 m3) 


_ sind die Decke und zwei Wande vollstandig mit Zylin- 


derabschnitten bedeckt, deren Achsen jedoch nur in 


' zwei zueinander senkrechten Richtungen verlaufen. Die 
- yon dem Autor beschriebenen Hallraume 2 (V =112 m?) 


und 3 (V=83m*) kénnen am ehesten mit dem IRT- 
Hallraum ohne zusatzliche Streukérper verglichen wer- 


_ den. Beide Réume sind rechteckig. Wahrend Raum 3 
- ebene Flachen hat, ist Raum 2 auf einer Seite vollstin- 


dig, auf einer anderen nur zum Teil mit Zylinder- 


abschnitten ausgestattet. Entgegen unserer Anordnung 


in einer Flache waren bei den Messungen von VENZKE 
¢ die Proben in drei Teilflachen angeordnet, wodurch die 


_ Diffusitat auch ohne zusitzliche Streukérper oder mit 


f N . . 
nur wenigen verbessert ist. 


| 
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sich wieder eine Gleichwertigkeit der Streuanord- 
nungen. Bei tiefen Frequenzen reicht die Streuflache 
der gebogenen Platten abermals nicht aus. Sowohl 
das Hinzufiigen des ruhenden Fliigels als auch des 
bewegten F'liigels verbessert die MeBgenauigkeit. 

Die Abhangigkeit der prozentualen Standard- 
abweichung der Nachhallzeit von der Streuflache S, 
beziehungsweise von der auf die Hallraumoberflache 
Sy bezogenen, prozentualen Streuflache wird in den 
Bildern 6 und 7 gezeigt. Dabei wurden die Fehler 
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Bild 7. Mittelwert der prozentualen Standardabwei- 


chungen der Nachhallzeiten bei hohen Frequen- 
zen (500 Hz bis 10kHz) in Abhangigkeit von 
der Streuflache beziehungsweise der auf die 
Hallraumoberflache bezogenen _ prozentualen 
Streuflache (a,) bis (b,) und (ay) bis (bg) 
siehe Bild 6. 


fiir alle dreiBig Messungen (iiber 3700 MeBwerte) 
innerhalb der beiden Frequenzbereiche 100 bis 
500 Hz beziehungsweise 500 Hz bis 10 kHz arithme- 
tisch gemittelt. 

Fiir die tiefen Frequenzen (f < 500 Hz) betragt 
der untere Grenzwert der prozentualen Standard- 
abweichungen im Leerraum 3%. Diesen Wert erhal- 
ten wir bereits mit einer prozentualen Streuflache 
von 16% bei den Anordnungen A und C. Der 
ruhende Fligel verschlechtert diesen Wert. Dagegen 
sind alle Streuanordnungen fiir hohe Frequenzen im 
Leerraum gleichwertig. Sie bringen kaum eine Ver- 
besserung (vergleiche auch Bild 4). 


6 Die Werte fiir die Kurven g in Bild 4 beziehungs- 
weise Bild 5 errechneten wir aus den Mefsprotokollen zu 
[9], die uns Herr Dr. Venzxe freundlicherweise zur Ver- 
fiigung stellte. Die Kurve g in Bild 4 gibt den arithmeti- 
schen Mittelwert von zwei Messungen wieder, die in 
Bild 5 von drei Messungen, wobei jede Messung wieder 
fiir 100 Hz < f <1kHz aus neun, fiir f >1kHz aus 
sechs Nachhallmessungen errechnet wurde. Die als 
Streuflachen verwendeten 60 Kugelsegmente bewirken 
im Leerraum eine Verbesserung der Mefgenauigkeit 
zwischen 200 Hz und 640 Hz (Bild 4), im Hallraum mit 
Priifmaterial zwischen 400 Hz und 2,5kHz (Bild 5) 
gegeniiber den Kurven e in den Bildern 4 und 5, wobei 
die Kurve g fiir eine Priifflache von 12 m? Sillan 100 
gilt, das heiSt bei einer Belegung von 5% der Hallraum- 
oberflache. Bei tiefen Frequenzen ist kein EinfluB8 der 
verwendeten Streukérper auf. die Mefgenauigkeit zu 
erkennen, da vermutlich die Abmessungen dieser Streu- 
k6rper zu klein sind. 


56 


Eine Voraussetzung fur so kleine Meffehler ist 
unter anderem eine besondere Schablone zur Aus- 
wertung der Nachhallkurven, die eine Ablese- 
genauigkeit von 1% ermoglicht. 

Kurve a, in Bild 6 zeigt, daB die Streuflache der 
Platten im Raum mit Prifmaterial fiir tiefe Fre- 
quenzen nicht ausreicht, um den Grenzwert von 4% 
zu erreichen. Dazu waren 25% notwendig, von uns 
realisiert durch Hinzufiigen des gedrehten Fligels. 
Da der Fligel, wenn er rotiert, in einem weiten 
Bereich den gebogenen Platten gleichwertig ist 
(Kurve c), aber ungiinstiger als diese wirkt, wenn 
er nicht rotiert (Kurveb), wird man in Zukunft 
ganz auf einen solchen Fliigel verzichten (Bild 6 
und 7). 

Baracuanpran [12] hat auch den Einflu8 von 
Streukérpern auf die MeBgenauigkeit untersucht. 
Aus seinen Messungen lassen sich aber keine kon- 
kreten Schliisse ziehen, weil der Einflu8 der verschie- 
denen Streukorper, deren Flache nicht angegeben ist, 
mit Wobbeltonen festgestellt wurde. Dabei war die 
MeBgenauigkeit mit allen Streukorpern immer noch 
geringer als ohne Streukorper bei Verwendung von 
Rauschen. Der Vorteil von Rauschen, mit dem alle 
vorhandenen Eigenschwingungen des Hallraums 
wirklich angeregt werden, gegeniiber Wobbeltonen, 
die nur wenige reine Tone enthalten, ist damit er- 
neut bewiesen. 

Die von BatacHanpRAN mitgeteilten Werte fur die 
zufalligen Fehler lassen sich auSerdem mit unseren 
Angaben nicht vergleichen. Die prozentualen Stan- 
dardabweichungen dort gelten nicht fiir Einzelmes- 
sungen, sondern fiir zwanzig Mefreihen zu je 200 
automatisch addierten Mefiwerten. 


8.2. Zufdllige Fehler des Absorptionsgrades 


Eine GroBe Y, die eine Funktion von 2;= s; ist, 
wie 


Veja Spe hoc sa cee Apmis 


hat nach dem Gaussschen Fehlerfortpflanzungsgesetz 
einen Fehler s, 


Creare (se) : 
3 


Analog ergibt sich der Fehler M,, (vergleiche 8.1). 

Obgleich wir zur Bestimmung des Absorptions- 
grades der Mineralfaserplatten die Eyrinesche Nach- 
hallformel benutzten, haben wir zur Berechnung des 
zufalligen Fehlers des Absorptionsgrades die SaBrnz- 
sche Nachhallformel angewandt. Das ist nicht nur 
weniger zeitraubend, sondern wir glauben auch, daB 
man dadurch leichter diese Ergebnisse mit den 
Resultaten anderer Untersuchungen vergleichen 
kann. 
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Fir die Fehlerrechnung gilt also 


] 1 
SEK tee 
y ie Tr) * 


wobei K = 0,162 V/S bedeutet. 

Dabei wurde der Absorptionsgrad % der Hall- 
raumbegrenzungsflachen und dessen Fehler nicht 
beriicksichtigt. 

Wenn wir fiir alle Abmessungen des Hallraums 
den Fehler + 0,02 m zulassen und die Konstante K 
mit 0,162+0,001 ansetzen, erhalten wir die in 


Bild 8 und 9 dargestellten prozentualen Fehler des — 


Absorptionsgrades in Abhangigkeit von der Fre- 
quenz. 

Bild 8 zeigt den Fehler 5a, der sich aus dem zu- 
falligen Fehler einzelner Messungen der Nachhall- 
zeit ergibt. In Bild 9 sind wir dagegen bei der Be- 
rechnung des Fehlers “a von den mittleren Fehlern 
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Bild 8. Prozentuale Standardabweichung des Absorp- 
tionsgrades in Abhangigkeit von der Frequenz. 
(a) nur Platten (Ay), 
(co) nur rotierende Fliigel (Co), 
(cy2) Platten und rotierender Fliigel (C,.), 
(d) Vergleichsmessungen an Hallonit aach 

[13], 

(e) Wiederholungsstreuungen nach [14]. 
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Bild 9. Prozentualer, mittlerer Fehler des Absorptions- — 


gradmittelwertes bei je sechs Nachhallmessun- 
gen, a, Cy , Cyo Siehe Bild 8. 
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M (vergleiche 8.1.) der Mittelwerte von jeweils 
sechs Nachhallmessungen ausgegangen. Die Kombi- 
nation von ruhendem Fliigel und gebogenen Platten 
(Mefreihe B) wurde hier nicht beriicksichtigt, da 
sie hinsichtlich der zufalligen Fehler ungiinstiger 
als die anderen Anordnungen ist. In beiden Bildern 


haben die Kurven a (nur Platten) ein Maximum 


bei 200 Hz, die Kurven cy (bewegter Fligel ohne 
Platten) bei 160 Hz. Erst die Kombination beider 
Streukérper (Kurve cj.) ergibt fiir tiefe Frequen- 
zen ({<360Hz) einen relativ glatten Verlauf. 
Offenbar hat die bei tiefen Frequenzen sicher stark 
frequenzabhangige Schallstreuung bei den angegebe- 
nen Frequénzen, abhangig von den jeweiligen Ab- 
messungen der Diffusoren, ein Minimum. Bei der 
Kombination verschieden grofer Platten werden 
diese Minima ausgeglichen 7. 

Es ist demnach unginstig, allen Streukorpern die 
gleichen Abmessungen zu geben, sondern besser, 


_Korper mit verschiedenen Abmessungen zu kom- 


binieren. 

In Bild 8 sind zum Vergleich auch Ergebnisse frii- 
herer Untersuchungen eingetragen, bei dem das 
Absorptionsmaterial auf drei zueinander senkrech- 
ten Flachen des Hallraums angebracht war, wobei 
die Diffusitaét ohne Streukorper erheblich besser ist 
als bei der Anbringung auf nur einer Flache. Bei 


_ Ersenserc [13], der iiber Vergleichsmessungen in 


mehreren Laboratorien berichtete, handelt es sich 


um Vergleichsstreuungen, bei Kiarpor [14] um 


Wiederholungsstreuungen in einem Laboratorium, 
im letzteren Fall um den Mittelwert der Messungen 
an zwei verschiedenen Materialien. Durch die Streu- 


_ korper wird die MeBgenauigkeit trotz der ungiinsti- 


geren Anbringung des Priifmaterials erheblich ver- 
bessert, insbesondere bei mittleren und hohen Fre- 
quenzen. 

Fiir Messungen mit Wobbelténen an Stelle von 
Rauschen gab tbrigens Etsenserc fir 100 bezie- 
hungsweise 200 Hz einen Fehler von +37% bezie- 
hungsweise + 35% an. 


8.2.1. Einflu8B der Leerraum-Absorp- 
tion des Hallraums auf den zu- 
falligen Fehler des Absorptions- 
grades 


Aus Gl. (4) ist zu erkennen, daf} der zufallige 
Fehler hauptsachlich von dem Ausdruck 


(5) 


abhangt. 


_ Um diesen Fehler klein zu halten, mu8 man ver- 


suchen, die Forderung 7) >T7, zu verwirklichen, 
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das: hei®t die Absorptionsflache A der Probe muf 
gro gegen die des Leerraums Ay sein. Nach der 
Neufassung von DIN 52212 [15] und der ISO- 
Empfehlung [16] ist die Probengré8e mit engen 
Toleranzen. festgelegt, so da8 nur ein Senken der 
Leerraumabsorption eine Méglichkeit bietet, den zu- 
falligen Fehler zu verringern. Das ist bei tiefen Fre- 
quenzen durch starre Raumbegrenzungsflachen, bei 
hohen Frequenzen durch unporése, zum Beispiel 
lackierte Oberflachen zu erreichen. 

In Bild 10 ist das starke Anwachsen des Faktors 
»F“ beim Absinken des Verhaltnisses A/A) dar- 
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Bild 10. Faktor F zur Berechnung der prozentualen 
Standardabweichung des Absorptionsgrades 
aus dem Fehler der Nachhallzeit in Abhangig- 
keit von dem Verhiltnis der Absorptionsflachen 
von Probe und Leerraum. 6, =F 67, ; 
u= Oty/OTy . 


gestellt. Aus dieser Darstellung kann man leicht die 
prozentuale Standardabweichung 6a des Absorp- 
tionsgrades berechnen. Im Gegensatz zu 0 wird bei 
6a nur der EinfluB der Nachhallzeitenfehler 67 be- 
rucksichtigt (vergleiche Gl. (4) und (5)). 
Ist der Fehler 67, der Nachhallzeit bekannt, so er- 
gibt sich mit 
A=agS=A,—Ay (6) 


die prozentuale Standardabweichung des Absorp- 
tionsgrades 


a =F Ory, (7) 


wobei der Faktor F dem Bild 10 entnommen wird. 
Als Parameter ist noch u=<7 9/671 aufgetragen, das 
heiBt das Verhaltnis des zufalligen Fehlers der Nach- 
hallzeiten im Leerraum zu dem im Raum mit Probe. 


In Tabelle III sind die Mindestforderungen [15], 


7 Bei den Untersuchungen von Meyer und Kurrrurr 
[4] (dort Bild 13 und 14) zeigte sich auch, daf die 
Wirkung der Streukérper in Abhangigkeit von der Fre- 
quenz, und zwar in diesem Fall im Hinblick auf das 
Erreichen des ,,wahren“ Absorptionsgrades, sehr yon 
ihrer GroBe abhangt. Die Platten der GroSe IIT wirken 
schon bei 1 kHz relativ gut, bei 1,5 kHz weniger und 
bei 2 kHz wieder gut, wahrend die doppelt beziehungs- 
weise halb so grofen Platten II und IV eine glattere 
Frequenzkurve ergeben. 
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[16] der Nachhallzeiten fiir einen Hallraum mit dem 
empfohlenen Mindestvolumen von V = 180m? und 
die sich daraus ergebende Leerraumabsorptionsfla- 
che zusammengestellt. AuSerdem haben wir dort die 
nach Gl. (7) zu erwartende prozentuale Standard- 
abweichung des Absorptionsgrades fiir eine Probe 
von 12 m2 Sillan 100 berechnet, indem wir die Feh- 
lerwerte der Nachhallzeiten den Bildern 4 bezie- 
hungsweise 5 entnahmen, wobei die Frequenzen mit 
dem Faktor (125/180)'* umgerechnet wurden. Ob- 
gleich die Forderungen der Mindestnachhallzeiten 
nicht streng sind (vergleiche Bild 11), wurden diese 
Bedingungen nach Einbringen von Diffusoren jedoch 
von zehn der sechzehn am Ringversuch 1959/60 8 
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Bild 11. Nachhallzeiten im Leerraum (a,), (b,) und im 
Hallraum mit Priifmaterial (aj), (bs) in Ab- 
hangigkeit von der Frequenz. 

(a) ohne Streukorper, 
(b) mit der gréBten bei den Versuchen ver- 
wendeten Streuflache. 


beteiligten Laboratorien nicht erfillt. Kin weiterer 
Hallraum hatte ein Volumen von nur 86m, das 
hei®t er entspricht nicht der Bedingung, da8 beim 
Weeglassen der tiefsten Meffrequenz 125 Hz wenig- 
stens V = 100 m? [15], [16] sein muf. 

Nach unseren Berechnungen erscheint es sinnvol- 
ler, die gro®te zulassige Leerraumabsorptionsflache 
an Stelle der minimalen Nachhallzeit festzulegen. 


Tabelle III. 
Mindestnachhallzeiten fiir das in der Norm empfohlene 
Mindestvolumen V =180m° und die sich daraus fir 
den Absorptionsgrad (S=12 m? Sillan 100) ergebende 

prozentuale Standardabweichung. 
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8 Die Nachhallzeiten der bei dem zweiten Ringver- 
such 1959/60 benutzten leeren Hallraume mit Streu- 
k6rpern sowie die Volumina der Hallraume entnehmen 
‘ wir dem Bericht [17] von Herrn Prof. C. W. Kosten, 


den er uns dankenswerterweise zur Verfiigung stellte. 
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Damit wiirde man fir alle Hallraéume ungefahr die 
gleiche Fehlergrenze erhalten. Die Nachhallzeit 
mii®te dann mit dem Faktor (V/180) steigen bezie- 
hungsweise fallen. Betragt die GroBe des Faktors 
nur (V/180)'* beziehungsweise (V/180)**, so sind 
die Fehler bei den grofen Hallraumen gréfer als 
bei den kleinen Raumen ®. 


9. Systematische Fehlerquellen 
_ bei Hallraummessungen 


9.1. Systematische Fehler des Absorptionsgrades in- 
folge unzureichender Diffusitét des Schallfeldes 


Bekanntlich wird bei unzureichender Diffusitat 
des Schallfeldes der Absorptionsgrad bei hohen Fre- 
quenzen zu niedrig gemessen. Bei unseren Messun- 
gen an Sillan 100 ohne Streukorper ist das fiir 
f>250Hz der Fall beziehungsweise bei der Mes- 
sung von Absorptionsgraden a>0,55 (Bild 3, 
Kurve b). Gleichzeitig findet man fiir tiefe Frequen- 
zen. (bei uns fiir f< 250 Hz) einen zu hohen Absorp- 
tionsgrad. 

Bild 11 zeigt die Nachhallzeiten des IRT-Hallrau- 
mes in Abhangigkeit von der Frequenz. Fir den 
Leerraum und den Raum mit Priifmaterial sind nur 
die Grenzwerte der Nachhallzeiten dargestellt, ein- 
mal ohne Streukorper (a, , a), das andere Mal mit 
allen Streukérpern (b, beziehungsweise bj). Wah- 
rend die Nachhallzeit im Leerraum beim Einbringen 
der Streukorper sinkt, ist dieses im Raum mit Priif- 
material fiir tiefe Frequenzen nicht der Fall, das 
heif®t fiir tiefe Frequenzen ist der bei unzureichen- 
der Diffusitaét gemessene Absorptionsgrad grofer 
als der bei ,,voller“ Diffusitét gemessene Wert ap, 
wenn das Priifmaterial nur auf einer Begrenzungs- 
flache angebracht ist. 

Diese Erscheinung kann man in den Arbeiten meh- 
rerer Autoren bestiatigt finden. Zunachst in Bild 3 
der vorliegenden Arbeit (f<250Hz), in Bild 11 
von [7] (f<800Hz), in Bild 1 einer eigenen frii- 
heren Veroffentlichung [6] (f<150Hz), in Bild 2 
und Bild 11 von (4) sowie in dem Bericht von C. W. 
Kosren [17], in dem diese Tendenz bei zehn der 
fiinfzehn verschiedenen Hallraume zu finden ist. 

Kine Erklarung fiir diese Erscheinung kann nicht 
gegeben werden. Bei Hinzufiigen der Streukérper 
nimmt also ein Teil der Gesamtabsorption zu, der 
andere ab. 

Fiir die MeBreihe A haben wir in Bild 12 die pro- 
zentualen Gesamtfehler des Absorptionsgrades dar- 
gestellt, und zwar wurde der auf a, bezogene Abso- 


® Bei der Neuausgabe von DIN 52212 [15] wird 
jetzt in den Erlauterungen (dort 2.14) dieser Tatsache 
Rechnung getragen. 
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‘Bild 12. Prozentuale Absolut-Abweichung des Absorp- 
tionsgrades vom Bezugswert ap in Abhangig- 
keit von der Frequenz. Einfluf der Streufla- 
chen. Ay bis Ay» siehe Tabelle II. 


lutwert der Abweichung von ag, das heift die pro- 
_ zentuale Absolutabweichung in Abhangigkeit von 
der Frequenz aufgetragen. Der Bezugswert ap 
wurde, wie schon unter 7. erwahnt, aus sechs Mes- 
_ sungen gebildet. Der Gesamtfehler enthalt sowohl 
den systematischen Fehler auf Grund unzureichender 
Diffusitaét als auch die zufalligen Fehler, die sich 
nicht voneinander trennen lassen. Analog haben wir 
fiir die Mefireihen B und C die Fehler betrachtet. Da 
sich beide MeBreihen wie die MeBreihe A verhalten, 
verzichten wir auf die Darstellung dieser Ergebnisse. 
_ Bild 13 zeigt den prozentualen Gesamtfehler der 
Messungen in Abhangigkeit von der Streuflache 
beziehungsweise prozentualen Streuflache. Dabei 
wurden wieder fiir jede Mefireihe die Fehler inner- 
halb zweier Frequenzbereiche gemittelt. Die Kur- 
ven a,, b,, c, geben die gemittelten Fehlerwerte 
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Bild 13. Mittelwert der prozentualen Absolut-Abwei- 
chung des Absorptionsgrades vom Bezugswert 
ag in Abhangigkeit von der Streufliche bezie- 
hungsweise der auf die Hallraumoberflache 
bezogenen prozentualen Streuflache. 

(a,) bis (c;) 100 Hz bis 500 Hz (az) bis (bs) 
500 Hz bis 10 kHz. (a) bis (c) siehe Bild 6. 


fiir die MeBreihen A, B, C in dem Frequenzbereich 
_ von 100 bis 500 Hz an, analog gelten az, b, fiir den 
_ Frequenzbereich von 500 Hz bis 10 kHz. Wie man 
sieht, wiirde bei hohen Frequenzen schon eine Streu- 
 fliche von 27m, das heif&t 16% der Hallraum- 
- oberfliche geniigen, um einen Fehler < 2% zu er- 
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halten. Um bei tiefen Frequenzen einen Fehler <4% 
zu erhalten, benétigt man eine Streuflache = 21% 
der Hallraumoberflache. 

Wir finden also bei der Gesamtfehlerbetrachtung 
des Absorptionsgrades die gleichen Bedingungen 
fiir die GroBe der Streuflache wie unter 8.1. bei der 
Betrachtung der zufalligen Fehler der Nachhallzei- 
ten. Dadurch ist bewiesen, da der systematische 
Fehler, der durch unzureichende Diffusitat hervor- 
gerufen wurde, durch das Anbringen von Streu- 
korpern mit ausreichender Streuflache vollkommen 
ausgeschlossen wird. Gleichzeitig ist damit die Rich- 
tigkeit der Mittelwertbildung von ap bestatigt, da 
die Absorptionsgradfehler der zur Mittelwertbildung 
herangezogenen Messungen nur durch die zufalligen 
Fehler bestimmt werden. 


9.2. EinfluB der Hallraumdimensionen auf den 
systematischen Fehler des Absorptionsgrades 


Wie wir schon unter 7. erwahnt haben, ist der 
von uns gemessene Absorptionsgrad ag, wenn das 
Priifmaterial auf der kleinsten Wand angebracht ist, 
fiir alle Frequenzen systematisch grofer als der zu 
erwartende Wert a,. Fir hohe Frequenzen finden 
Wir sogar 4 >1,0, worauf wir schon hinwiesen 
[8], was aber im Widerspruch zu der bisherigen 
Auffassung (zum Beispiel [7], [18]) steht. Die 
Moglichkeit, Absorptionsgradwerte tiber 1 zu mes- 


sen, selbst wenn man eine geringe Erhohung durch 


Beugung rechnerisch beriicksichtigt, wurde in der 
letzten Zeit von Venzke [9] bestatigt. Unsere letzte 
Messung, bei der das Priifmaterial am Boden an- 
gebracht war (Kurved in Bild3), ergab einen 
scheinbar zu geringen Absorptionsgrad. Die Ursache 
fir die positive beziehungsweise negative Abwei- 
chung von dem zu erwartenden Wert sind nicht der 
MeSvorgang bzw. die Auswertung der Nachhallkur- 
ven, sondern u. a. die Anbringung des Priifmaterials 
auf verschiedenen Flachen eines quaderformigen 
Raumes. Das Verhiltnis der Dimensionen des Rau- 
mes spielt aufSerdem eine Rolle bei der Abweichung 
der mittleren freien Weglange /,, von dem gebrauch- 
lichen Wert 4 V/S,). Diese Einfliisse werden von kei- 
ner der bekannten Nachhallformeln beriicksichtigt. 


9.2.1. Einflu®8 der Hallraumdimensio- 
nen auf die mittlere freie Weg- 
lange 

Bei der Herleitung seiner Nachhallformel verwen- 
det Eyrine den Begriff der mittleren freien Weg- 
lange 1,, =4 V/S). Nach einer statistischen Methode 
haben Atirep und Newnovuse [19] mit Hilfe elek- 
tronischer Rechnung /,, fiir fiinfzehn verschiedene 

Verhaltnisse der Dimensionen eines quaderférmi- 

gen Raumes bestimmt. Dabei ergeben sich Abwei- 


Mg * ba Vy ary set 
piped iucinns: -~: 
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chungen bis zu 5% von dem Wert 4V/S). Indem 
wir die berechnete freie Weglange /, in Abhangig- 
keit von dem Verhaltnis des quadrierten geometri- 
schen Mittels zu dem arithmetischen Mittel der Hall- 
raumabmessungen, das hei®t tiber 3V"/x+y+z 
aufgetragen haben (Bild 14), ist es moglich, den 
genauen Wert l,, fiir jeden in der Praxis vorkom- 
menden Quaderraum einer Kurve zu entnehmen. 
Das Verhaltnis |, :4V/So laBt sich nicht in Form 
einer einzigen Kurve darstellen. Fiir verschiedene 
Verhaltnisse der Dimensionen erhalt man dann viel- 


7 


a3 


iz 3 4 5 6 dh 8 


fe 9m 10 
3(xYy2) X42 — 


Bild 14. Berechnete mittlere freie Weglange in Hall- 
raumen [19] (V=2-y-z) in Abhingigkeit von 
dem Verhaltnis des quadrierten geometrischen 
Mittels zu dem arithmetischen Mittel der Hall- 
raumabmessungen. 


mehr eine Kurvenschar. Fir den IRT-Hallraum er- 
gibt sich aus Bild 14 eine Weglange, die um 1,3% 
kleiner ist als 4V/S). Der entsprechend korrigierte 
Absorptionsgrad ist also um 1,3% kleiner als bei 
der Anwendung der Eyrinc-Formel. Nach unserer 
Meinung diirfte ein Einflu8 der Streukorper auf J, 
nicht vorhanden sein, da der Schall bei der An- 
wesenheit von Streukérpern zwar auf anderen 
Wegen lauft als ohne Streukérper, die Raumbegren- 
zungsflachen jedoch nicht haufiger getroffen werden. 
Mit /,, wird aber im Sinne der Eyrineschen Nach- 
hallformel der mittlere Weg des Schalls zwischen 
zwei aufeinanderfolgenden Reflexionen an den 
Raumbegrenzungsflichen bezeichnet. Wir haben 
daher bei der Berechnung des Absorptionsgrades ay 
in der Eyrineschen Nachhallformel den Wert 


0,162 V 


—_—— 8 
ee (8) 
benutzt und nicht 
0,162 V 
SETTER 8b 
Dosing GE 


wie es Kurrrurr [3] getan hat. Der Unterschied der 
beiden Rechnungen macht sich bei dem von uns 
untersuchten Priifmaterial nur bei hohen Frequen- 
zen starker bemerkbar. Die Absorptionsgradwerte, 
die nach Gl. (8 a) bestimmt werden, sind etwas klei- 
ner als die nach Gl. (8b) errechneten, und zwar bei 
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160 Hz um etwa 1%, bei 1 kHz um etwa 1,5% und 
bei 10 kHz um etwa 3%. 

Eine endgiiltige Klarung der Frage, ob die Dif- 
fusorenoberflache in die Nachhallformel eingeht 
oder nicht, mu wohl einer statistischen Untersu- 
chung mit Hilfe einer elektronischen Rechenmaschine 
iberlassen bleiben. Unsere MeBergebnisse lassen die 
Auswertung nach Gl. (8a) als richtiger erscheinen. 


9.2.2. EinfluB der Hallraumdimensio- 
nen auf das effektive Absorp- 
tionsvermogen auf den verschie- 
denen Hallraumoberflachen 


Autrep und Newnouss haben in ihrer Arbeit [19] 
auch die Zahl der Schallreflexionen an den einzelnen 
Raumbegrenzungsflachen fiir die fiinfzehn Dimen- 
sionsverhaltnisse berechnet, indem sie das Verhalt- 
nis n/n bildeten, wobei n, die nach der Flachen- 
groBe zu erwartende Anzahl der Reflexionen, n die 
berechnete Zahl bedeuten. : 

In Bild 15 haben wir nun den reziproken Wert 
dargestellt. Dabei bedeuten S, die Hallraumbegren- 
zungsflache, auf der das Priifmaterial angebracht ist, 
1, die Dimension des Raumes senkrecht zu der Flache 
S, 1°, Die stark gezeichneten Kurven geben die Werte 
fiir Hallraume wieder, fiir die die Raumdiagonale 
d<1,9 V') ist, wie es fiir die Hallraume nach [15], 
[16] gefordert wird. Die gestrichelt gezeichneten 
Kurven beziehungsweise Kurventeile gelten fiir sehr 
flache oder sehr langgestreckte Hallraume, die diese 
Bedingung nicht erfiillen. 

Wahrend man fiir die Werte l,///S,<1 die Ergeb- 
nisse der Rechnungen, die mit einer Genauigkeit von 
etwa 2% durchgefiihrt wurden, durch einen Kurven- 
zug darstellen kann, ergibt sich fir 1,VS,>1 eine 
Kurvenschar, die zur Interpolation zwingt. Die drei 
senkrechten Balken in Bild 15 geben die Grenzen 
des gesuchten Wertes bei einer zweidimensionalen 
Naherungsrechnung entsprechend einer ganz ande- 
ren, noch unverOffentlichten Betrachtung des: Pro- 
blems von WestpHaL und Kartu an. Die Werte 
l,/VSa=1,5 beziehungsweise 2,0 stimmen trotz ganz 
anderer Voraussetzungen gut mit der Kurvenschar 
tiberein. 

Die Punkte 1 und 2 geben die Korrekturwerte fiir 
unsere Messungen an. Bei der Anbringung des 
Absorptionsmaterials auf der kleinsten Flache des 
IRT-Hallraums ergibt sich also ein um 7,7% 
(Punkt 1) zu groBer MeSwert. Liegt dagegen das 
Priifmaterial auf dem Hallraumboden, so mi&t man 
einen um 5% zu geringen Absorptionsgrad (Punkt 2). 


10 Statt 1,/V/S_ kénnte man als Abszisse auch S,/6 Sa 
wiahlen, das heiBt das Verhaltnis des arithmetischen 
Mittelwertes einer einzelnen Raumbegrenzungsflache 
zu a « 


: 


i 
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Bild 15. Zur Korrektur des in quaderformigen Hall- 
raumen gemessenen Absorptionsgrades. 
Ordinate: Verhaltnis der berechneten Zahl (n) 

der Schallreflexion an einer Raum- 


begrenzungsflache [19] zu der 
nach ihrer Groéfe erwarteten An- 
zahl (ns). 


Abszisse: Verhaltnis der Hallraumabmessung 
(2,), die senkrecht zu der mit Prif- 
material belegten Raumbegren- 
zungsflache (S,) verlauft, zu dem 
geometrischen Mittel der Abmes- 


sungen dieser Flache. 
Stark ausgezogene Kurventeile gelten fiir Hall- 
eee 
raume, deren Raumdiagonale d< 1,9 VV ist, 


die diinn gezogenen Kurventeile gelten fir 
Raume, die diese genormte Bedingung nicht 


erfillen. 
Verhaltnis der 
Kurven Hallraumabmessungen 
(x: y:2) 
a 10:10:10; 8;6;4;2 
b 10: 8: 8;6;4;2 
e 102-6. 63/4; 2 
d 10: 4: 432 


n/n; in Abhangigkeit von dem Verhiltnis 1,///S, 
Die Indizes der Kurvenbezeichnungen gelten 
fir 1: 2/Vyz; 2: y/Vxz; 3: x/Vxy. Die 
senkrechten Balken geben die Grenzwerte der 
Ergebnisse einer zweidimensionalen Nahe- 
rungsrechnung von Wesrpxat und Karu an. 
Punkt 1 und 2: Reziprokwerte der Korrektur- 
faktoren fiir den IRT-Hallraum. Priifmaterial 
auf.der kleinsten (1) beziehungsweise groften 
(2) Hallraumflache. 

Punkt 3 bis 5: Verhaltnis des effektiven Ab- 
sorptionsgrades zu dstat. (Mittelwert fiir hohe 
Frequenzen bei Messungen an Sillan 100, An- 
derung der mittleren freien Weglinge rech- 
nerisch beriicksichtigt.) 

3 und 4; nach Meyer und Kurrrunr [4], 

5: nach Venzxe [9], Beugung nach [6] elimi- 
niert. 


Die Punkte 3—5 sind MeSpunkte verschiedener 
Autoren (vergleiche Legende Bild 15) und gelten 
fiir den nach Eyrine ausgewerteten Grenzwert des 
Absorptionsgrades fiir Frequenzen iiber 1 kHz bei 
ausreichender Diffusitat, wobei in einem Fall [9] 


7 
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(Punkt 5) die Beugung rechnerisch nach [6] elimi- 
niert und in jedem Fall die mittlere freie Weglange 
(vergleiche 9.2.1.) eingesetzt wurde. Die restlichen 
Abweichungen vom theoretischen Wert liegen unter 
1,5%. Bild 15 besagt also folgendes: Bringt man das 
Absorptionsmaterial auf dem Boden eines flachen 
Hallraumes an, so ist der gemessene scheinbare Ab- 
sorptionsgrad ap kleiner als o,¢a,. Beim Anbringen 
auf der Langswand wird er grofer und ist auf der 
kleinsten Wand, der Querwand, am grof&ten. Ob- 
gleich dieses Ergebnis von AttreD und NewHouse 
fur einen ,,diffusen“ Hallraum nicht ohne weiteres 
zu erklaren ist, haben wir doch durch Rechnung und 
MeB8ergebnisse eigener und fremder Untersuchungen 
die Giltigkeit innerhalb der MefSigenauigkeit be- 
statigt. 


9.3. EinfluB der Schalldruckerhohung an den Kan- 
ten und in den Ecken des Hallraumes auf den 
Absorptionsgrad 


Es ist bekannt, da8 fur Absorptionsmaterial un- 
abhangig von der Schallbeugung unterschiedlich 
groke Absorptionsgradwerte gemessen werden kon- 
nen, je nachdem ob das Prtfmaterial auf einer 
Wand, an einer Kante oder in einer Ecke des Hall- 
raumes angebracht ist. Diese Erscheinung, die auf 
die Schalldruckerhohung an den Kanten und in den 
Ecken eines Hallraumes zuriickzufuhren ist, wurde 
durch Rechnungen und Messungen von mehreren 
Verfassern [20] — [24], [6], [8] behandelt, sehr 
ausftihrlich von Warternouse [24]. Der Verfasser 
berechnet die Verhaltnisse p?/po?, v?/vo? und E/E, 
in Abhangigkeit von dem Wandabstand und besta- 
tigt seine Ergebnisse durch Messungen. py, vy und 
E, bedeuten Schalldruck, -schnelle beziehungsweise 
-energie in der Raummitte, p, v und EF die entspre- 
chenden Werte in der Nahe der Raumbegrenzungs- 
flachen. Da bei unseren Untersuchungen pordses 
Absorptionsmaterial verwendet wurde, interessiert 
hier nur der Verlauf von p?/po”. 

Bei harten Begrenzungsflachen ist der Schalldruck- 
pegel in einem ,,diffusen“ Schallfeld an den Wanden 
um 3 dB, an den Kanten um 6 dB und in den Ecken 
um 9dB gegenttber dem Pegel in der Raummitte 
erhoht. Wird eine Raumbegrenzungsflache vollstan- 
dig mit Absorptionsmaterial belegt, so wird der 
Schalldruck auf dieser Wand, abhangig von dem 
Absorptionsgrad des Priifmaterials, gegentiber dem 
auf der harten Wand sinken. Bei einer vollig absor- 
bierenden Wand betragt die Pegelerhohung gegen- 
uber der Raummitte an Wandkante und Ecke nur 
noch 0 dB, 3 dB beziehungsweise 6 dB. Die Pegel- 
differenz zwischen Kante und Wand, beziehungs- 
weise Ecke und Wand bleibt aber unabhangig von 
der Absorptionsfahigkeit des Priifmaterials 3 dB 
beziehungsweise 6dB. Die von Warernouse fir 
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starre Wande berechneten Werte lassen sich also bei 
allen Prifmaterialien und Frequenzen verwerten. 


Nach [24] gilt in einem diffusen Schallfeld fiir den 


Schalldruck mit po? = 1 und k=2 x/2 
ka 
1 sin 2 kx 
2—]4 : ee ey 9 
‘ eed Qhe ) 
Ky 


wobei x den Abstand senkrecht zur Raumbegren- 
zungsflache, k, und ky die Wellenzahlen der Grenz- 
frequenzen f{,, f, eines Frequenzbandes bedeuten. 
Die mittlere Erhohung von p? ergibt sich also in 
dem Bereich von 0 bis 2, durch 


iy (10 a) 


Wenn x, >A, kann die obere Grenze des Integrals 
| naherungsweise durch o ersetzt werden. Dann gilt 
zum Beispiel fiir ein ager le von Terzbreite 


hee [fee Pps oes 
4) con Siti 


(10b) 


wobei /,, die Wellenlange der Mittenfrequenz von 
f, und fy bedeutet. Der gegentiber der Hallraum- 
mitte entsprechend Gl. (9) erhdhte Schalldruck auf 
einer senkrecht zum Absorptionsmaterial stehenden 
Wand erhéht an der Kante auch den auf das Ab- 
sorptionsmaterial wirkenden Druck und damit die 
Schallabsorption, und zwar in einem mit der Wellen- 
lange wachsenden Bereich. Wenn der Schalldruck 
auf der belegten Wand, dessen Verhaltnis zu py vom 
Absorptionsgrad abhangt, in der Mitte («>4A,,/8) 
konstant und an der Kante in einem streifenformi- 
gen Bereich mit der Breite x(0 <2 <4,,/8) um 
den Faktor /2 erhéht ware, so hatte dieser einfache 
Druckverlauf auf die Erhohung der Absorption die 
gleiche Wirkung wie der tatsachliche, komplizierte 
Druckverlauf. 

An den Raumecken tuberdecken sich die beiden 
senkrecht aufeinanderstehenden Streifen. Diese klei- 
nen Flachen werden also doppelt gerechnet, womit 
die Erhdhung des Pegels um 6dB gegeniiber der 
Mitte der Wand und die Vergroéerung des Absorp- 
tionsgrades um den Faktor 4 (an Stelle des Fak- 
tors 2 an den Kanten) beriicksichtigt wird. Bei einer 
vollig bedeckten Raumbegrenzungsflache mu also 
die Schalldruckerhdhung und die dadurch verur- 
sachte scheinbare VergroBerung des Absorptions- 
grades durch eine aquivalente VergréBerung der 
Priifflache AS beriicksichtigt werden. Ist U der Um- 
fang der vollig mit Priifmaterial bedeckten Raum- 


begrenzungsflache, so ergibt sich AS durch 


AS=U 2 
8 


(11) 
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9.4. Korrektur der systematischen Fehler 

In Bild 16 ist die Differenz zwischen den gemesse- 
nen und statistisch zu erwartenden Absorptionsgrad- 
werten in Abhangigkeit von der Frequenz aufgetra- 
gen. Die Kurven a, und b, geben die Werte ap — atat 


0,16 


-0,04 


-0084, 
10 2 4 


L 
6 810° 2 4 
f —.> 


6 Hz 10* 


Bild 16. Abweichungen des gemessenen Absorptions- 
grades a von dem aus der Impedanz berech- 
neten Absorptionsgrad dgtat in Abhangigkeit 
von der Frequenz. Priifmaterial auf der klein- 
sten (Kurve (a)) beziehungsweise auf der 
groBten (Kurve (b)) Raumbegrenzungsflache 
angeordnet. 

(ay), (b;) MefSwerte ag nicht korrigiert, 

(as), (by) wie (a,), (by), jedoch Schalldruck- 
erhéhung an Kanten und KEcken rechnerisch 
beriicksichtigt. 

(ag), (bs) wie (as), (by), jedoch Abweichung 
der mittleren freien Weglinge von 4V/S und 
Absorptionsvermogen auf der belegten Raum- 
begrenzungsflache rechnerisch beriicksichtigt. 


der Messungen bei Befestigung des Absorptions- 
materials auf der kleinsten beziehungsweise grof- 
ten!! Raumbegrenzungsfliche des IRT-Hallraumes 


1 Bei dieser Messung wurde der Hallraumboden 
mit 42m? Priifmaterial véllig ausgelegt. Da wir mit 
einer Streuflache von Sp=52 m? bei hohen Frequenzen 
fiir den Absorptionsgrad zu niedrige Werte erhielten, 
erhéhten wir Sp auf 90 m?, das heiBt Sp/Sy=52%. Ver- 
groBert man die Streuflache im gleichen Verhaltnis wie 
die Priifflache, dann hatte fiir hohe Frequenzen 
(f 2 500 Hz) analog zu den Messungen des Absorp- 
tionsmaterials auf der kleinsten Raumbegrenzungs- 
flache (Kurve a, in Bild 13) schon eine Streuflache von 
71,5 m? geniigt, um keine durch mangelnde Diffusitat 
hervorgerufenen systematischen Fehler und die gering- 
sten zufalligen Fehler zu erzielen. Fiir diese Kontroll- 
messungen mit 42 m? Priifmaterial geben wir die Ab- 
sorptionswerte nur fiir die Frequenzen f > 500 Hz an, 
da bei den Frequenzen 160 Hz < f< 500 Hz eine Er-— 
hodhung des Absorptionsgrades bis zu 0,09 gegeniiber 
Gstat auftrat. Diese Erhéhung war durch ein Hohlliegen 
und Mitschwingen der Sillan-Platten bedingt, die zum 
Schutz ihrer Kanten und Ecken mit einem diinnen 
Blechstreifen eingefaft waren. Beide Fehlerursachen 
wurden durch Versuche an kleineren Priifflachen be- 
statigt und durch Auflegen von Massen verhindert. Da 
der Hallraumboden an einigen Stellen nicht ganz eben 
ist und mit drei Wanden keine scharfen Kanten bildet 
(Hohlkehlen), war es bei der Messung mit 42 m® Priif- 
material nicht méglich, diese Fehlerursachen zu ver- 
meiden. 
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an (vergleiche Bild 3, Kurven cy und d). Nach Kor- 
rektur der Druckerhohung (9.3.) ergibt sich der 
Verlauf a, beziehungsweise b,. Wird die Anderung 
der mittleren freien Weglange (9.2.1.) und die 
Lage des Prifmaterials im Hallraum (9.2.2.) rech- 
nerisch beriicksichtigt, so betragt die Differenz zwi- 
-schen dem ,,richtigen“ Absorptionsgrad a, und agtat 
weniger als + 0,025 (Kurve a, beziehungsweise bs). 
Die Abweichung des gemessenen und korrigierten 
Absorptionsgrades von dem aus der Impedanz 
berechneten liegt somit innerhalb der nach Bild 9 
za erwartenden Grenze der zufalligen Fehler. Dabei 
mu beachtet werden, das auch die Rohrmessung 
und damit a,,, mit einem zufalligen Fehler behaftet 
ist. Systematische Abweichungen zeigen sich nicht 
‘mehr. 


10. SchluBfolgerungen 


Die im Hallraum zur Verbesserung der Diffusitat 
des Schallfeldes benotigten Streukorper verringern 
die zufalligen Fehler bei Nachhallmessungen. Bei 
unseren Untersuchungen zeigten sich in ihrer Wir- 
kung gekrimmte Platten einem bewegten Fligel 
gleicher Streuflache gleichwertig beziehungsweise 
leicht tiberlegen, einem stehenden Fligel weit tiber- 
legen. Dabei ist es giinstig, Streukorper verschiede- 
ner Abmessungen zu kombinieren. Die Diffusoren 
diirfen nicht zu klein bemessen werden, damit auch 
bei tiefen Frequenzen ihre schallzerstreuende Wir- 
kung erhalten bleibt. Die giinstige Wirkung der 
Streukorper auf die zufalligen Meffehler ist be- 
grenzt, da mit der VergroBerung der Streuflache 
gleichzeitig eine VergroRerung der Leerraumabsorp- 
tionsflache verbunden ist, wodurch der Fehler wieder 
anwachst. 
In unserem 125m? groBen Hallraum benotigen 
wir im Leerraum fiir Frequenzen f = 500 Hz keine 
und fiir tiefe Frequenzen 27 m? Streuflache (Vorder- 
und Riickseite zusammen), das hei®t 16% der Hall- 
raumoberflache, um das Minimum der zufalligen 
Fehler zu erreichen. Bei einer Priifflache von 16 m? 
hochabsorbierenden Materials (das heif®t 133% der 
Norm-Prifflache beziehungsweise 10% der. Hall- 
raumoberflache) sind 43 m? Streuflache (das ent- 
spricht 25% der Hallraumoberflache) notig, um auch 
bei tiefen Frequenzen die kleinsten Abweichungen 
zu erhalten. Fiir Frequenzen f = 500 Hz hatte dafiir 
schon eine Streuflache von 36m? ausgereicht. Bei 
der angegebenen Grose der Streuflachen ist der 
durch mangelnde Diffusitét verursachte systemati- 
sche Fehler beseitigt. Die benotigte Streuflache gilt 
fiir die angegebene Absorptionsflache. Verkleinert 
_ beziehungsweise vergréBert man diese, so kann 
_ beziehungsweise mu man die Streuflache ebenfalls 
" variieren. 
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Im einzelnen ergeben sich mit den oben geforder- 
ten Streuflachen fiir den mittleren Fehler des Ab- 
sorptionsgrades folgende Werte: 5% fiir f= 125 Hz; 
2% fiir s = 500 Hz. 

Damit haben wir gezeigt, da} es in einem 125 m3 
groBen Hallraum moglich ist, Absorptionsgradmes- 
sungen auch bei tiefen Frequenzen mit guter Genau- 
igkeit durchzufthren. 

Die von Meyer und Kurrrurr im Modellraum 
fiir ein relativ sehr dickes Material bestimmte untere 
Grenzfrequenz, bei der der mittlere Fehler des Ab- 
sorptionsgrades 10% betragt, wurde bei dem von 
uns untersuchten Priifmaterial mit 5 cm Dicke nicht 
gefunden. 

Um systematische Fehler klein zu halten, ist es 
ratsam, das Priifmaterial weder auf der gro ten 
noch auf der kleinsten Hallraumbegrenzungsflache 
anzubringen. Im ubrigen ist es mit den in der vor- 
liegenden Arbeit angegebenen Korrekturen moglich, 
die durch die Hallraumdimensionen bedingten syste- 
matischen Abweichungen vollig auszuschalten. Bei 
Vergleichsmessungen in verschiedenen Hallraumen 
hat man nur noch mit den kleinen zufalligen Fehlern 


zu rechnen. 
(Eingegangen am 21. Juli 1960.) 
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Mitteilung 


Eine internationale Tagung mit dem Thema: ,,Der 
audiogene Krampf betrachtet vom psychophysiologi- 
schen, neuropharmakologischen und biochemischen Ge- 
sichtspunkt“ findet im November 1961 in Gif sur Yvette 
(Seine et Oise) statt und wird vom franzésischen Natio- 
nalen Zentrum fiir Wissenschaftliche Forschung 
(C.N.R.S.) abgehalten. 

Personen, die an diesem Symposium interessiert sind, 
wenden sich bitte an Herrn Prof. Dr. R. G. Busnet, Di- 
rektor des Laboratoriums fiir physiologische Akustik, 
(C.N.R.Z.), Jouy-en-Josas, (S. et O.), France. 


«Psychophysiologie, neuropharmacologie et biochimie 
de la crise audiogéne», un colloque international qui 
aura lieu du 6 au 9 Novembre 1961 a Gif sur Yvette 
(Seine et Oise), France, organisé par le Centre Natio- 
nal de Recherches Scientifiques frangais. 


“Psychophysiology, neuropharmacologie, and _bio- 
chemistry of the audiogenic seizure”, an international 
symposium held on November 6—9th 1961 at Gif sur 
Yvette (Seine et Oise), France, by the French National 
Center for Scientific Research. 
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